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Abstract——Investigations on water flow in Raschig ring packings were carried out employing 
columns constructed of wire gauze to minimize wall effects. A study of the mean flow through 
the ring interiors revealed that this decreased from 29 to 12 per cent of the total flow as the 
packing size increased from | in. to 1} in. Examination of cross-sections cut through a model of 
a ring packed bed suggested a cause of the observed variation. 


Résumé —Pour étudier Pécoulement de dans des garnissages en anneaux Raschig auteur 
utilise une colonne en toile métallique afin de réduire les effets de paroi. L’étude de écoulement 
moyen a travers lintérieur des anneaux montre que cet écoulement décroit de 29 a 12°, de 
lécoulement total quand la dimension du garnissage augmente de 4 in, a 1} in. L’examen de la 
section perpendiculaire dun modéle de lit de garnissage cn anneaux a suggéré une cause de la 
variation observée. 

Zusammenfassung — Die Wasserstromung in Packungen aus Raschigringen wurde untersucht, 
wobei zur Verminderung des Wandeffektes eine Kolonne aus Maschendraht verwendet wurde. 
Die Untersuchung der mittleren StrOmung durch das Innere der Ringe ergab, dass diese von 
29 auf 12°, der GesamtstrOmung abtiel, wenn der Ringdurchmesser von 12 auf 32 mm erhdht 
wurde. Eine Nach priifung von Querschnitten durch cin Modell einer Ringpackung deutete auf 


die Ursache der beobachteten Verinderung. 


INTRODUCTION 


In the past, several attempts have been made to 
determine the extent of surface wetting occurring 
when liquid flows over Raschig ring packings. 
The methods used fall into two distinct categories, 
namely direct and indirect. Of the direct methods, 
Mayo and others [1] used paper rings with dye 
solution as the wetting fluid, while SuuLMan and 
others [2] employed specially fabricated naphtha- 
lene rings, measuring the concentration of 
naphthalene in the issuing gas stream for wetted 
and for dry rings. The indirect methods basically 
consist in comparing mass transfer data for the 
system in question with data for a similar system, 
in which it was known that the whole of the ring 
surfaces was utilized. For example, SHerwoop 
and Picrorp [3] and Weisman and Boni [4] 


used the data of Tarcker and Hovcen [5], 
obtained for the evaporation of water from rings 
with a porous surface which ensured complete 
utilization. This was compared, by the use of 
Cuivron and CoLBuRN’s jp) mass transfer factor 
[6], to the mass transfer data of several workers, 
and the effective area thereby determined. 
SuuLMAN [2, 7] calculated the jp factor for the 
sublimation of naphthalene from dry rings and 
compared this with the mass transfer data of 
FeLuincer [8]. 

The results of these investigations are somewhat 
contradictory. Thus, Mayo found that the wetted 
area varied with liquid rate but was independent 
of gas rate, whereas WreisMAN and found 
substantially the reverse. Making this one 
exception, however, most workers are in agreement 
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with Mayo’s conclusion. Similarly, Jesser and 
Excin [9], investigating the related phenomenon 
of liquid hold-up, were able to show that this also 
was independent of gas flow. Considerable 
disagreement is, nevertheless, apparent in regard 
to the effect of such factors as ring size and liquid 
rate. In such a complex system as an irrigated, 
packed bed many causes can be suggested for 
such disagreement. For instance, Bonp [10] 
showed that the absorption of a smal! quantity 
of ammonia was suflicient to rupture a hitherto 
stable water film. Furthermore, Baker et al. 
[11] showed that unless the ratio of column 
diameter to packing diameter exceeded eight, a 
high percentage of the liquid would flow adjacent 
to the wall. Many of the wetted area determina- 
tions were carried out under conditions not in 
accordance with this ratio. 

In view of the complexity of these systems, it 
was considered that it might be profitable to 
attempt to resolve them into two distinct and 
perhaps thereby less complex systems. Liquids 
flowing through a ring packed bed must, at any 
one point, be flowing over the exteriors or through 
the interiors of the rings. In flowing over the 
exteriors there is no constraint on the direction 
taken by the liquid, beyond that its resultant 
should be vertically downwards. Through the 
ring interiors, however, the flow is constrained to 
a stream running parallel to the ring axis. No 
previous work having been reported along these 
lines, an what 
proportion of the liquid flows through the ring 
interiors. 


investigation was made into 


APPARATUS AND EXPERIMENTAL 


The four ring sizes used in the work were of 
}in., fin., Lin. and 1}in. diameter. Two columns 
were employed, each of length 3ft, one having 
a diameter of 3 in., for use with the } and 1 in. 
rings, and one of 4 in. for use with the 1 and 1} in. 
rings. Since Baker's limitation of 8:1 tower to 
packing diameter was not observed with these 
sizes, it was decided to minimize wal! effects by 
minimizing wall area. In accordance with this, 
both columns were constructed of } in. copper 
mesh. Initially it was feared that this mode of 
construction would result in liquid spurting out 
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of the sides of the column. In the event, this fear 
proved to be groundless, an occasional drop of 
liquid only being ejected by this means. 

The distributors consisted of cylindrical shells, 
arranged to be push-fits into the tops of the 
columns. At right angles to the axes of the 
evlinders, distributor plates were fitted, drilled 
with concentric circles of } in. holes to secure an 
even distribution. To keep the impingement 
velocity of the liquid passing into the column 
to a minimum, small domes were mounted just 
above the plates, so that the liquid, leaving the 
circulatory system nozzle (q.v.), struck the dome 
and was projected horizontally against the vertical 
shell wall, where its kinetic energy was sub- 
stantially destroyed. 

The columns were suspended above a reservoir 
which collected the liquid leaving the column base. 
A small pump took suction from the base of this 
reservoir and discharged, via an orifice plate, 
through a nozzle which projected the liquid 
vertically down onto the distributor dome. 

To determine the liquid flow through the 
interiors of individual rings, four large apertures, 
1} im. square, were cut at varying levels in the 
column sides. The ring coinciding with the 
column centre line and the centre of one of the 
apertures was then removed from the column 
after removing those rings immediately beneath 
it, care being taken not to disturb the surrounding 
rings. A specially tapped ring of the same size 
was then inserted in the place of the originally 
selected ring and the other rings replaced. The 
tapping consisted of a toy balloon, cut across the 
diameter parallel to the mouthpiece, and stretched 
over one of the mouths of the ring. The mouth- 
piece led, via a length of rubber tube, to a measur- 
ing cylinder, A plate was provided around the 
tube to prevent liquid from reaching the cylinder 
via the tube exterior. The only limitation in 
replacing this tapped ring was that the tapped 
mouth should be lower than the free mouth. 
Fig. 1 illustrates the tapped ring in_ position 
in the column. 

Interior flows were determined (by means of 
the measuring cylinder and a stop-clock), for 
varying values of the overall flow rate for a series 
of individual rings at cach ring size. The number 
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Fic. 1. Tapped ring assembly in position. 
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of rings in each series investigated in this manner 
ranged from twenty to forty. During the indivi- 
dual determinations, it was noticed that normally 
the flow through the ring interiors was smooth 
and uniform, becoming rather more erratic at 
very low interior flows. This was not considered 
to be of any great importance, since its effect on 
the average value would be negligible. 

Water was used as the circulating fluid through- 
out. 

ReEsULTS 

The interior flows at overall flows of 0-5, 1-0, 
1-5, 2-0 and 2-5 ml /sec cm*® of bed cross-section 
were obtained for each ring size by interpolating 
between the experimental values of overall 
flow. These interpolated values were averaged to 
obtain mean interior flows in ml sec. The results 
are summarized in Table 1. 

The forty individual determinations for the 
jin. rings were analysed statistically but no 
significant difference could be detected between 
the results of each of the four apertures. Were 
any wall effects present, a falling off in the interior 
flow would be expected from the central rings as 
the liquid descended. It was therefore concluded 
that the method of column construction adopted 
fulfilled its purpose of minimizing wall effects. 

In order to compare interior flow with overall 
flow, it is first necessary to convert the former 
from ml/sec into ml/cm*sec units. For this 
purpose the packed bed can be considered as a 
series of layers of rings, the rings having a definite 
average cross-section of this layer per ring. From 
the geometrical similarity of Raschig rings of all 


Table 1. 


sizes it can be shown that the ring diameter should 
be inversely proportional to the cube root of the 
number of rings per unit volume. This was 
verified for the rings used in the present work, 
the constant of proportionality being such that 
the average volume of bed taken up by a ring 
was equivalent to a cube of side 6-9 per cent in 
excess of ring diameter. The average area of 
cross-section taken up by a ring would be the 
square of this side. This square was therefore 
divided into the results of Table 1 to obtain 
mean interior flows in ml cm? sec, these being 
plotted in Fig. 2. 


Cross-secTions THroucnu Packep Beps 


Figure 2 indicates the variation of interior 
flow with ring size and overall flow. To obtain 
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Fic. 2. Variation of interior flow with overall flow. 
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further insight into this variation, a model of a 
packed bed was constructed with the object of 
ascertaining the appearance of plane cross- 
sections through these beds. Solid } in. cylinders 
were cut from } in. thick slabs of black plasticine 
by means of a cork borer. These cylinders were 
packed at random into an open topped tin-plate 
box, 12in. high with a six inch square base. 
Molten paraffin wax was poured into the box 
which had just previously been heated in an 
oven to a temperature slightly in excess of the 
wax melting point in order to preclude any 
premature wax solidification. After thorough 
cooling, the tinplate was stripped off and the 
block thus obtained sliced in planes parallel to 
its base by means of a guided wire. The resultant 
cross-sections were then traced, a typical tracing 
being reproduced in Fig. 3. 


Examination of these cross-sections reveals, as 
would be expected, that all the profiles appear as 
conic sections of various types. Straight lines 
must represent sections thorough the ring mouths 
or sections through the cylindrical surfaces of 
rings whose axes are parallel to the plane of 
cross-section. If liquid is considered to be flowing 
through one of these cross-sections it is apparent 
that the interior flowing liquid must be moving 
perpendicularly to the mouth; in half of the 
cases leaving the ring, while in the remaining 
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half, entering. It was postulated that the total 
length of the ring mouth profiles, when halved, 
corrected for wall thickness and divided by the 
total cross-sectional area, represented a measure 
of a bed’s capacity to accept liquid into the ring 
interiors, This measure was accordingly denoted 
as the “inlet availability’ of the bed, and 
on being computed for the mode sections 
(? in. rings) was found to have a value of 0-101 
em/em*?. From the geometrical similarity of 
Raschig rings, it follows that inlet availability is 
inversely proportional to ring size and hence can 
be calculated for all the ring sizes used in the 
present work. 

At given values of the overall flow rate the 
interior flow rate is directly proportional to the 
inlet availability, values of the proportionality 
constant being plotted against overall flow in 
Fig. 4. This plot is linear; its failure to pass 
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Fic. 4. Variation of the ratio of interior flow to inlet 
availability with overall flow. 


through the origin possibly indicating that a 
threshold value of overall flow must be reached 
before interior flow commences. With the aid of 
this plot, and a knowledge of the inlet availability, 
interior flows can be estimated. 


CONCLUSIONS 


The interior flow in ring packed beds has been 
shown to vary linearly with overall flow through 
the bed. Figure 2 demonstrates that as ring size 
increases there is a definite falling off in interior 
flow. (This constitutes about 29 per cent of the 
total flow through the bed for } in. rings, 19 per 
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cent for ? in., 14 per cent for 1 in. and 12 per cent 
for 1}in. rings). With ring sizes in excess of 
1}in it can be anticipated that this proportion 
will be 10 per cent or less and that, as a conse- 
quence, the interior surface can contribute little 
to the mass transfer in the bed. A model of a 
packed bed has been used to explain the observed 
variations and to derived a method of predicting 
the interior flow. 

Further work demonstrates how a knowledge of 


the interior flow can be used to resolve liquid 
hold-up into interior and exterior components, 
thereby obtaining better insight into the flow 
through ring packed beds. 
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Abstract—The progress of a reaction catalysed on a tube wall is analysed for an assigned radial 
pattern of axial velocity. A procedure is given for mathematically transforming an observed axial 
profile of cross-section average reactant concentrations into a plot of the instantaneous reaction 
rate vs, the wall concentration at which the reaction is proceeding. It is also shown how to calculate 
the average conversion at any contact time for an arbitrary reaction rate function, 


Résumé—L‘évolution dune réaction catalysée sur la paroi d'un tube est étudiée pour une 
distribution radiale fixe de la vitesse longitudinale, Un moyen est donné pour transformer 
mathématiquement un profil observe longitudinal de concentration moyenne de réactants pour 
une section déterminée, en fonction de la concentration de paroi a laquelle la réaction se développe. 
I/auteur montre aussi comment calculer la conversion moyenne pour un temps de contact 
queleonque pour toute fonction arbitraire de la vitesse réactionnelle. 


Zusammenfassung Das Fortschreiten ciner Reaktion, die an einer Rohrwand katalysiert ist, 
wird fiir eine angenommene radiale Verteilung der Axialgeschwindigkeit untersucht. Ein Verfahren 
wird angegeben zur mathematischen Tranformation eines beobachteten axialen Profils der tiber 
den Querschnitt gemittelten Konzentrationen der Reaktanten in eine Beziehung zwischen der 
momentanen Reaktionsgeschwindigkeit und der Wandkonzentration, bei der die Reaktion 
fortschreitet. Es wird auch gezeigt, wie man die mittlere Umsetzung bei irgend einer Kontaktzeit 
berechnen kann fiir eine beliebige Funktion der Reaktionsgeschwindigkeit . 


1, INTRODUCTION analyses, although certain new forms of answer to 


in some reactive material, flows 
down the inside of a tube, At the tube wall, the 
reactive molecules are converted into a new 
substance at a rate which is some function of their 
concentration at the wall. The reactive molecules 
diffuse out to the tube wall just rapidly enough 
to balance this reaction rate. 

Two questions suggest themselves. The first 
is a question of kinetic analysis: without making 
any prior assumption about the form of the rate 
function, what accessible measurements can one 
transform mathematically into a knowledge of the 
instantaneous reaction rates and of the corres- 
ponding wall concentrations of the reactive 
subtance. The second is a design question: 
knowing the rate function and the other physical 
parameters of the system, how can one calculate 
the contact time, say, needed to achieve a desired 
conversion. The present study is aimed at kinetic 


A Gas, rich 


design questions will appear as a by-product. 

The design question has been considered in a 
variety of contexts. For example, TowNsenp 
[1] studied a gas flow with a parabolic velocity 
profile, and solved an ion-diffusion problem 
with an absorbing tube wall. This corresponds 
chemically to the case of infinitely fast reaction. 
Tayior [13, 14] studied a similar flow, and gave 
an approximate solution to a diffusion problem 
with a reflecting tube wall. This corresponds 
chemically to the case of infinitely slow reaction. 
DamkOuLer [2] and Baron and his co-workers 
[3] solved corresponding problems with first 
order reactions of finite speed taking place at the 
tube wall, but with a flat velocity profile (plug 
flow) in the tube. 


Cuampre [15] formulated this class of problems 
ina general way that accounted for the hydrodyna- 
mics of the boundary layer which develops at the 
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wall, as well as for the reaction catalysed on the 
wall. He gave, for certain classes of flows, integral 
equations expressing surface reaction rate as an 
integral operator on surface reactant concentra- 
tion. Also, he and Acrivos [4] gave an exact 
(in terms of boundary layer theory) analytical 
solution of their enlarged design problem for a 
limited class of rate functions. Further, by 
approximating the gas velocity profile near the 
catalytic surface, they arrived at integral equations 
of the same kind for the general boundary layer 
situation, expressing surface reaction rate as an 
integral operator on wall reactant concentration. 
These integral equations furnished in principle a 
method of carrying out kinetic analyses, in that 
if a reactant concentration profile were measured 
at the catalytic surface, it could be transformed 
into a surface reaction rate by applying a suitable 
integral operator, 

The present study sets itself the more limited 
objective of accounting for a reaction catalysed 
on the wall of a cylindrical tube, where the process 
gas undergoes an axial flow with a preassigned 
radial gradient. In this limited context, it is able 
to indicate a much easier experimental way to 
make kinetic analyses, by permitting the experi- 
menter to work with cross-section average con- 
centrations, rather than forcing him to measure 
concentrations right at the wall. A solution to the 
design problem is given, in terms of an integral 
equation in the wall concentration of the reactant. 
The integral equation gives the wall reactant 
concentration as an integral operator on the 
surface reaction rate. This inversion of the integral 
equations of [4] and [15] permits the direct calcula- 
tion of reactant concentration at the wall from a 
knowledge of the instantaneous surface reaction 
rate. It is shown further that if the cross-section 
average reactant concentration be known as a 
function of contact time, a suitably scaled 
differentiation with respect to contact time will 
give a profile of instantaneous reaction rates. 
Thus, without making any prior assumption about 
the form of the rate function, one can obtain 
reaction rate as a function of wall concentration 
by making suitable measurements of the cross- 
section average concentrations alone. 

The numerical procedures involved in setting 


up the integral equations needed to handle either 
the design problem or the kinetic analysis problem 
are quite formidable. They can, however, be scaled 
so as to depend only on the shape of the velocity 
profile, and not at all on the actual speed of the 
gas, the diameter of the tube, the size of the 
diffusion coetlicient, the inlet condition of the gas, 
or the form of the rate function. These preliminary 
set-up calculations can thus be performed once 
and for all for any physical situation involving, 
say, a parabolic velocity profile. The calculations 
that then have to be made in any particular case 
of design or data analysis are, in comparison, 
quite modest in size. 

For our present purposes, it seems to be 
necessary to evaluate the kernel in our integral 
equation anew, despite the wealth of mathe- 
matically very similar problems already treated 
in the literature (see, for example, [17]). Perhaps 
the closest to our needs is the numerical solution 
to the Graetz problem for a parabolic velocity 
profile given by Seviars and his co-workers [16]. 
Their solution, however, is given as an expansion 
in eigenfunctions vanishing at the tube wall, and 
so, because of Gibbs’ phenomenon dilliculties 
(see, for example, [5]) cannot, as we require, be 
readily evaluated at the wall. It has, as a result, 
limited utility for our kinetic analyses, although 
it can be used as a kernel in the integral equations 
of [15] for design calculations. 

The present mathematical analysis holds for 
quite general velocity profiles. Since, however, 
axial diffusion is neglected compared to the main 
convective effect of the gas stream, the range of 
application will likely be for Poiseuille flow 
(parabolic profile) and a correspondingly small 
molecular (laminar) diffusivity. It should be 
noted also that the mathematical analysis is all 
for a steady state, isothermal situation. An 
experimental study that meets these conditions 
is in progress, and will be reported in the near 
future, along with the explicit numerical schemes 
for analysing the data. 

Section 2 below formulates the mathematical 
problem corresponding to the physical system in 
question, and Section 3 carries out the reduction 
to an integral equation. Section 4 sketches the 
use of the integral equation in kinetic analyses, 
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and Section 5 its use in the design problem. 
Section 6 is essentially a mathematical appendix 
on the development in a suitable eigenfunction 
expansion of the kernel of the integral equation 
given in Section 3. 


2 Tur Dirrerentiat Equations 


In a cylindrical tube, establish a co-ordinate 

system : 
z — axial co-ordinate (measured in the direc- 
tion of gas flow from = — 0 at inlet), ft 


r = radial co-ordinate (measured from r = 0 
at tube centre to r — R at tube wall), ft 


Gas flows down the tube at a steady velocity 
U = average axial velocity (in the direction of 
increasing axial co-ordinate 2), ft/sec 
and with a velocity profile p(r), so that 
Up(r) = axial velocity at radial co-ordinate r, ft/sec 


The velocity profile p must then be taken to be 
normalized on the cross-section area of the tube 


(1) 


R 
| Pind 1. 


A general power law class of profiles p vanishing 
at the tube wall is 


where » = 1 corresponds to the Poiseuille para- 
bolic flow, »y = 2 to the turbulent fourth power 
profile, and infinite » to plug flow. The flow is 
taken to be at constant temperature. 

In the process stream, let the concentration 
of the reactive molecules be 


e(r, 2) = reactant concentration, mols ft* 

with the uniform value at the tube inlet 

€ = inlet reactant concentration, mols, ft® 
Let the reactive molecules have diffusivity 

D = radial diffusivity of reactant, ft®/ sec 
(axial diffusivity is neglected), and let them be 
consumed at the tube wall at the rate 

k (c) = surface reaction rate, mols, ft® sec 


Katz 


It is this function k whose form will be left 
unspecified in what follows. The assumption of 
constant temperature permits us to take & to 
depend only on the reactant concentration ¢. 

A differential material balance on the reactive 
molecules at any point in the tube then gives 


= = @; 
rw Ptr) 
z>o0 


At the tube inlet, the reactant concentration is 
specitied 
C=; O<r< R, 


(4) 
At the tube wall, a differential material balance 


R, >0 (5) 
Finally, for mathematical completeness, a regu- 
larity condition at the tube centre is imposed 


¢ bounded; r= 0, (6) 


Equations (3)-(6) represent the mathematical 
specification of the concentration function ¢ (r, z). 

Of particular interest in what follows will be 
the cross-section average reactant concentration 
weighted by the local gas velocity, 


K 
é (z) | 2) (7) 


This is a true average, in view of the normalization 
(1) of p. The definition (7) of @ imagines a 
measurement as made by collecting the total 
flow past a cross-section of the tube for a given 
time, and assaying the material so collected for 
reactant. It is adopted as conforming most 
closely to the experimental work mentioned in 
the ing section, but it is by no means the 
only reasonable definition. For example, if one 
viewed a measurement as made by focussing 
attention on a differential slice of the tube and 
instantaneously assaying the material in that 
slice for reactant, the simple area average 
K 
1 


= 


e(r, 2) 2ardr 


would be the more appropriate. 
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It will be more convenient to carry out the 
subsequent analysis in terms of certain dimension- 


less variables. Introduce the dimensionless co- 
ordinates 
r 


-~ dimensionless radial co-ordinate 
_ 

== dimensionless axial co-ordinate 


(8) 


It will be noted that @ is a contact time, normalized 
on a diffusion time. Define in terms of the 
velocity profile p the weighting function 

= 2ap(Re) 


with, from (1), 


[ = 


(9) 


The power law form of p (2) corresponding to (2) 


Is 
(v + 1) 


= poh —z”); 


v(l (10) 


where, as before, vy — 1 corresponds to parabolic 
flow, ete. Finally, define a dimensionless concen- 
tration function 


J 


2D 
-— dimensionless reactant concentration 


and a dimensionless rate function 
R 
h(v) = cod) (eo?) (12) 


— dimensionless surface reaction rate 


In terms of these new quantities, the boundary 
value problem in equations (3)-(6) becomes 


d dv dv 7 
v= 1 
> (13) 

v bounded 
av 

h(v) = 0 
(v) 


Chemical reactions catalysed on a tube wall 


205 


The equations (13) represent the mathematical 
specification of the dimensionless concentration 
function v(a,@). The parameters R, cy, D, U 
have been sealed out, and all that remains is the 
velocity profile shape p (x), and the dimensionless 
rate function h(v). Detine finally an average v, 


v (0) (7, 0) p(x) dx 


(14) 


By virtue of the normalization (9) of p, this is a 
true average. Also, 


v (A) | 


(15) 


RU 
) 


by (7) and (11). 
The following section will discuss the reduction 
of (13) to an integral equation. 


3 REDUCTION TO AN INTEGRAL 
EQUATION 


The reduction of (13) to an integral equation 
begins by replacing the assigned function of 
concentration h(v) by an assigned function of 
“time” {(@). The non-linear boundary value 
problem (13) is thus replaced by the linear 
problem 


ow 
— = p(t) — 0< 
awe \ de) 

s@<2<10=80 


(16) 


zw bounded 


+f(@)—o 


in the new function w (a2, @). There is a certain 
arbitrariness in the choice of the problem (16). 
One could, for example, replace the assignment 
of dw Qw at # = 1 by an assignment of w. This 
would lead however in Section 6 to an eigenvalue 
problem (28) with ¢ = 0 at # = 1, and the series 
expansions made in these eigenfunctions would be 
subject to Gibbs’ phenomenon difficulties (see, 
for example, [5]) at the wall x = 1. 


Ss. 


Introducing Laplace transforms (see, for 
example, [6]) suggests the form of the solution 
of (16). Setting 

7 
w* (2, = w (x, 0) dé 


f* (0) = | (8) 


the equations (16) give 

d 

—j}> — 1}; O 
w* bounded 


=0 


The equations (17) can be solved in the form 
w* 0) = I* (a, 2) f* (18) 


where I* has the nature of a transfer function 
connecting f with w and must itself satisfy the 
equations 


> (19) 


Now I (x, @), the inverse Laplace transform of 
I* (x, 
(a, a) = ic °° I (x, 0) de 


is essentially the response of w to a unit impulse 
f, since (18) can be translated back into the time 
domain to give the solution of (16) in terms of 
as 


6 


Taking the Laplace transform inverse of (19) 
gives a boundary value problem for J 


Katz 


1-0 
bounded =-0,0>0 
_ 


(21) 


In (21), 8(@) is a Dirac delta function (see, for 
example, [7]), strictly, a right-sided delta function 
with 


8 (7) — | 
@>o0 


0; @=0 


for suitable testing functions . 


The question of solving (21) will be taken up 
in the concluding section of this study. It should 
be noted now, however, that J depends only on 
on the velocity profile shape p, and not at all on 
any other specific physical or geometrical charac- 
teristic of the physical system. Proceeding 
then as though 7 were known, recall from (13) 
and (16) that w reduces to the dimensionless 
concentration v if f(@) is replaced by h {v (1, 0)}. 
Making this identification in (20) gives 

6 
A) = 1 - — r)hf{v(l, r)}dr (22) 


Taking « — 1 in (22) gives the desired integral 
equation in the tube wall concentration v (1, @) 
v(1,0)=1 | we r)h{o(l, r)}dr (23) 
where 


M (0) =1(1, @) (24) 


Equation (23) contains the root of the solution 
of the boundary value problem (13). In kinetic 
analyses, it will be used to translate a knowledge 
of the reaction rate h{v (1, @)} to a knowledge 
of the corresponding wall concentration v (1, @). 
In making designs, (23) will be solved for v (1, @), 
and correspondingly h {v (1, @)}, and the solution 
substituted into (22) to give the detailed per- 
formance of the reactor. These matters will be 
discussed in the two sections following. 
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4. AppLicatTion to Kinetics Srupies 


Preliminary to the application of the integral 
relation (23) to rate data lies the determination 
of these rate data from cross-section average 
concentrations. 

To see how this determination is to be made, 
integrate the first equation of (13) with respect 
to z. Recalling (14) gives 


1 
dv(@) dv (a, A) 
[ow dx 
1 


which, on applying the boundary conditions in 

(13), becomes 

(0) 


hf{v (1, 


(25) 


A kinetics experiment might then proceed 
along the following lines. A suitably fine grained 
profile of cross-section average reactant concen- 
tration @ (z), as defined in (7), would be measured. 
(More realistically, measurements of @ might 
be made on a single apparatus for a single value 
of z over a suitably fine grained range of the 
average gas velocity U), The measured @ would 
then be scaled, as in (8), (15), to a profile of 
dimensionless average concentrations (@). 
Differentiating 0(@) numerically according to 
(25) would give the corresponding profile of 
dimensionless reaction rate h {v (1, @)}, and carry- 
ing out the integration in (23) numerically with a 
precalculated M (@) would give the corresponding 
profile of dimensionless wall concentrations 
v (1, 0). These dimensionless wall concentrations 
and reaction rates could then be scaled back to 
the physical quantities ¢(R,2) and kfe(R,2)} 
according to (8), (11), (12). This sealing back 
and forth would, of course, require a knowledge 
of the gross physical parameters of the system: 
the inlet reactant concentration co, the reactant 
diffusivity D, the mean gas velocity U, the tube 
radius R. 


The end product of the analysis would be a 
table (plot) of the instantaneous reaction rate 
k(c) against the reactant concentration ¢ for 
whatever range of wall concentrations ¢ (R, 2) 
happened to be implicit in the experimental set-up. 
Such information could then serve to test assump- 
tions about the reaction mechanism, exactly as 
if it had been obtained from, say, suitable measure- 
ments made right at the catalytic surface of a 
differential reactor. 

The same sort of analysis could of course be 
made for homogeneous reactions, but there, as is 
clear (and well known), measurement of average 
concentration will in general give knowledge 
only of average reaction rates, It is just the 
specific character of the surface catalysed reaction 
that permits measurements of average concentra- 
tion to lead to knowledge of the reaction rates 
at specific places where the working concentrations 


can be deduced as well. 


5. APPLICATION Reactor DeEsIGn 


Suppose that here the surface rate function k (ce) 
and the velocity profile Up (r) have been specified, 
as well as the other parameters of the physical 
system: the inlet reactant concentration ¢9, the 
reactant diffusivity D and the tube radius R. 
Suppose further that the dimensional scaling in 
(8), (11), (12) has been carried out, leaving only 
the dimensionless reactant concentration v (x, @), 
the dimensionless rate function h (v), the velocity 
profile shape p (wv), and of course the precalculated 
integral kernels, J (2, 0) and M (@). 

The first step in the design would be to solve 
the integral equation (23). This would give the 
wall reactant concentration v(1,@), and con- 
currently the instantaneous surface reaction rate 
h{v(1,0)}. Then, if one wanted to predict the 
detailed performance of the reactor, one could 
enter (22) with h {v (1, @)} as a known function of 
@, and calculate v (x, @) for any desired x and @. 
If, however, it were suflicient to know the average 
conversion at any cross-section, one would 
proceed differently. Averaging the initial condi- 
tion of (13) according to (14) gives 


v(0)=1 


and with this, (25) can be integrated to give 
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dz 


v(@)=1 t)} dr 


(26) 


In (26), the integrand can be regarded as a known 
function of @, following the solution of (23), and 
one running integration would give a plot of 0 
against @ from which any desired conversion 
could be picked off. 

It remains to discuss the solution of the integral 
equation (23). This is a non-linear integral 
equation of Volterra type, and can be solved, 
formally at least, by the iterative scheme 


6 
(1,0) = 1 [a (0 — r)hfv, (1, dr; 
(27) 


== @,1,2,... 
starting from any first guess at v, say, 
Uo (1, 0) = 1 


Convergence of this scheme is discussed for linear 
integral equations in, for example, [8]. The linear 
case, where h (v) is proportional to v, corresponds 
physically to the case of a first order reaction. 
Convergence is proved in the linear case under 
certain regularity conditions on the kernel M, 
encompassing the weak singularity to be expected 
in M (see the example in Section 6) because of the 
wall condition in (21). The extension of the 
proof to the non-linear case at hand is straight- 
forward, and requires only certain mild regularity 
conditions on the rate function h (v). For instance, 
that h, as a function of v, have a continuous 
second derivative would be enough not only to 
guarantee the convergence analytically, but to 
permit an integration by parts in (23) to remove 
the singularity of M. Such an integration by 
parts would furnish a means of carrying out the 
iterative solution (27) by repeated numerical 
integration, using a precalculated M (@) (or rather, 
its indefinite integral). The final integration in 
(26) would of course also be carried out numerically 

The design method sketched above (perhaps 
with an improved first guess at v (1, @) made along 
the lines of [18]) assumes that one has _pre- 
calculated the kernel M(@) in the integral 
equation (23). This, as discussed in Section 6, 
is itself a very formidable calculation, and would 
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likely be made only in cases where extensive 
kinetic analyses were to be carried out along 
the lines of Section 4. If an isolated design 
calculation had to be made, one would use instead 
the integral equation of [15], for which, in the 
case of Poiseuille flow, the kernel has been 
calculated in [16]. Or possibly, it would be more 
economical to work directly from the partial 
differential equations (3)-(6), or their dimensionless 
form (13). One would, as is commonly done, 
introduce a network of points in the 2-0 plane, 
replace the differential equations by the appropriate 
difference equations, and solve these difference 
equations numerically (see, for example, [11]). 


6 EVALUATION OF THE KERNEL 


This section will present expansions of the 
integral kernels J and M of Section 3 in series 
of certain eigenfunctions. The relevant ecigen- 
value problem arising out of (21) is to find 
constants A and functions ¢ (7) for which 


d dd\ 


bounded 


Problems very similar to this (with ¢ rather than 
dé dx vanishing at x = 1) were solved numerically 
by Townsenp [1] and by Seviars and his co- 
workers [16] for a parabolic velocity profile, 
v = Lin the p (2) of (10). These workers made their 
calculations in terms of analytic approximations 
to the solutions ¢(7). It is possible that more 
economical computing schemes can be arrived 
at by introducing, as is commonly done, a chain 
of points on the 2-axis, replacing (28) by the 
appropriate difference equations, and carrying 
out the calculation of the eigensolutions by what 
are essentially algebraic means (see, for example, 
[12}]). We shall, however, reserve for future 
report the question of numerical solutions of (28), 
and shall give here only an expansion, in series 
of these solutions, of the impulse response I (2, @) 
of (21) and its boundary value M (@) of (24). 

The problem (28) is an example of an eigen- 
value problem of Sturm-Liouville type with a 
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singular endpoint, such as is discussed for example 
in [9]. The singularity (the vanishing of the 
coefficient of d?¢/dz? at = 0) is the same as 
occurs in Bessel’s equation, which is indeed the 
particular case of (28) with p (a) = 2a, correspond- 
ing to infinite v (plug flow) in (10). There is only 
a sequence of values of A for which (28) has 
solutions ¢(#) not identically 0. In the cases 
of physical interest where p (x) > 0, these eigen- 
values A are all > 6, and may accordingly be 
arranged in the order 
(29) 
To each A,, corresponds a single eigenfunction 


= $, (2); (20) 


satisfying (28), and quite arbitrary functions f (7) 
can be expanded in series of the ¢,. Further, 
the eigenfunctions 4, (@) are orthogonal with the 


weight p (2) 


= @, 1, 2,...- 


and can clearly be normalized to make 


! 
0 


Thus a Fourier expansion 


(38) 


has the coefficients 
1 
= [ £02) (2) p(x) dn 


Note finally that the leading eigen solution (29), 
(30) is 


(34) 


= 0. dy (a) = 1 (35) 


as may be seen by inspecting (28). 

In solving (21) for the impulse response J, 
it turns out to be convenient to work through 
the step response 


s(a, 0) = t)dr 


Inserting the Laplace transform of s into (19) 
and inverting the result shows how to integrate 
(21) to obtain a boundary value problem for s 


8 
ar 


or 


s=0 


s bounded r= 0, 


ds 


1 = I, > 0 
ow 


from which J may be recovered as 


8 (2, 0) 


I (x, @) = 


(37) 
We rewrite the boundary value problem (36) by 
setting 

(38) 


8 (a, @) = 8, (7, + 8, (a, @) 


7 
p(t); 
0 


;0<2<1,0=0 
;%=—0,0>0 


°@<2<1,0>0 
;@<2<1,0=80 


8, bounded 6, 


ds 
or 


The term 2? 2 in (38) is taken to cancel the inhomo- 
geneous boundary condition in (36). It is taken 
as x? /2 rather than as, for example, x, in order to 
preserve the expected flatness of the solution at 
r= 0. 

The boundary value problem (39) may be 
solved in terms of the eigen solutions (29), (30) 
by direct separation of variables. Let 


(36) 
OL, 
10 
with 
oom « 
or Q 
n=0 d d = 
= 
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1 
a = = 6, 1,1 
0 
so that, recalling (33), (34), we have 
= a, 4, (2); 


n=0 
Then 
s, (z, 0) = — a,, d, (7) 


which, taking specific note of the leading eigen- 
solution (35), may be written 


8, (z, 0) = — a, — a,, 
nel 


The problem (40) is solved somewhat less 
directly. Let 
by = [ dy (0) de o=w@1,8.... 


so that, by (33) and (34), 


( 43) 


and note that, by (35), 
= 1 


Then assuming a solution for s, in the form 


8, (2,0) 4, (x) (0) (44) 


leads to the following ordinary differential 


equations in the 


dy, 
dé +A, = 5,3 


0 ; 


which have the solution 


Substituting (45) into (44) and recalling (35) gives 


i (46) 


8, 0) = 0 4 b, +, (2) 
nel 


Substituting (42) and (46) back into (38) and 
applying (37) gives the solution of (21) as 


I(x,0)=1 + (Aa, (47) 
nel 
and its boundary value M of (24) as 
M(@)—1 + 5 aa, + b,) (48) 
nel 


where a, and b, are given by (41), (48). 

The analytical behaviour to be expected of 
(47), (48) can be seen by considering the special 
case of the problem corresponding to plug flow 

p(z) = 2r 


obtained by taking infinite v in (10), The eigen- 
value problem (28) then has the Bessel function 
solutions 


by = 1 


(Xm 2). 
#0 (2): (%) 


ortho-normalized as in (31) and (32), with the eigen- 
values 


4 = 0 

A, an =1,2,... 

where the «, are the positive zeros of J, 
J, (%,) =O; 


O< ag<... 


| 


as tabulated, for example, in [10]. The kernels 
IT and M turn out to have the simple forms 


(a, 2) 
M(0)=14+ 
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The shrinkage forces developed in porous substances by capillary effects 
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Abstract Shrinkage forces in porous solids are of importance in filtration, drying and granulation. 
The geometrical method of calculated shrinkage forces is limited to regularly packed beds of 
spheres with small amounts of moisture, whereas the study of energy changes has provided an 
approach capable of dealing with porous beds of widely differing configurations and moisture 
contents, 


Résumé—Les forces de contraction dans les solides poreux sont importantes en filtration, 
séchage et granulation. La méthode géomeétrique pour le caleul de ces forces est limitée aux lits 
réguliers de garnissage avec sphéres et faibles humidités, L’étude des variations énergétiques a 
permis dobtenir une méthode de calcul des lits poreux de configurations et teneur en eau trés 
différentes. 


Zusammenfassung Schrumpfungskrifte in pordsen Feststoffen sind bei der Filtration, beim 
Trocknen und bei der Granulation von Bedeutung. Die geometrische Methode zur Berechnung 
von Schrumpfungskriiften ist auf regelmiissig geschiittete Kugelbetten mit kleinen Feuchtigkeiten 
begrenzt, wihrend die Untersuchung der Energiciinderungen auch fiir porése Betten von sehr 
verschiedenen Formen und Feuchtigkeitsgehalten anwendbar ist. 


1. INTRODUCTION 1 1 


P = A. + B.T. sin 
where A is the wetted area of each plate 
B is the wetted perimeter 
T is the surface tension of the water 
z is the angle of contact of the liquid 
ry is the smaller radius of curvature of the 
film and 
R, is the bigger radius of curvature. 


(1) 


Wuewn a wet porous substance is dried the drying 
process creates within that substance forces 
which tend to make it shrink. The amount of 
shrinkage exhibited on drying depends on the 
magnitude of the shrinkage forces and the com- 
pressibility of the substance. Shrinkage forces 
arise from the surface tension effects associated 
with the numerous, highly curved liquid-gas 
interfaces in the pores. Another more studied 
manifestation of these surface tension effects 
is the suction developed by a wetted porous 
substance and it is the purpose of this paper to 
examine the relationship between the suction 
and the shrinkage forces. 


Thus the expression for P consists of two terms. 

The first, A.T. (* A. x) describes the pull due 
OM 

to suction in the water and will be called the 

“suction” term. The second, B.T. sin x, 

describes the pull due to the tension in the 


water surface and will be called the “ surface 
tension " term. In the special case of flat parallel 


2. Tue Putt Betrween Two 
Dur tro tue Presence or Water 


Betrween 


Before proceeding with a discussion of the 
previous work done on shrinkage forces it is 
instructive to examine the pull developed between 
two plates due to the presence of a water droplet. 
The magnitude of the pull, P. is given by the 
following equation [1}. 


plates the “ surface tension "* term will disappear 
when the angle of contact is zero. If the plates 
are curved then this term will reappear, but the 
angle involved will depend on the geometry of 
the plates as well as the angle of contact. As the 
meniscus is stationary the components of forces 
due to suction and surface tension, acting parallel 
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to the direction of movement of the meniscus, 
will be equal in magnitude and opposite in direc- 
tion. Therefore when the equilibrium of the 
meniscus is being considered the two terms must 
be equated and when the pull between the two 
plates is required they must be added. The 
principle illustrated in this simple case is an 
important one and should be kept in’ mind 
when we come to discuss the work of Haines 
[2, 3, 4) and Fisner [5, 6]. 

Consider the two straight parallel 
illustrated in Fig. 1 (a). If the water is withdrawn 
through a small hole in the centre of the right- 
hand plate then the suction will remain constant 


plates 


(c) 


and the pull, P, will fall, due to the decline in 
the wetted area, A. Suction vs. moisture and 
pull vs. moisture graphs may be constructed 
for this system and these are given in Fig. 1 (b). 
The system may be elaborated a little and the 
plates be bent near their centres (see Fig. 1 (c)) 
so as to reduce the local distance between them ; 
the suction vs. moisture relation will then corres- 
pond closely to the drainage curve for beds of 
uniform particles. 

Figure 1 (d) shows the suction vs. moisture and 
pull vs. moisture graphs for the pair of curved 
plates. 

Such simple models cannot be used too freely 


Amount of water 
between the pilates 
‘ increasing 


| =2AT(* +87 sna 


between 
pilates 


Amount of water 
between plotes 


ng 


Armmount of water 
between pictes 
increasing 


Amount of woter _ 
Detween picates 
increasing 


Fic. 1. The pull between two wetted plates. 
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when discussing capillary effects in porous beds, 
but, as the mathematical procedures used by 
Harxes and Fisner are essentially the same as 
those just used, these two workers might have 
saved themselves much argument if they had 
returned to such fundamental Fisnen 
would certainly have quickly discovered that his 
basic ideas were untenable. 


Cases, 


Wonk 


Haines [2] calculated the pull between two 
spheres caused by a water droplet between them. 
This pull was then related to the corresponding 


3. Previous 


face area on the outside of a bed of regularly 
packed spheres and so the pull per unit area of 
bed exterior, or cohesion, was evaluated. The 
method used by Hates when calculating the pull 
was essentially the same as that discussed above 
and illustrated in Fig. 1, but the calculations 
were more complicated because the plates were 
not flat but parts of the surfaces of two spheres. 
He neglected the “ surface tension term in his 
calculations that the 
contact of the water was zero. 
Haines was able to estimate the cohesion only 
when the water was present as discrete droplets 
at the points of contact of the spheres (the 
pendular state) and when the water completely 


and assumed angle of 


filled the pores (saturation). He had also to 
assume that the bed consisted of uniform spheres 
in regular rhombohedral or cubical packing. 
His results are given in Fig. 2 together with an 
experimental curve of cohesion vs. moisture 
content for a bed of ignited silt particles. The 
device used to estimate cohesion was devised by 
Arrernerc [7] and measures the force required 
to split a block of test material by driving a wedge 
into it. Haines concluded from his investigations 
that the cohesion would rise continuously with 
an inercase of moisture and reach a maximum at 
saturation. In his first paper [2] Haines con- 
sidered that this maximum value would equal 
the entry suction. The entry suction is defined 
as the suction that must be applied to a bed of 
particles in order to cause air to begin to enter 
the pore spaces in the bed. 

Fisner [5] added precision to Haines’ approxi- 
mate caleulations by including the “ surface 
tension term omitted by him, but considered 
that the cohesion would fall with increasing 
moisture to reach zero at saturation. The follow- 
ing equation was produced [6] in an attempt to 
disprove Haines’ results, it is applicable only to 
uniform spheres. 


( (4, — a;) 


oheson ct soturation for 
spheres in closest packing 


£ 


Cohesion ot 
saturation for 


spheres in most 


open packing 


Bed saturated 


Moisture content of bed 


Experimental values obtained 
with Atterberg apparatus 


Fic. 2. Haines’ results on cohesion (2). 
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Where ¢ is the cohesion, or shrinkage force per 
unit area of bed exterior, 

éq is the suction in the water films, 

r is the radius of the spheres composing 
the bed, 

A, is the projection of the drained areas 
on the outside of the bed per unit area 
of bed exterior and 

a, is the projection of the drained areas on 
the inside of the outermost layer of 
spheres, per unit area of bed exterior. 


Fisuer proceeded to argue that the suction, dy. 
must increase continually with decreasing moisture 
for a series of stable states (this had been shown 
experimentally by Hatnes) and therefore to 
obtain a rise in cohesion with the addition of 
water (4, — a;) must increase quickly with the 
addition of water, which he considered to be 
unlikely. 

Hates [3] interpreted the theorem in a different 
manner ; he discussed the case of saturation in a 
bed of spheres in closest or rhombohedral packing 
with air just about to enter the top layer of 
pores. The bed was then subject to its entry 
suction of 12-9 T/r and a; was zero. Hatnes 
estimated A, to be 0-75 (he did not indicate how 
he obtained his figure) and on substituting for 
the various quantities in equation (2) the cohesion 
became 107'/r. The acceptance of this value by 
Haines showed that he had abandoned his 
original concept of the equality between cohesion 
and applied suction at saturation. 

Fisuer also refuted the experimental results of 
Hanes by asserting that the Atterberg apparatus 
measured the work required to break the test 
sample rather than the cohesion and he showed 
that this work increased with increasing moisture 
over the pendular range. 

The results of the controversy between Haines 
and Fisier was to cast doubt upon the ability 
of the Atterberg apparatus to measure cohesion 
and to persuade Hatnes to abandon his original 
idea of the equality between cohesion and 
applied suction at saturation, Although Haines 
upheld the concept of falling cohesion with falling 
moisture, a method of calculating the cohesion 
had still to be found for the wide range of moisture 
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contents between the pendular state and satura- 
tion. 


4. A RE-EXAMINATION OF FIsSHER’sS 
EQuATION WHEN APPLIED TO THE 
SATURATION CASE 


It will be useful to restate the equation 


*) (4, — (2) 
r 


(44 + 
where A, is the projection of the drained areas 
of the outside of the bed per unit area 
of bed exterior, and 
a, is the projection of the drained areas on 
the inside of the outermost layer of 
spheres per unit area of bed exterior. 


Consider the top view of a bed surface as 
depicted in Fig. 3(a). The arbitrary boundary 
encloses an area of unity and a drained area of 
A,. As the bed is saturated a; must be zero. 
Within the boundary the projected water surface 


Draned geo A 


Water area (1-A,) 


Boundary enclosing unit 
area 


\ 
< 
Unit cell 
Fig, 8. Air-water interfaces in beds of uniform spheres at 


saturation, 
8 (a) Uniform spheres in random packing. 
8(b) Uniform spheres in rhombohedral 
packing, projected area of spheres 
per unit cell area = 0-91 (A,, on 
air entry = 0-86). 
8 (c) Uniform spheres in cubical packing, 
projected area of spheres per unit cell 
area = 0-77 (A,, onair entry = 0-71). 
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area is (1 — A,) and the “surface tension” 
force acting normal to the bed surface, according 
to Fisner’s argument [6], is (2. 7. A,r). This 
force also upon the air water interface 
within the boundary and, since a state of equili- 
brium is assumed, it must be balanced by the 


acts 


suction in the water acting over the projection of 
the water surface area, (1 — A,). The condition 
for a stationary interface is therefore : 


by (1 A,). (3) 


At saturation Fisuenr’s equation becomes : 


c¢= ba A, 
and so we see that 
e by A, T ba (1 


Thus the cohesion is equal to the applied suction 


A.) = b4 


at saturation, 

It is interesting to examine the case of spheres 
in closest packing when air is about to enter the 
bed. The entry suction as given by Haines [3] 
is 12-0 Tr, Substituting this in equation (3) we 
may obtain A 


A, 


bac 


This is higher than the value of 0-75 quoted by 
Haines. As Haines gave no indication of his 
method of calculating A, the source of his error 
is unknown. Figs. 3(b) and 3(c) illustrate 
the cases of air entry for rhombohedral and 
cubical packing. 


RELATIONSHIP 
AND Sucrion DEVELOPED 
A Porous Bep 


5. «Tur 
COHESION 
BY 
Fisner’s equation can only be used to estimate 
cohesion when the configuration of the water 
films within the bed is accurately known, i.c. 
with uniform spherical particles in regular packing 
and with the water present as discrete droplets 
at the points of contact between the particles. 


When the bed is composed of particles of varying 
size and shape and in random arrangement 
Fisuer’s equation is inadequate and some other 
method for calculating cohesion must be found. 


Much insight into the nature of the water films 
within porous beds may be gained by studies of 
their suction vs. moisture relationships. For 
example the suction developed by a porous bed 
immediately gives the curvature of the water 
films within it, and the work done when extracting 
a given quantity of water from the bed provides 
a clue as to the change in surface area of the 
water films accompanying this withdrawal. Before 
considering these area changes and their implica- 
tions in more detail we must discuss the important 
work of Haines [4] in connection with the 
suctions developed by porous beds. 

Haines prepared saturated beds of particles 
and subjected them to steadily increasing suctions. 
Some typical relationships applied 
suction (called suction potential by Hares) and 
moisture content are shown in Fig. 4. The mois- 
ture content is denoted by a ratio, M, where 


bet ween 


Moisture cantent 
of bed, M 


(a) 


Fic. 4. Typical suction vs. moisture content curves for 
porous beds. 
4(a) Curve for a bed having uniform pores. 
4(b) Curve for a bed having non-uniform pores, 
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ul total volume of water in the bed pores 
total volume of the bed pores 


thus at saturation M is equal to unity. The curve 
in Fig. 4 (a) is for uniform random packed spheres 
and it illustrates the phenomenon of hysteresis 
which is discussed shortly. Fig. 4 (b) refers 
to sand of non-uniform particle size and shape. 
Haines was the first to point out that the shape 
of the suction-moisture curve indicated the pore 
size distribution of the bed. A flat curve, Fig. 4 (a), 
shows that there are a large number of pores of 
similar size and entry suction, whereas a sloping 
curve, Fig. 4(b), suggests that there is a wide 
variation in the size of the pores. 

Haines also calculated the suctions developed 
in some of the water films occurring in beds of 
regularly packed spheres. ‘Two important values 
are the entry suctions, ¢,,. for rhombohedral 
and cubical packing. He estimated these by 
calculating the size of an air bubble about to 


Ry = 


(b) 


Fie. 5. Air entry for beds of spheres in regular packing. 
5(a) Air bubble about to pass between spheres in 
rhombohedral packing. 

5 (b) Air bubble about to pass between spheres in cubical 
packing. 


force its way between the spheres as illustrated in 
Fig. 5. The radii of the air bubbles were 0-155 r 
and 0-43 r for rhombohedral and cubical packing 
respectively. The entry suctions were thus 
calculated as 12-9 T r and 46 Tr. For random 
packed beds of uniform glass spheres, called 
‘glistening dew, of diameter 360, Hatnes 
found by experiment that the entry suction was 
around 67 

Different suctions are possible at the same 
moisture content depending on whether the 
water is being drained from the bed or imbibed 
into it. This phenomenon is known as hysteresis 
and it is illustrated in Fig. 4(a). The suction 
under drainage is always greater than the suction 
during imbibition for a given moisture content. 
There are two types of hysteresis loop: the 
largest or major loop, XNYZX, and the minor 
loops, of which there can be any number, xyzx. 


The free energy of the surfaces within a porous bed 


The increase in total surface energy of a system 
accompanying the isothermal creation of unit 
surface area of that system is 

dT 


Z=T — R — 5) 

dR 
where T is the free energy of unit surface area 
and may be identified with the mechanical work 
performed, 


RdT dR is the heat absorbed from the 
surrounds to maintain the temperature 
of the system and 

R is the absolute temperature. 


All drainage and imbibition processes to be 
discussed are assumed to take place isothermally, 
therefore the mechanical work performed on or 
by the bed will be associated with the term T in 
equation (5) and R. dT dR will be neglected. 

There are three types of surface within a 
porous bed : solid—gas, solid—liquid and liquid—gas, 
and the free energies of unit area of these surfaces 
are denoted by the symbols T,,, T,, and T,, 
respectively. At saturation the free energy of a 
porous bed is 8S,’ 7T,,*, where S,’ is the solid 

*This neglects the relatively small solid-gas and liquid— 
gas areas present at saturation. 
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surface per unit volume of bed. After air entry 
has taken place suppose that a fraction n’ of 
solid surface is drained and that the water 
surface is H then the increase in free cnergy is 


F = n' So —n' Sy Tu r HT,, 
If the angle of contact of the liquid is zero then 


and 


Ty, (n So (6) 


At saturation F will be zero and at complete 
dryness F will have a maximum value of T,, S,’. 

In future the term, F, describing the increase 
in the free energy of the surfaces within a porous 
bed accompanying the removal of water will be 
referred to as the energy of the porous bed. 
Changes of energy accompanying 
content may be deduced from 


changes of 
moisture the 
capillary characteristics of the bed. 

The areas under the drainage and imbibition 
lines of a suction vs. moisture content graph are 
a measure of the work performed during these 
processes. As pointed out by Hairs [4] the 
presence of an hysteresis loop indicates that work 
has been wasted in the drainage cycle and so the 
energy of the bed must always be less than the 
drainage work down to a given moisture content, 
The converse will be true during imbibition. A 


unique relation between the energy of the bed 
and the moisture content is wholly consistent 


with hysteresis but the water films have different 
suctions and thus different shapes during drainage 
and imbition and so an energy loop is more 
likely. 


Constructing the energy loop 

Fig. 6 (a) shows a graph of suction against 
moisture content for a bed of random packed 
uniformly shaped spheres. Fig. 6 (b) shows a 
graph where units of work or energy are plotted 
as ordinates and moisture content as abscissae. 
The energy of the porous bed and hence the work 
done on it or capable of being done by it will be 
zero at saturation (M equal to unity). 

When the porous bed is drained from saturation 
to a moisture content M, the work done on the 
bed is equal to the area under the drainage curve 
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(a) 


Moisture content, M 


Fic. 6, Calculating the energy changes in wetted porcus 
beds. 

Suction vs. moisture relation. 

Energy vs. moisture relation. 


6 (a) 
6(b) 


of the major hysteresis loop down to the moisture 
content .M, and is represented by the point Wa, 
in Fig. 6(b). The existence of hysteresis loops 
indicates that energy is wasted in the drainage 
cycles and is suflicient reason to assume that the 
energy of the bed under the suction $4, is less 
than the work done up to that suction, Wa, 

Now we can also reach the suction $4, from the 
suction ¢, via the drainage curve of ‘the minor 
hysteresis loop between them (see Fig. 6a). Let 
us assume that the energy of the bed when the bed 
is imbibing water at a suction ¢;, is Fy which 
will be below the drainage work. vs. maktese 
content line 0-W,. (see Fig. 6b). If the bed is 
drained from the moisture content M, down to 
the moisture content M, the suction will increase 
to ¢g, and the work performed on the bed will be 
given by the area under the drainage curve of 
the minor hysteresis loop and will be represented 
in Fig. 6b by the jump from Fy to Wa, Both 
Wa, and W, ° correspond to the same configura- 
tions of the ‘water films but W,.’ is the smaller 
term and comes closer to the actual energy of the 
bed under the suction $4,. If the bed is allowed 


§ 
Moisture content, M 
be Wo. 
10 
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to return to the moisture content M, and the 
suction ¢; it will perform the amount of work 
given by the area under the imbibition line of the 
minor hysteresis loop se ag in Fig. 6 (b) as 
the drop from Wie ‘to F;. Because the minor 
hysteresis loop is narrow the energy wasted 
during the cycle from 4; to 4, and back again 
is small compared to the total energy changes, 
and so it can be assumed, with but little error, 
that the increase in energy of the bed accompany- 
ing a change in moisture conditions represented 
by the extremities of a minor hysteresis loop is 
equal to the arithmetic mean of the areas under 
the drainage and imbibition curves of that same 
minor hysteresis Energy increments 
caleulated in this way can be put together to 
construct a complete energy loop and Fig. 7 
illustrates the necessary steps. 


loop. 


M, Me My Mg Ms M, 
Moisture content 


M, M, ™M. 
Moisture content 


Fic. 7. Constructing the energy loop. 
Starting at the origin the first part of the 
drainage curve is found to be reversible (as 
indicated by the double arrow) and so the incre- 
ment in the energy of the bed is given by the 
drainage work line (Wy) up to F,. The minor 
hysteresis loop between 1 and M, may now be 
used to fix the point Fa, as explained in Fig. 6 (b). 


Having positioned Fy, we may next draw a 
straight line between F, and Fz, and assume that 
it represents the energy vs. moisture relationship 
for drainage. If it is further assumed that the 
line F, - Fy, may be extrapolated then the 
position of F a, may be determined. The accuracy 
of such an ext rapolation depends on the distances 
involved and will be greater the closer are the 
minor hysteresis loops. By using the minor 
hysteresis loop between M, and M,, F; may be 
positioned from Fy and F; from FP, by extra- 
polating the line 0- P. . Fy may be obtained from 
F,, from the third minor loop and if we continue 
in this stepwise manner the complete energy loop 
can be constructed when sullicient hysteresis data 
are available. As explained previously when the 
hed becomes completely dry the energy will have 
reached a maximum value of 


A general cohesion vs. moisture relationship 

Until now we have used only the moisture 
content, M, as the abscissae in all the graphs. 
When discussing energy changes in porous beds 
we should, to be rigorous, use the volume of 
water withdrawn from the bed as abscissae and 
not moisture content, but the two quantities are 
directly proportional to one another in beds of 
constant porosity and so the previous arguments 
are quite valid. However during the following 
derivation of the cohesion equation we shall have 
to consider the effects of small changes in bed 
porosity and so, for the moment, we must replace 
the moisture content, M, by the volume of water 
withdrawn, Vi. The two quantities are related 
by the expression 


V, = ¢ — M) (7) 


where ¢ is the porosity of the bed i.e. the ratio: 
total volume of pores/total volume of the bed 
and J; is the total volume of the bed. 

In this paper we are primarily concerned with 
shrinkage forces and so the cohesion equation 
will be derived assuming that the bed is under a 
drainage suction. For mathematical purposes 
the drainage energy characteristics of a porous 
bed are conveniently represented as a plot of 
(dF, dV.,), as ordinates against V, as abscissae. 
(dF, dV,,) has the dimension of suction and as it 
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is closely related to the drainage suction it will be 
given the symbol ¢,’. 


Vw 


f withcrown, 


water 
(b) 
Fic. 8. Stretching a wetted porous bed. 

8(a) The two porous beds. 
Energy characteristics of the two porous beds. 


Voumne 


(b) 


Two porous beds in slightly different capillary 
and energy states are depicted in Fig. 8 (a). 
They are identical in every respect save porosity. 
They each contain the same number of particles 
the same volume of water and the pore con- 
figurations are such that bed 1 may be made 
identical to bed 2 by an elemental increase in 
volume 6 vy. The shrinkage forces acting on the 
beds may be deduced in the following way. Both 
beds are saturated with water and then subject 
to a series of drainage and imbitition cycles which 
enable graphs of ¢,° against Vi to be prepared 
for each bed (see Fig. 8(b). The beds are then 
saturated with water again and drained to the 
conditions illustrated in Fig. 8 (a) and represented 
by the abscissae + in Fig. & (b). 
From the areas under the curves of 44° against 


V. we can deduce the difference in energy 
between the two beds, 6 Fy, and this will equal 
the mechanical work which must be done against 
the shrinkage forces when bed 1 is given a rever- 
sible isothermal stretch so that it becomes 
identical to bed 2. If the beds are assumed to 
have unit cross-sectional areas perpendicular to 
the direction of stretching then the mechanical 


work done will be ¢ dV, 
From Fig. 8 (b) we see that 
5 Fy areaO CB(V, + — aread AV, (8) 
As the area 0 CV, is common to both curves 

6 F, areal CB(V, + 


Both 
therefore 


area 0 AC 0 (9) 


beds have the same volume of water, 


and so 


area V,.C.B.(V, + 5V,) (10) 


The area 0 AC 0 is equal to the mean value of 
5¢4 over the range 0 to Vi, multiplied by V,. 
The problem is now to calculate 46 4,’ at a constant 
value of the water withdrawn. Consider the graph 
in Fig. 9, the two vertical lines represent the same 
moisture content Win both beds and the left 
hand line is the position at which we must 
estimate 5¢,'. From Fig. 9 we see that 6 ¢,' 
consists of two parts: that caused by an increase 


Bed | 


($a) 


oats 
#s," ™) 


Volume of water withdrawn 
Fic. 9. Estimating 5¢,’. 


‘ 
ac 
Air Air 
Water | Water 
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of porosity d5¢ at constant moisture content and 
that caused by a reduction in the water withdrawn 
at a constant porosity e + de. Let us _ first 
consider the effect of changing the porosity. 
It will be assumed that when two beds have the 
same moisture content the shape of the water 
films in the pore spaces will be similar and that 
the increase in the radius of curvature of the 
water films will be similar to the increase in the 
mean pore radius accompanying the increase in 
porosity. It has been stated [8] that for a given 
particle size and shape 


constant 


therefore (“ constant $a 

de ly e(1 — e) 
Thus an elemental increase in porosity of d¢ will 
reduce the suction values by a_ fraction 
dee (1 — e). It is reasonable to assume that the 
values of 4, will be similarly reduced and so 


( de —¢) (Ht) 


and the first part of 5¢,' is dy ¢(1 —e). The 
second part is given by the product of the gradient 
of the lower curve at the abscissa V, and the 
separation between the two vertical lines. The 
gradients of the upper and lower curves are nearly 
the same and so 


(39)... (302). 


The separation between the vertical lines is 
(e + de) (1 — M) (V, + —e(1— M) Vy 


now 5 V,/Vy, equals de/(1 —e) and so the 
separation is 


+800 - M) Vi, (1 


) ~e(1 — M) Vy 


which can be simplified to (de, 1 — (1 — M) Vz 
by neglecting terms of the second degree. There- 
fore the second part of 5¢,' is (dd,'/d V,),. 
de/(1 — e) (1 — M) Vy and so 


be 
5d, = by e(1 e) 
_be a 


— M)V, (12 
aVJ. l—e Vy (12) 


now equals (1 — M) Vyand so V,), 
can be replaced by — (d dy'/dM),.1 ¢ Vy. 
Therefore 

e(1 — e) 


d bq be 
— M) (18) 


bbq) = by 


Taking mean values over the moisture content 
range 1 to M and denoting these by the sub- 
script m we have 


be 
e(1 —e) 


d bq be 
iM) * — ¢) (14) 


and the area 0 AC 0 is equal to (5¢,'),,.. Vj, and 
since equals e(1— M)V, and Vy, Vy 
equals 6¢ 1 — e we have 


From equations (9), (10) and (15) we find that the 
difference in energy between the two beds is 


bF, ba Vz (1 M) 8 Vp + 


(1 


(15) 


dM !, 


which is identical to the mechanical work per- 
formed stretching the bed by a volume 4 Vy Le. 
c 6 Vy, where ¢ is the cohesion in the porous bed. 
Thus the cohesion — moisture content relationship 


is 
= — (bam (1 — M) + 
+ a — (1—M) (17) 


where ¢, equals (dF, or 
(dF,/dM), and the subscript m denotes a 
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Fig. 11. _ 12. The cohesion in a bed having uniform pores. 
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mean value over the moisture content range 
from 1 to M. 


Discussion of the cohesion moisture content equation 


At saturation M equals one and so ¢ becomes 
equal to ¢,'. Over the earliest part of the suction 
moisture content curve the drainage process is 
reversible and therefore ¢,' is identical to ¢, the 
drainage suction. Hence at saturation the 
cohesion is equal to the applied suction. 

In order to demonstrate the effects of the 
various terms in equation (17) three sample 
computations are illustrated in Fig. 10. As the 
og versus M curves bear a resemblance to their 
parent é, versus M curves the drawings illustrate 
the effects of pore size distribution on the cohesion. 
With beds having uniform pores the last term in 
in equation (17) is usually of little importance 
because the é, versus M graph is flat (see Fig. 10 
(a)). However, if pores of widely varying size 
are present, as is the case in Figs. 10(b) and 
10 (c) then the last term in equation (17) must be 
included ; the error incurred by omitting it is large 
at the lower moisture contents. 


Graphs of F against M and ¢,' against M were 
estimated for a sample of glass spheres of diameter 


360 from suction moisture content data 
provided by Haines. The original graph of 
dg and ¢,; versus M has been corrected for the 
effect of bed thickness. Unfortunately only three 
minor hysteresis loops are given and so the 
determination of F and ¢,' is bound to be approxi- 
mate. The results are shown in Figs. 11 and 12. 
When calculating the cohesion-moisture content 
relationship from the data in Fig. 12 (a) the last 
term in equation (17) was omitted. The omission 
was thought to be justified because over the 
majority of the moisture range the gradient of 
the ¢4° versus M graph was zero and just after 
air entry when the gradient had a large but 
uncertain value the factor (1 — M) was small 
and appeared raised to the power two in equation 
(17). A plot of cohesion versus moisture content 
calculated from Haines’ data is given in Fig. 12 (b) 
together with Haines’ own measurements of 
cohesion taken from [2]. Also given in Fig. 12 (b) 
are measurements of cohesion taken with 
the special tensile testing apparatus depicted in 


Fulcrum of tensile testing 
apporotus 


Hard fibre gasket 
retracted downwards 


Wegnt to amply 
tensile bad to bed 


Gasket retracting screw 


Gasket 0-0625 in thick 


Fic. 13(a) Tensile testing device. 
Fre. 13 (b). Details of mould for porous bed. 


Fig. 13. The tests were made upon moist beds of 
uniform glass spheres of 170 «* diameter formed 
in a horizontal dumbell shaped mould which was 
constructed in two halves. Whilst water was 
being sucked from the bed of spheres the two 
halves of the mould were clamped together with 
a hard fibre gasket between them. When the bed 
had been partially dried the suction was stopped, 
the clamping pressure released and the gasket 
retracted. The bed of spheres was then stretched 
and broken by placing weights on the scale pan 
(see Fig. 13). Except with saturated samples, 
where some necking was observed before fracture, 
the beds broke cleanly across the centre. There 


, *Glass spheres of this size were conveniently available 
at the time of the experiments. 
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is wide disagreement between the two sets of ex- 
perimental measurements presented in Fig. 12 (b) 
but no fruitful comments can be made about 
this at present because different types of porous 
bed and different pieces of apparatus were used 
in each set of experiments. 

The results in Fig. 12 (b) shows that the theore- 
tical cohesion- moisture content equation (equation 
17) is basically of the right form; but clearly 
much more experimental work is necessary to 
obtain more accurate measurements of cohesion, 
and to obtain more information about the 
hysteresis loops in the suction potential moisture 
relationships of beds with uniform and non- 
uniform pore size distribution. 


CONCLUSIONS 
As a result of this study we may conclude that : 


(i) the cohesion in a wetted porous bed is 
equal to the applied suction, when the 
bed is saturated, 

(ii) the cohesion will fall below the applied 
suction after air entry into the bed, and 


(iii) the precise relationship between cohesion 
and moisture content will depend on the 
pore size distribution in the bed. 


NOTATION 


a, = the projection of the drained areas on the inside 
of the outermost layer of spheres, per unit area 
ot bed exterior 

A = the wetted area of one of a pair of parallel plates 


A, = the projection of the drained areas on the outside 


B 
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Exact forms of the unrestricted Gibbs=-Duhem equation 


H. C. Van Ness* 
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Abstract—The exact forms of the Gibbs Duhem equation in terms of fugacities and activity 
coeflicients, applicable to any multicomponent phase, are developed in a straightforward manner. 
Reduction to restricted equations valid for constant temperature, constant pressure, or constant 
composition becomes a simple matter of dropping terms. 


Résumé Loauteur développe de fagon directe les formes de Véquation de Gibbs-Duhem en 
fonction des coefficients de fugacite et @activite, applicable & mimporte quelle phase plusieurs 
composants. Le passage aux équations restreintes, valables pour une température constante, 
une pression constante, ou une composition constante, se résume a une simple suppression de 
termes, 


Zusammenfassung —Die exakten Formen der Gibbs-Duhem-Gleichung in Ausdriicken der 
Fugazititen und Aktivitatskoeflizienten, die auf cine Phase mit mehreren Komponenten anwendbar 
ist, werden in einfacher Weise entwickelt. Die eingeschriinkte Gleichung fiir konstante Temperatur, 
konstanten Druck, oder konstante Zusammensetzung erhalt man durch einfaches Streichen von 


Ausdriicken. 


The much-quoted paper by In. and Doncr [3] 
which appeared in this journal a few years ago 
presented derivations of the exact forms of the 
Gibbs-Duhem equation in terms of fugacities 
and activity coefficients for binary systems in 
phase equilibrium at (a) constant temperature, 
and (b) constant pressure. The object of this 
paper is to generalize these equations to apply 
to any phase, whether in equilibrium with another 
phase or not and without restriction as to the 
number of components or constancy of tempera- 
ture or pressure. Moreover, the derivations are 
direct and without complication. Thus they 
allow a uniform, yet simple, treatment of all 
cases, 

We start with the basic equation for an open 
system (see, for example, Dennicn [2}): 


Vdp —SdT (n; (1) 


This equation applies to the whole phase ; 
indeed it applies to the whole of any multiphase 
system, provided that equilibrium exists through- 


out the system. However, for the purposes of 
this paper it will be applied to single phases, and 
will be put on a mole basis by dividing through by 
the total moles in the phase : 


vdp —SdT = (a; (2) 


where v is the molal volume and S is the molal 
entropy of the solution. The x,’s represent mole 
fractions, and jy; is the chemical potential of 
component 7 in solution. It is of course identical 
with the partial molal Gibbs free energy. 

The defining equation for the fugacity of a 
component in solution is: 


du, = RT dinf (constant temperature) (3) 


This equation takes into account variations of ,; 
as a result of pressure and composition changes 
at constant temperature, and hence cannot be 
substituted in equation (2) in the general case 
where du; is the total derivative of yu; as a result 
of changes in temperature as well as in pressure 
and composition. Rather, we integrate equation 


*During 1958-59, Visiting Fulbright Lecturer in Chemical Engineering at King’s College, University of Durham, 
Newcastle upon Tyne. 
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(3) at constant temperature from the low-pressure, 
ideal-gas reference state of the pure component 
to the actual state of the component in solution : 


= RT \n f, RT \n 
where the asterisk designates the reference state 
for which, by definition, /;* = p*. Hence 

RT In f, RT In p* (4) 
Since p* is the constant reference-state pressure, 
general differentiation of equation (4) gives 
dy; : dF * RTd\n f, + Rinf,dT Rinp* dT (5) 
In this equation dy; is the required total deriva- 
tive. Equation (5) may be simplified by the 
following substitutions. Since F;* is a function 


of temperature only 


S* dT 


and by equation (4) 


* 


Therefore 


dy; —S*dT 4 


dT RT dif, 
or since S* H* T from the detining 
equation for free energy, 


du, —(H#/T)dT + (nu, T)dT + RT dlnf, (6) 
If equation (6) is multiplied through by 7, and 


summed over all components /, 


* 

= — ap 

T 
But 2(7,H*) = H*, where H* is the molal 
ideal-gas enthalpy of the solution, and 
2 (x; 4) = F. where F is the actual molal free 

energy of the solution. Thus 


(a, dp;) = —(H*/T)dT + (F/T) aT 
t RT S(a,dinf,) (7) 


dT RT 2 (a, d In fi) 


Combining equation (7) and (2), 


vdp —SdT — —(H*/T) aT + 
(F/T) aT + RT (a,d nf) 


H. ©. Vaw Ness 


or since F 
v (H* — 
RT dp + RT? dT = (a,d In (8) 


Equation (8) is the general equation relating 
pressure, temperature, and the fugacities for any 
liquid or vapour phase, whether in equilibrium 
with another phase or not. Reduction to the 
special cases of constant pressure or constant 
temperature or both is just a matter of dropping 
the appropriate terms. 

Equation (8) may also be used as the starting 
point for the further development of both specific 
and general equations. Particularly useful is an 
equation relating pressure, temperature, and the 
activity coeflicients for any phase. Equation 
(8) is of course valid for a pure material as well as 
for a solution. If it is written for pure i at the 
same temperature and pressure and in the same 
phase (liquid or vapour) as the solution in which 
i is a component, we have : 


(H* —H) 
' / d 
rT! re 


Multiplication by 2; and summation over all 
components of the solution gives : 


gy. — Hd op 
RT RT? 


din f; 


(9) 


Subtracting equation (9%) from equation (8) 


v — 2 dp - 


But v — Ya; 2,) is the volume change of mixing 
Av, H — (a, H,) is the heat of mixing AH, and 
[ad In (f, — (ad In y,). The last equality 
results from the fact that In 1/2) = 0. 
Thus the final result is 


dp - as dT = ln 


RT? 
Equation (10) applies to any phase provided 
that Av, AH, and the y,’s are all taken with 
respect to the pure components at the same 


(10) 


| 
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temperature and pressure and in the same 
physical state as the solution (often regarded 
as standard states). Applied to phases in vapour- 
liquid equilibrium, this usually requires the 
consideration of fictitious states for at least one 
of the pure components. In this case equation 
(10) is ordinarily written for the liquid phase. 
Since the properties of liquids are usually insensi- 
tive to changes is pressure, values of Av and AH 
measured at pressures high enough so that all 
components are stable as liquids can be used in 
equation (10) at the equilibrium pressures in 
question without introducing sensible error. 
The equations resulting when equation (10) is 
restricted to conditions of constant pressure or 
temperature or both are obvious. 
One additional step is of interest. Since 


d ay ln dy, + In 


then =d L(a, In y,) — (in 


But AF®/RT y,) where AF® is the 
excess free energy of mixing. Therefore, equation 
(10) may be written : 


A ) Av AH 
dp — aT + (11) 


There are several uses of this equation. Note 
first that d(AF® RT) is an exact differential. 
Hence the various reciprocity equations may be 


applied to give relationships between — the 


differential coefficients. For example, it is 
easily shown that : 
iT 


This equation permits the determination of the 
effect of pressure on the heat of mixing from 
measurements of the volume change of mixing. 

Another use of equation (11) is to derive the 
various “area” tests which may be used to 
examine data for thermodynamic consistency. 
Since d(AF*®/RT) is an exact differential, 


integration gives : 


Exact forms of the unrestricted Gibbs-Duhem equation 


(ier).~ (rer), 


2 

Av AH 


The limits 1 and 2 are here taken to represent 
any two states of an open system, and integration 
is of course carried out for some specific path 
for which data are available. This way of writing 
the equation illustrates the infinite possibilities, 
but in practice integration is ordinarily carried 
between limits and along paths which result in a 
considerable simplification of the equation. For 
example, if the composition of a binary solution 
is varied from x4 =0 to v4 =1 at constant 
temperature and pressure, equation (13) reduces 
to the Repitca and Kisrer [4] equation : 


If data are taken for a phase at constant 
composition, the last term of equation (11) is 
zero, and the equation may be solved for AH: 


E 
AH TAv ( (AF 
aT 


This is the equation used by Asuiey and Brown 
[1] to caleulate heats of mixing from phase- 
equilibrium data for the HCl C,H, system. 

The main object of this paper has been the 
development of the general relationships given 
by equations (8) and (10). In addition, a few 
illustrations of their use have been pointed out. 
However, no attempt has been made to exhaust 
the possibilities of their application to specific 
cases. 

NOTATION 

Jf, = fugacity of pure component i at the T and p of 
the solution and in the same physical state 

J, = fugacity of component i in solution 

F = molal free energy (Gibbs) of the sotution 

F, = mola) free energy of pure component i at the T 
and p of the solution and in the same physical 
state 


H. Van Ness 


molal free energy of pure component ¢ in the 
reference state of an ideal gas at p* and the T 
of the solution 

excess free energy change of mixing 

molal enthalpy of the solution 

= molal enthalpy of pure component ¢ at the T 
and p of the solution and in the same physical 
state 

molal enthalpy of pure component i in the 
reference state of an ideal gas at p* and the T 
of the solution 

enthalpy change of mixing at constant T and 
p: Le., integral heat of mixing per mole of 
mixture formed 

moles of component i 

pressure 

reference-state pressure where all gases are ideal 
universal gas constant 


REFERENCES 


S-= total entropy of a phase or system 
S = molal entropy of the solution 


* 


molal entropy of pure component in’ the 
reference state of an ideal gas at p* and che 
T of the solution 

absolute temperature 

total volume of a phase or system 

molal volume of the solution 

molal volume of pure i at the T and p of the 
solution and in the same physical state 
volume change of mixing at constant T and p 
per mole of solution formed 

mole fraction of component ¢ in solution 


activity coeflicient of component 7 in solution; 
equal to aj fj 

chemical potential or partial molal Gibbs free 
energy of component ¢ in solution 
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Abstract 


Mathematical expressions are derived for the residence time and the distribution of 


residence times of falling particles in a turbulent flow system. 


Résume 


L’auteur déduit des expressions mathematiques pour le temps de résidence et la 


répartition du temps de résidence de particules tombant dans un systéme a écoulement turbulent. 


7 


@ -Fiir die Verweilzeit und die Verteilung der Verweilzeit fallender Teilchen 


in einem turbulenten StrOmungssystem werden mathematische Ausdriicke abgeleitet. 


INTRODUCTION 


Tue mixing and staging of flowing tlhuids are of 
interest in many physical and chemical processes. 
Various articles in this journal [2, 3, 4, 5, 6, 9) 
testify of significant progress in the analysis and 
understanding of these phenomena, In all these 
studies, however, only the mixing of a continuous 
phase in the flowing systems was considered. 
In many two-phase processes not only the mixing 
of the continuous phase is of interest but also 
the behaviour of the dispersed phase. 

In the mathematical approach of the longi- 
tudinal mixing of continuous tlow systems, the 
diffusion equation has been used. This equation 
can also be applied in dispersed flow systems 
[7, 8]. From this equation expressions can be 
derived for the residence time and for the distribu- 
tion in residence times of dispersed particles. 
In particular the case will be considered where 
particles are falling in a fluid which is moving 
countercurrently, The particles are introduced 
into the system at the top and withdrawn at the 
bottom, below the continuous phase inlet. The 
rate of withdrawal can be influenced by physical 
means, for instance by controlling the settling 
rate and the rate of coalescence. 


Tue EQUATION 


A vertical, one-dimensional flow system will 


DIFFERENTIAL 


*See list of symbols used. 


be considered. The diffusion equation then can 
be written :* 


d d 
pit _ 
d 


oY (1) 


It has been shown by Soo [7] that in turbulent 
systems the diffusion coetlicient D of the dispersed 
particles is about equal to that of the continuous 
Thus of the diffusion 
coellicient obtained by tracer injection in’ the 


phase. measurements 
continuous phase, also can be used for determina- 
tion of the dispersed phase diffusion coefficient. 

For simplification of equation (1) two assump- 
tions will be made : 


1. D is independent of co-ordinate y, which 
will be true for homogeneous turbulence. 
The particle velocity v is independent of 
co-ordinate y and also independent of 
particle concentration n. This is a rather 
sweeping assumption since it means that 
hindered settling effects are neglected, which 
is only permissible for a very low hold-up 
of the dispersed phase and low continuous 
phase velocity. 


The differential equation for a system with the 
above mentioned assumptions will then be : 


dy? (2) 


OL. 
LO 
on d 
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When we consider a state of equilibrium, 
where there is no variation of the concentration 
with time, equation (2) reduces to 

on 
dy or (3) 
dn 


D— — rm = 


dy constant —¢ 


The constant ¢ is the current of transportation of 
the particles (first boundary condition). 


Averacre Restvpence Time 


Integration of (3) yields: +Cew? (4) 
The constant C can be found from the second 
boundary equation, which will depend on the 
manner in which the particles are removed, 
Since this can be done in different ways, it may be 
possible to influence the hold-up of the dispersed 
phase by the manner in which the removal of the 
particles is realized. Derivation of a mathematical 
expression for the second boundary condition, in 
the case of a countercurrent operation, will be 
given with reference to Fig. 1. 


CONCENTRATION 
OF PARTICLES 


COEFFICIENT 
VELOCITY 


DIFFUSION 


SECTION 


! . 
2 
— 


Fie. 1. Schematic situation of boundary conditions, 


Dispersed phase is introduced at the top of the 
column (y= 0). The continuous phase is in- 
troduced at y= L. The particle velocity v for 
0 <y <L, will be the particle settling velocity 
k minus the fluid velocity. The value of the 
diffusivity will be D over the length of the sections 
land 2(0<y<L’'). Aty >L’ (section 3) the 
diffusivity coetlicient will be zero. 


At the interface y = L’, it follows from the 
condition : 


D —kn —kn, 
(y@L'—) 


that = 0. 
dy 


Hence for L <y <L', dn/dy = 0 and the con- 
centration equals q/k. At the interface y = L, the 
boundary condition becomes : 


—kn = —kn 


dy 


or D hen (5) 
dy (weL-) 

We have assumed above that at L and L’, the 
concentrations on each side of the boundary are 
equal, This assumption has been made, since we 
feel that continuity in the concentration will be 
ensured by the existence of diffusivity. 

An analogous boundary equation can be derived 
for the case of coalescence of the dispersed phase 
at y —L. The value & will then represent the 
sedimentation velocity of the coalesced particles, 
Solving equation (4) with boundary condition 
(5) we obtain : 


For y 


It can be seen that from (7) that np < q/v, if 
k > v. For large values of k, ny, +0. This con- 
dition may be considered as a limiting case. 
The particle concentration distribution along the 
axis then becomes : 


n=! — (8) 


v 


The average concentration of n over the whole 
length L of the column will be : 


v 


VOL. 
10 


dn 
dn 
D 
| 
v—k 
a=! “| (6) 
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where Pe = Lv D. The average residence time 
of the particlse 0 is given by 0 = (L/q) n,, 


or (10) 


In Fig. 2 the relation between the residence time 
6 and Pe, according to equation (10) is given. 


FLOW (key) 
~ 
10)(k=00) 


DISPERSED PHASE FLOW 
_ 


rT. 


| 


3 4 


| 
| 
| 
Pe/2=1v/20 


L=LENGTH OF APPARATUS 
v=PARTICLE VELOCITY 


Dimensionless residence time as a function of 
Péclet number. 


| 


Fig. 2. 


or RESIDENCE 
VARIANCE 


DistTkIBUTION 
Times, THE 

The variance o? will be used to characterize 
residence time distribution. The variance is 
calculated by VAN Laan’s method [5]. This 
is done with the aid of the transient differential 
equations. For section 1 (0 < y —L) we get: 


(2) 


For section 2 (L <y<L’) the differential 
equation becomes : 

ou 
dy? 


on 
=D (2a) 


k/v . 


For section 3 (y > L’) the diffusion coeflicient D 
equals zero ; hence 


on dn 

—==—k— 2b 

vt dy 

At the entrance of the column (y = 0) a pulse 

of tracer mathematically described by the Dirac 

impulse function is introduced, The total quantity 

of injected material per unit cross-sectional area 

equals q*. 

The initial and boundary equations now be- 
come : 

o- n=O 


(11) 


vu — 6 (t) 


on 
va — D 


on 
= —~DpD — g* 


L+ 


on 
kn 


My — = Mp4 


We are interested in the response R(t) at y = L 
to the 6(f) input. For the sake of simplicity we 
seek the solution for the case where L’ —L 
approaches zero, If R(p), the Laplace transform 
of R(t), given by 


a 


R(p) = fem 


is introduced, the solution of equations (2) 
taking into account the initial and boundary 
conditions (11) can easily be constructed. The 
solution R (p) reads : 


(12) 
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(13) 


The average residence time 90 and the variance 
o can be calculated from (12) by means of simple 
differentiation [1, 5] 


= Lim — (14) 


poo P 


+ = Lim R(p) 


p>o dp* 
In this way we find : 


6 Pe 


= Pe —1 (15) 


o® Pet 


i) +2 ) Pe 
(16) 


It is now possible to make a comparison with a 
system consisting of N ideal mixers. The com- 
parison is based on the assumption of equal 
variance of the two distributions of residence 
times. For a system of N ideal mixers in series 
we may write : 


N (17) 


Dividing the square of (15) by (16) gives us the 
number of equivalent number of ideal mixers. 
For finite k, limiting cases are : 


In the case of continuous phase flow, where 
k = v, (15) and (16) lead to the same solution as 
given by van per Laan [5]. 

It is of interest to note that the number of 
stages N approaches the value one, if Pe +0. 
The case of k infinite, which would lead to a value 
of 15 if Pe is trivial. Fig. gives 


Ve Pe=V2 iv/0 
=LENGTH OF APPARATUS 
v=PARTICLE VELOCITY 


Fic. 3. Relation between (@ ¢)* and diffusivity 
(dimensionless). 


(9 o)? asa function of the Péclet number. Devia- 
tions from the curve for the continuous phase 


flow are relatively small except for small values 
of Pe. 
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Pe 


a 


Pe=0 


a 


WS 


Biv k) + [2 — (av Pe 


_ 
Pe 


VOL, 
10 


where 
s= 
JP; 
Le 
Pe = — 
| 
| | / 
> 
> 
| 
vO vw. Pe v 
y “its 
vO 
Lk 
Pe 1 Iv ve 
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current of transportation of particles 
total amount of injected tracer material (L~?) 
time (T) 
particle velocity (LT) transform: R (p) fo R (1) dt. 
length co-ordinate (L) 0 
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NoTaTION 


longitudinal diffusion coeflicient (12 T-!) (t) = impulse function 
particle settling velocity (LT) R(t) = response function 
= length of tube (L) o® = variance of the spread in residence times 
= number of mixing stages — of the particles 
= particle concentration (L-) 0 = average residence time of particles 


= parameter in Laplace transformation 
Pe Péclet number = 


bar over function designates the Laplace 
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Abstract—The diffusivities of acetone, carbon tetrachloride and water vapours in air have been 

determined experimentally using the method first developed by WinkeLMANN and values have 

been obtained over a range of operating conditions to an accuracy of within + 2 percent. The 

evaporation of two-component liquid mixtures, consisting of a volatile and a non-volatile material 

has also been been studied in the same apparatus, It has been shown that, when the non-volatile 

liquid is the denser, the liquid remains completely mixed and the evaporation rate can therefore 

be simply calculated for an ideal mixture. When the non-volatile component is the less dense, VOL. 
convective mixing is negligible and the mass transfer in the liquid is a process of unsteady state 10 
molecular diffusion. 1959 


Résumé—L ‘auteur a détermine expérimentalement les diffusibilités dans lair des vapeurs 
d'acétone de tetrachlorure de carbone, et d'eau, en utilisant la méthode de Winker LMANN. Les 
valeurs obtenues dans certaines conditions d’opération ont une précision a ~ ou — 2%. Le méme 
appareillage sert a étudier lévaporation de mélanges liquides a deux constituants comprenant 
un produit volatil et un non volatil. 

Quand le produit non volatil est le plus dense, Pauteur montre qui le liquide reste entiérement 
a état de mélange, et la vitesse d’évaporation peut done étre calculée simplement comme pour 
un mélange idéal. Quand le corps non volatil est le moins dense, le mélangeage par convection est 
négligeable et le transfert de masse dans le liquide est un procédé de diffusion moléculaire en 
régime transitoire. 


Zusammenfassung Die Diffusionskoeflizienten’ von Azeton, Tetrachlorkohlenstoff und 
Wasserdampf in Luft werden experimentell in einer von WINKELMANN angegebenen Apparatur mit 
einer Genauigkeit von 4 2°, bestimmt. Die Verdampfung von flissigen Gemischen aus zwei 
Komponenten, einer flichtigen und einernichtflichtigen, wurden ebenfalls in der gleichen Apparatur 
untersucht. Wenn die nichtfliichtige Flissigkeit die gréssere Dichte hat, dann bleibt die Flissigkeit 
volistindig gemischt und die Verdampfungsgeschwindigkeit kann daher auf einfache Weise fiir 
ein ideales Gemisch berechnet werden. Wenn die nichtflichtige Komponente die geringere Dichte 
hat, kann man die konvektive Vermischung vernachlissigen und den Stofftransport in der 
Fliissigkeit als einen nichtstationiren Vorgang der molekularen Diffusion auffassen. 


INTRODUCTION will sink into the bulk of the liquid and considerable 
Wuerras there have been many investigations —™Xing will oceur, whereas in the latter, the residue 
relating to the evaporation of a pure liquid into a will tend to concentrate in the upper layers and 
gas stream, little attention has been paid to the the more volatile liquid must diffuse through this 
evaporation of a liquid consisting of two or more — TeS!0" before it can evaporate. 
components of differing volatilities. The more For the present study, two simple two-com- 
volatile components will then evaporate prefer- ponent systems have been selected. Each con- 
entially leaving a residual liquid which may be of — tained a virtually non-volatile liquid (dibutyl 
greater or lower density than the original liquid, phthalate) ; in the first mixture the volatile liquid 
according to the values of the densities of the was of lower density (acetone), and in the second 
components. In the former case the heavy residue — the volatile liquid was of higher density (carbon 


tetrachloride). The method of investigation which 
was adopted was to allow the mixture to evaporate 
at constant temperature in a vertical tube across 
the top of which was passed an air stream in 
which the vapour was carried away. The rate of 
evaporation was measured by observing the rate of 
fall of the liquid surface. 

A series of preliminary experiments was carried 
out using the pure volatile liquid, in order to 
examine the reproducibility of results and = the 
effect of experimental conditions. As a result of 
these experiments, values of the diffusivity of three 
vapours in air—acetone, carbon tetrachloride and 
water —have been obtained. 


APPARATUS AND EXPERIMENTAL 
METHOD 


The method adopted was based on that origin- 
ally used by Sreran [1] and [2], 
and subsequently modified by a number of workers 
including Lex and Wiikr [3]. The liquid to be 
evaporated was contained in the vertical limb of 
the tube shown in Fig. 1, the upper section A, 
serving merely for the introduction of the liquid. 


Al «Tube for filling 


( 


Liquid 


Tube diameters 
4,6 and 8mm 


+———25cm 


Fic. 1. Evaporation tube. 


The tube was immersed in a_ thermostatically 
controlled bath. Air from a cylinder of compressed 
gas passed through a copper coil immersed in the 
bath and its temperature was raised to that of the 
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bath. This air stream entered through the horizon- 
tal section B of the evaporation tube and the air 
leaving the tube was metered. The walls of the 
bath were of glass and the level of the liquid in the 
tube was measured by means of a cathetometer to 
within 0-002cm., At the beginning of an experi- 
ment the liquid was introduced into the tube and 
allowed to reach the temperature of the bath 
before the passage of air was commenced. Read- 
ings of liquid level were then obtained as a function 
of time. 

The effect of the operating variables was studied 
using acetone as the liquid. At a given tempera- 
ture, the two important variables were the rate of 
passage of air and the diameter of the tube. The 
air rate was varied between 0-5 and 2-1 cm*/sec 
and three tube diameters were used (4, 6 and 
8mm). The lower permissible limit of air rate is 
set by the necessity of maintaining an approxima- 
tely zero concentration of vapour at the mouth of 
the tube, and the upper limit by the production of 
excessive turbulence in the upper parts of the 
diffusing vapour. In tubes of too small diameter, 
the overall rate of diffusion will be affected by the 
velocity profile over the cross-section, and in 
excessively large tubes, turbulence is induced too 
readily. 


EvaPporation oF Pure Liquips 
Theory 


Under conditions of molecular diffusion, the rate 
of mass transfer per unit area of liquid surface is 
given by Stefan’s law as :— 


where h is the distance between the air stream and 
the liquid surface at time t, Pg,, Pg, are the partial 
pressures of air in the air stream and at the surface, 
p is the molar density of the liquid, D, is the 
diffusivity of the vapour in air, and C,, is the total 
molar concentration of gas and vapour in the 
diffusion region ; this will be constant and equal 
to the number of moles of gas plus vapour per 
unit volume in the tube. In a given experiment, 
h and ¢ are the only variables and therefore on 
integration : 


IL. 
59 
dh P 
N=p—= 
Pa (1) 


J. F. Ricnarpson 


(2) Thus, if plotted against h — hg, a straight 
0 


where hg is the initial value of h. line should be obtained. The diffusivity D, can 
Now the effective distance. hb. over which then be calculated from the slope s ; 

molecular diffusion is occurring cannot be estim- 

ated accurately because of end-effects ; the position — Le. D, = oC ~~. p (4) 
of the liquid surface cannot be defined precisely od =e 

because of the curvature of the meniscus, and at 
the upper end of the tube turbulence is induced in 
the diffusing vapour. On the other hand, values of 
h — hy can be measured very accurately and it is 
therefore convenient to rearrange equation (2) as: 


(h — hy) (h hy + 2hy) 2D.Cor In P Results 
t Pr 


and the effective initial distance over which 
diffusion is occuring can be calculated from the 
intercept on the ordinate axis. A typical curve is 
shown in Fig. 2. 


The effect of tube diameter and air tlow rate 
(h — he) was studied using acetone at 30°C. In each case, 
the results were corrected to standard atmos- 


t VOL. 


ie. = 
h—h,” 2D, Pm 


bo p hy (3) pheric pressure by multiplying by the factor 
Cop Pa Py, 
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Fic. 2. Experimental results. Acetone at 50-0 C in 6 mm tube. 


Table 1. Effect of tube diameter and air flow rates on apparent diffusivity 


Tube Air flow | Pressure 
diameter rate Temp. curve (s) dD, P 
(mm) (em? /sec) | Cc) (hr em?) (em? sec) (mm Hg) 


— 


300 0-957 O-105 
1-10 30-0 0-902 oli 
210 
O42 OAT O-105 
1-41 20-3 0-120 
O-830 
1-41 30-4 O-878 


t Pp, 
10 
| 
| 
§ 
c 
7 
| 
D, 
(em? sec) 
6 70 
767 
| 0-108 
752 | 
4 753 
8 74 
236 


It will be seen that the results do not vary by 


more than + 2 per cent from the mean value of 
0-112 em? sec, and it is therefore concluded that 
neither tube diameter nor air rate has a critical 
effect for the ranges investigated. The succeeding 
experiments were all carried out in the 6 mm tube, 
with an air flow of about 1 em*® sec. The variation 
of diffusivity Dr with temperature @ is seen in 
Fig. 3, and the diffusivities of acetone and carbon 


x 


| »="Carbon tetrachioride 


° 


0 20 30 40 50 60 70 
Temperature, 
Fic. 3. Diffusivities of acetone and carbon tetrachloride in 
air as function of temperature. 


tetrachloride vapours in air are given to within 
approximately 2 per cent by the relations : 
Acetone D, = 0-087 0-00086 (5) 
Carbon 
tetrachloride D, — 0-076 


0-00082 0 (6) 


Gusu [5] obtained a value of 0-111 em?,sec for 
the diffusivity of acetone for temperatures 
between 30 and 50°C, 
No values for the diffusivities are quoted in the 
literature, but approximate values have been 
calculated from Gilliland’s equation and compared 
below with those calculated from equations (5) 
and (6) at 42°C, the temperature at which the 
experiments with 2-component mixtures were 
carried out, 
Acetone tetrachloride 
From equations (5) and (6) 0-123 0-690 em*/sec 
From Gilliland’s equation 0-113 0-081 em?/sec 


An isolated measurement for water at 46 C gave 
a value of 0-258 cm? sec. 
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In a number of the experiments in the 6mm 
tube, the intercept on the ¢ (hk — hg) axis was 
measured, and the effective position of the top of 
the region of molecular diffusion was determined. 


Mixtures —Vonatite Liquip 


Less DeNsE 


Two-CoMPONENT 


The evaporation of a mixture of equal volumes 
of acetone and dibutyl phthalate in the 6 mm tube 
was studied at 42°C. Measurements were taken in 
the same way as for single component liquids. 
Dibutyl phthalate can be considered as entirely 
non-volatile (vapour pressure at 42°C is 
2x 10-4 mm Hg), and on the assumption that the 
mixture is ideal, the relation between the partial 
pressure at the interface and the concentration of 
acetone in the liquid at the interface is given by 
Raoult’s Law. 
denser than the evaporating liquid, it will be 
assumed that it mixes completely with the bulk of 
the liquid which will therefore be of uniform 
composition throughout at any given time, so that 
the partial pressure of acetone vapour at the inter- 
face, and hence the evaporation rate, can be cal- 
culated in terms of the amount of evaporation 
which has taken place up to that time. These 
various assumptions will be justified if the meas- 


As the residue is considerably 


ured and calculated amounts of evaporation 
correspond during the whole course of the experi- 
ment, The calculated amount of evaporation is 
obtained as follows. 

After time t, the number of moles of the more 
volatile component which have evaporated 

(h — ho) p, A. Thus, if n, and ny are the num- 
bers of moles of the volatile and non-volatile 
components initially, the mole fraction of the 

- (h — ho) p, A 

ny + —(h — hy) p, A’ 
The molar concentration of vapour at the surface 
of the liquid is therefore 


volatile liquid at time ¢ 


= n, — (h — hy) p, A (7) 


Ny + Ny — (h— hg) py A. 
where C,9 is the corresponding value for the pure 
volatile liquid at the same temperature. 

Now the rate of diffusion through the gas phase 
is strictly given by equation (1). However, the 


= 
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ratio P/P,,, was equal to 1-18 at the commence- 
ment of the experiment, and fell rapidly during the 
initial stages to a value of approximately unity. 
Since an algebraic relation between h and ¢ can be 
obtained only when the ratio is unity, this sim- 
plification was introduced, because it was con- 
sidered that the resulting errors would be less 
than those introduced by a graphical solution of 
the problem. 
Thus, from Fick's Law, 


dh ny 


(h — he) p, 


n, + my — (h A 


Rearranging, this gives 


dts p& 
dh DC, 


On integration : 


Ns 


Pr 


2D.C.6 
Ny ny 
DC AN\p.A 


(h? — h,?) (h — ho) — 


+ hy) in (1 — 


(h — (9) 


The experiments were carried out at 42°C with 
mixtures consisting initially of equal volumes of 
dibutyl phthalate and acetone. Under these 
conditions, the values of the quantities in equation 
(9) were as follows : 


O-764 
= = 0-0132 moles 


0-123 em® see 443 em? hr (from equa- 
tion 5) 

2 0-764 — 0-0263 ¢ moles. 

5s 
1-048 
27s 
A 0312 em? 
hg 

315 22,400 

== 2-32 x 10-5 moles/cm*. 


= 2 0-00753 moles, 


On substitution and simplification, the relation 


between ¢ and / (in hours and centimetres respect- 
ively) becomes : 
t — 2-40h 35-6 


— 41-8 logy, (7-56 —h) (10) 


— ° Calculated points 


t, mr 
Fic. 4. Experimental and calculated results for evapor- 


ation of equivolume mixtures of acetone and dibutyl 
phthalate at 42°C, 


In Fig. 4 is shown the experimentally deter- 
mined relationship between A and ¢ together with 
that calculated from equation (10). The agree- 
ment is seen to be within 4 per cent over the whole 
range, with the theoretical curve slightly below 
the experimental one, This small difference could 
be attributed to small inaccuracies in experiment 
or tothe assumption that P P,,, is equal to unity. 


Two-Component Liquip 
Drnser 

When the volatile liquid is denser than the 
non-volatile component, the residue after evapor- 
ation is of relatively low density and therefore 
shows little tendency to mix with the bulk of the 
liquid. Under these circumstances, therefore, 
diffusion must take place through the bulk of the 
liquid phase before further evaporation can occur. 
Since the capacity of the vapour phase is very 
much less than that of the liquid phase, it will be 
assumed that in the gas phase the rate of transfer 
is proportional to the overall concentration 
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gradient and that in the liquid phase the transfer is 
an unsteady state diffusional process. The liquid 
depth was assumed to be effectively infinite. 
Measurements of the rate of evaporation of a 
mixture of equal volumes of carbon tetrachloride 
and dibutyl phthalate in the 6 mm tube at 42°C, 
should then enable the liquid phase diffusivity, D, 
to be calculated. If the calculated value of D 
over the whole duration of the experiment remains 
constant the above assumptions concerning the 
mechanism of the process would appear to be 
justified. 

The rate of diffusion in the gas phase is given 
by : 
where C, is the molar concentration of vapour at 
the liquid surface at time ¢. If the mixture is ideal, 
Raoult’s Law will be applicable and the molar 
concentration of the volatile component in the 
liquid at the surface C, will be proportional to 
the concentration in the vapour. 


ie. Cy — (12) 


kp, h dh = bh (say) 


Thus Cy = D. dt 
v 


(13) 


For diffusion in the liquid : 
(14) 


where C is the molar concentration of the volatile 
liquid at a depth y below the surface. 

This equation must be solved to give the con- 
centration gradient at the surface for the following 
boundary conditions : 
t=0, 0<y< C=C, (uniform constant value) 

oY 
C=C, 
Cc == C, 
C=C, et 
.C,, are values of C, calculated 
from equation (13), using the experimentally 
determined values of h and dh dt at regular time 


intervals, This procedure was necessary because 
it was not found possible to express C, as a simple 
algebraic function of ¢. 

The solution of equation (14) is given in Appen- 
dix 2 as: 


(fort >t,) (15) 


At time ¢, the molar rate of evaporation per unit 
area is given by 


(16) 


Then from equations (13), (15) and (16), at time f¢, 


D (h dh dt), , —(hdh dt), 
4) 


h dh dt), —(h dh ‘dt dh 


From equations (12) and (13), 
k=C,/C, for pure carbon tetrachloride, 
and b = k p,/D,. 


Under the conditions of the experiment —-atmos- 
pheric pressure and 42°C-—-the vapour pressure 
of carbon tetrachloride is 251 mm Hg, so that 


1-54/155 
k=; 22,400 273,315 251 760 = 


2-4 hr em? 


[D, = 0-090 cm*/sec = 324 cm? /hr from equation 


(6)] 
Substituting in equation (17), 


dh 
D=2 (i) 


vt—t, 


(h dh dt), —(h dh 
(units of centimetres and hours) (tf > ¢,) 


+ (18) 
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hr 


Fig. 5. Experimental figures for the evaporation of equi-volume 
mixture of carbon tetrachloride and dibutyl at 42°C. 


(19) 


22 (say) 


Experimental measurements of a mixture of 
equal volumes of carbon tetrachloride and dibutyl! 
phthalate were made over a period of 4400 hr and 
values of o were calculated at various times from 
the experimental results, The results are given in 
Fig. 5 and Table 2. 


Table 2. Calculation of liquid phase diffusivity for 
carbon tetrachloride in dibutyl phthalate 


dh ‘dt e dh at 


Time, tf | 
(em hr) (em? 5) 


(hr) 


060 
127 
O-108 
120 
O-1038 
0-076 


0095 
0-060 
O22 
0-00085 
O-O0052 


138 
0585 
O-00824 
O-00768 


Over the period from 50 to 2600 hr, the value of 
(dh dt) o remained substantially constant. For 
lower values of time, the gradient dh/dt was 
changing very rapidly so that accurate measure- 
ments were diflicult to obtain. After more 
prolonged periods than about 2600 hr, the concen- 


tration was starting to change rapidly at the 
bottom of the tube and it was no longer justifiable 
to consider the liquid depth as infinite. A rough 
calculation has shown that the concentration at 
the bottom of the tube would have changed by 
about 20 per cent after 1000 hr. The mean value 
of (dh dt) o over the period was 0-113 hr? em. 
Hence the mean value for the diffusivity for the 
carbon tetrachloride — dibutyl phthalate system 
was given by equation (19) as 0-028 em*, hr. 


CONCLUSIONS 


The WiInkeLMANN method has been used with 
varying experimental conditions to determine the 
diffusivity of a vapour in air, It has been con- 
cluded that the results obtained are substantially 
independent of both air rate and tube diameter, 
over the range investigated. Values of diffusivities 
have been obtained for acetone, carbon tetra- 
chloride and water vapours to within an estimated 
accuracy of — 2 per cent. 

The evaporation of a mixture of acetone and 
dibutyl phthalate has been studied in the same 
apparatus and it has been shown that the heavy 
residue sinks and completely mixes with the 
remaining liquid. The levels of the liquid surface 
calculated on this assumption were found to agree 
closely with the experimental values throughout 
the course of the experiment. 

When a mixture of carbon tetrachloride and 
dibutyl phthalate evaporates, the evaporation rate 
falls off very rapidly with time because the residual 
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liquid is of a lower density and therefore remains 
close to the surface. The value of the liquid phase 
diffusivity has been calculated at various times 
during the process, on the assumption that the 
amount of mixing in the liquid due to convection 
is negligible and that the transfer process in the 
liquid approximates to that which would occur 
under conditions of molecular diffusion in a liquid 
of infinite depth. The value of the diffusivity so 
calculated was constant for a period of over 2500 hr, 
At the beginning of the experiment, accurate 
values were diflicult to obtain because of the rapid 
change in the evaporation rate. Towards the end 
of the experiment, low values of diffusivity were 
obtained because the concentration was changing 
appreciably in the liquid at the lower end of the 
tube. 

In this study, the behaviour of two simple ideal 
systems has been elucidated, Further investi- 
gations are needed in order to determine the be- 
haviour of the more complex mixtures which are 
normally encountered the petroleum and 
chemical industries, both under conditions of 
natural evaporation and during the combustion 


of the liquid in bulk. 
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Arrenpix—l 
Puysican Properties or 
Mareniats Usep 
Carbon tetrachloride 
Density (gy em*) 
1463255 — 10110 x 10°40 — 6-90 x 
Vapour pressure (P?, — mm Hg) is given by : 
1769 
logy, Py = TH875 — 
Acetone 
Density (g em*) 
= 081248 — 11x 10°39 — 8-58 x 
Vapour pressure (P?, — mm Hg) is given by : 
1660 
logy, P, = 7-930 — 
Dibutyl phthalate 
Density (g/cm®*) at 42. C = 1-048 


Arrenpix—2 

SoteTion or THe Liquip Puase Dirrustion 
EQuATION 


When ¢ 
Pute 


Then 


and when? = 0,¢ (O<y< @) 
ay 


When y 0, C= Cyande =e, att=t, 


When y 


c=C, and atl = ete. 
Mean value of ¢ over time interval to 4, = > 
a,, (say). 


Taking Laplace Transforms of both sides of equation (20), 


D (22) 
de 
When y = @, y = 0, and y 
B,=0 
when y = 0.¢ = By, 
d¢ 
and B,. yp D 
x 
Thus B, = dt 


ty 


dt ac e~Pt 4 


P 


a 
+ n-1 — Pn) +... 
Pp 


(20) 
| de ‘ 
eh! dt D wt dt 
(— + pe D 
Me op 
Le. dy? p* 0 
t 
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1 1 


a,,) 


1 
—e Mn 
+ (a 


1 


1 
+— 


1 1 


Substitution in equations (21) and (22) gives : 


— (C, — D 


1 
r ad "ap Cy 
nel 
ae 


and — = — —- 
2y pp 


= @ 


1 


d 
Taking the inverse transform [4] of equation (24), ° 


non 
net L2y 


(fort t,) (25) 
When ¢ = ¢,, the corresponding term is zero for all value 
of y. 
Thus, when y = 0, 
wy vit—t) vi | 


(fort (26) 


NOTATION 


A & cross-sectional area of evaporating surface 
a, mean value of ¢ over internal to 
B,, By = integration constants 
b = kp,/ D, 
C = molar concentration of volatile liquid at time ¢, 
depth y below surface 
uniform constant value of C att = 0 
= molar concentration of volatile liquid at surface 
value of C att = 
molar concentration of vapour at liquid surface 
molar concentration of vapour at surface of pure 
liquid 
total molar concentration of gas and vapour 
liquid phase diffusivity 
vapour phase diffusivity 
distance over which vapour is diffusing at time / 
initial value of h 
Henry's constant C, C, 
molar rate of transfer per unit area 
- number of moles of volatile component initially 
number of moles of non-volatile component 
initially 
total pressure 
= partial pressure of air 
= logarithmic mean partial pressure of air 
Vapour pressure 
parameter in Laplace transformation 
slope of curve in single component experiments 
absolute temperature 
time 
= distance below surface 
molar density of single component liquid 
molar density of more volatile component in 
liquid form 
= function of h, t, dh/dt as detined by equations 
(18) and (19) 
@ — temperature in degrees centigrade 
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Abstract-——Experimental frequency-response results and their interpretation by the diffusion 
equation lead to a correlation of liquid mixing with operating conditions for large-scale bubble-cap 
and sieve plates, in terms of a modified Péclet number and the liquid hold-up and froth density. 
Other mixing concepts are critically discussed. A simple method of predicting the effect of liquid 
mixing on plate efficiency is proposed, assuming complete mixing of the vapour between plates. 


Résumé Les résultats expérimentaux fréquence-réponse et leur interprétation par Péquation de 
diffusion conduisent a une relation entre le mélange liquide et les conditions d’opération pour les 
barbotteurs a champignons et les plateaux perforés & grande échelle en fonction du nombre de 
Péclet modifié, de la retention liquide et de la densité de la mousse. L’ auteur discute et critique 
d'autres concepts du mélangeage. Il propose une méthode simple pour prévoir Vinfluence du 
mélange liquide sur leflicacité du plateau, en supposant un mélange complet de la vapeur entre 
les plateaux. 


Zusammenfassung— Versuchsergebnisse des Frequenzganges und ihre Interpretation durch 
die Diffusionsgleichung fiihrt zu einer Bezichung zwischen der Durchmischung der Fliissigkeit 
und den Arbeitsbedingungen von grossen Glocken-und Siebbéden in Ausdriicken einer modifizierten 
Péclet-Zahl des fliissigen Inhalts und der Schaumdichte. Andere Vorstellungen iiber die Durch- 
mischung werden kritisch diskutiert. Um die Wirkung der Durchmischung der Fliissigkeit auf 
den Bodenwirkungsgrad vorauszusagen, wird eine cinfache Methode vorgeschlagen, bei der 
volistiindige Durchmischung des Dampfes zwischen den Boden angenommen ist. 


INTRODUCTION experimental determination until recently [4, 5). 
The need for such a parameter in a rigorous 
treatment of the non-steady state behaviour of 
distillation columns was the original reason for 
this investigation. In this paper a simple method 
is presented for characterizing liquid mixing on 
large plates in terms of an assumed eddy-diffusion 
mechanism, for predicting the value of the mixing 
parameter from the operating conditions, and for 
using this parameter to predict the true plate 
efliciency. Complete mixing of the vapour 
between plates is assumed. 


Avrnovcn methods have long been in use [1] 
for calculating the effect of vapour mixing on 
plate efficiency, the effect of longitudinal mixing 
in the liquid has been relatively neglected. In 
small laboratory columns perfect mixing can be 
assumed with negligible error; however, the 
effect of liquid mixing increases in importance 
with increasing path length, and only on very 
large plates is it justifiable to assume no mixing, as 
is usually done when calculating the vapour- 
mixing effect [1]. It has been shown [2] that on 
plates having a length of liquid path greater than 
about 5ft, the liquid-mixing effect tends to 
become the more important. Most previous work on mixing on plates has 

The existence of the liquid-mixing effect has aimed at a correlation of point and overall 
long been recognized [3], but the lack of a suitable | Murphree plate efliciencies in terms of a quantita- 
quantitative mixing parameter prevented its — tive mixing parameter. Plates have been regarded 
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as consisting of a number of perfectly mixed 
stages [2] which has been related to the point 
and overall efliciencies, but which has not been 
related experimentally to the plate size or the 
operating conditions, Data which are based on 
column overall-etliciency determinations and 
correlation with point efficiencies calculated from 
relative volatility [6], indicate that for liquid 
paths less than about 10 ft the equivalent number 
of perfectly mixed stages is roughly equal to the 
number of feet. Some step-response experiments 
performed with a single row of bubble caps [7], 
however, indicate that the number of stages is 
approximately equal to the number of caps — or 
rows of caps — traversed by the liquid. 

Other mixing parameters have been detined 
[8, 9] as functions of the concentration jump at 
the inlet weir, Their experimental determination, 
unfortunately, depends on concentration samples 
taken from the unstable “entrance length” 
of the plate close to the inlet weir; their use in 
relating overall plate efliciency to point efliciency 
is discussed by Crozier [10}. 

Dispersion on bubble-cap and sieve plates is 
more realistically treated as a continuous-flow 
phenomenon. In recent years considerable 
thought has been given to longitudinal dispersion 
in continuous-flow systems [11, 12, 13, 14, 15), 
specifically in pipes [16, 17], in fluidized beds [18}, 
in packed beds [19, 20, 21, 22, 23, 24, 25, 26) 
and in reactors [27, 28}. 

Several workers under GiLuiLanp at the 
Massachusetts Institute of Technology have 
supposed liquid mixing on bubble-cap and sieve 
plates to take place by an eddy-diffusion 
mechanism. — Longitudinal 
eflicients have been measured over a limited 
range of operating conditions [29] and a relation 
between point-and overall Murphree plate 
efficiencies has been developed [30] in terms of 
an eddy diffusivity. Similar relations have also 
been derived elsewhere [31, 32, 33). 

More recently, in a significant contribution by 
Foss [4, 5], liquid mixing and its effect on plate 
efliciency have been characterized in terms of 
liquid residence-time distribution functions based 
on the response to step changes in concentration 
of a tracer, The residence-time distribution 


eddy-diffusion —co- 
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approach yields the maximum useful information 
about continuous-flow systems, and can be used 
even when channelling and by-passing occur ; 
it is considerably more complicated than the 
assumption of an eddy-diffusion mechanism. If 
the mixing can be shown to obey the diffusion 
equation, this simplities both the experiments and 
the interpretation and subsequent use of the 
results. 

The diffusion model is likely to be adequate 
so long as the turbulent intensity is uniform and 
as long as no channelling or stagnation of liquid 
occurs on the plate. It would be expected to 
break down where the eddies are comparable in 
seale with the vessel in which they are formed. 
On a distillation plate of conventional design it 
can be assumed that the eddies are small compared 
with the vessel dimensions and that no concentra- 
tion gradient exists in a direction perpendicular 
to that of liquid flow. The effects of velocity 
and concentration gradients resulting from the 
changing cross-sectional area perpendicular to 
the direction of flow on circular plates are un- 
likely to be significant : the mean width differs 
very little from the width of the weirs. 

In the present work it was therefore decided 
to assume a diffusion model and to use frequency 
response, which, of the possible methods of 
dynamic analysis, is in some ways the most 
attractive [34,35]. [t has been shown by Kramers 
et al, {12} that a characteristic mixing parameter 
can be determined for a continuous-flow system 
by comparison of the experimental and theoretical 
frequency responses, provided it is valid to assume 
that the mixing process can be represented 
mathematically by a combination of piston flow 
and longitudinal diffusion. 


Frequency Rersronse 
From the unsteady-state material balance in 
time dt about an element of length ds of unit 
cross-sectional area, assuming a constant long- 
tudinal diffusivity D, it follows that ; 
Ne 
at (1) 
The transfer function of exit, with respect 
to inlet concentrations for this system, is derived 
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and its significance explained elsewhere 


If the concentration vectors measured are those 


upstream of the inlet (X;) and at a point distance s from the inlet (X,)* where 


s=aJ, 


the transfer function can be shown to become : 


X, expaz[p(2/a — 1) 


+ 1] + (p — +1) 


Xx, i[(p + exp 22p — (p— 


where p = 


and 


This was solved for attenuations with the aid 
of a digital computer, and the computed curves 
are shown in Figs. 2 and 3. The quantity “z,” 
which is a measure of the departure from perfect 
mixing, can be regarded as the ratio of the 
diffusion time-constant, or mixing time, to the 
mean residence time. It can be shown [12] that 
it is numerically approximately equal to the 
number of equivalent perfectly mixed stages 
(when this exceeds five) and is thus analogous 
to the number of transfer units. 

The dimensionless group ul D is a form of 
modified Péclet number. 


APPARATUS 

Several requirements were recognized : 

(a) The plates should be typical of full-scale 
plant, but a rectangular design facilitates inter- 
pretation of the results. 

(b) A tracer is necessary for following con- 
centration changes. It should be well dispersed 
on reaching the inlet weir, where uniform lateral 
distribution is required. It was decided to use 
electrical conductivity as an inferential measure 
of concentration because of its convenience and the 
ease with which it can be measured continuously. 
Sodium thiosulphate was selected as the tracer : 
its resistivity is nearly linear with concentration 
over the range used, and it was readily available. 

(c) The boundary conditions must correspond 
to those assumed in the theoretical treatment, 
but both upstream- and downstream measurements 
should be made using cells of identical design to 
enable measuring lags to be neglected, since these 


tri |} 


ul 


are not always small compared with the residence 
time of the liquid on the plate. This can be 
achieved by building both cells into the plate. 
The down-stream cell then replaces a cap or a 
number of perforations; the upstream cell, 
however, is required to sample the solution from 
above the inlet weir, with which it must therefore 
be supplied to the exclusion of the surrounding 
froth. 

(d) The sampling time-lag should be reduced 
to a minimum consistent with the necessity for 
deaeration prior to the continuous measurement 
of conductivity. This lag places a lower limit on 
the liquid residence-time on the plate, which 
should be greater than the cell time-constant. 
Conversely, as was found later, a lower limit of 
liquid-flow rate is reached when the output 
signal is masked to such an extent by “ noise ”’ 
as to be unrecognizable as a sine wave. The 
velocities of the large-scale eddies are then pre- 
sumed to have become comparable in magnitude 
with the mean liquid velocity, the flow being no 
longer predominantly in one direction but 
becoming increasingly random as the mean 
velocity tends to zero, 

Fig. 1 is the flow diagram. The plates are 
rectangular sections 1 ft wide and of 3 ft active 
length. Mains water mixes in the downcomer 
with a slow stream of concentrated sodium 
thiosulphate solution, the flow rate of which is 
modulated sinusoidally by applying a sinusoidally- 
changing air pressure through a valve positioner 
to a control valve. The amplitude, mean point 
and frequency of the oscillations are adjustable. 
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Fic. 1. Flow diagram. 


Stirrers in the inlet weir box reduce segregation 
to a minimum before the solution crosses the 
weir onto the plate. A small quantity is con- 
tinuously diverted from the weir through a cell 
containing two platinum electrodes. A similar 
conductivity cell is placed centrally at a convenient 
point downstream about three-quarters of the 
way towards the end of the active bed. Several 
sampling arrangements were tried before a 
satisfactory system was found. Both cells are 
connected to a fast and sensitive strip-chart 
galvanometer-recorder which also supplies a high- 
frequency excitation voltage to the electrodes, thus 
precluding polarization. 

The layout of both plates is conventional : 
7 rows, consisting alternately of 2 and 3 slotted 
caps, 3 in. diameter, mounted with a skirt height 
of ¥& in. on 2 in. diameter risers on an equilateral 
pitch of 4}in. The sieve plate has 41 rows 
consisting alternately of 14 and 15 holes ¥ in. 
diameter drilled in jin. thick perspex on an 
equilateral pitch of } in. Both plates are installed 
in a tank with perspex walls which enable visual 


tests to be made with traces of dye to ensure 
that all dead space has been eliminated over 
the range of liquid velocities used, and that 
interference of the cells with the flow pattern on 
the plate is negligible. 


EXPERIMENTAL Procedure anv ReEsuLts 


Experiments with both plates were done at water flow 
rates of approximately 80-320 ft® hr. ft width, superficial 
air velocities between 14 and 3 ft sec, and weir heights 
between 1 and 4in. For every combination of weir 
height and liquid rate a series of frequency responses was 
recorded for each of two gas rates. 

Owing chiefly to the large magnitude of the phase lags 
involved (many cycles long), it was found to be more 
accurate and more convenient to measure only the attenua- 
tions, the amplitude being taken to be the average value 
for a number of oscillations. The flow rates through the 
cells were equalized and adjusted to prevent air bubbles 
from being drawn down through the electrodes. 

It was found to be important to know accurately the 
residence time of the clear liquid on the plate. In a 
separate series of experiments this was determined for 
each set of operating conditions by simultaneously shutting 
off the water supply and blocking the exit weir, thus 
trapping the froth on the plate, and then collecting and 
measuring the volume of the clear water (still supported 
by the gas) by draining the plate through the holes left 
after removal of the samplers. Loss by spray and leakage 
were reduced to an acceptable minimum, Pressure probes 
were believed to be subject to large errors and were not 
used to measure clear-liquid depths. At the same time 
it was found that froth heights could also be measured 
with consistency and this was done for each combination 
of operating conditions. 

Calculations of mean liquid velocity and residence time 
are based on the volume of clear liquid on the plate, 
To show the dependence of this quantity on the gas-and 
liquid rates and the weir height (which for the purpose 
of this paper will be referred to hereafter as the operating 
conditions), an acration factor“, (defined here as the ratio 
of “h,” the clear liquid height after collapse of the foam, 
to “hy + h*,” the sum of the weir height and crest over 
the weir as deduced from the Francis formula) is plotted 
in Figs. 4 and 5 against liquid rate for different combina- 
tions of weir height and gas rate (expressed as a superficial 
air velocity). The ratio f, is in fact the ratio of the depth 
of liquid after collapse of the foam to the depth which the 
liquid would attain when flowing across the plate at the 
sume flow rate in the absence of aeration, 

The results derived from the bubble-cap and sieve plate 
runs are given in Tables 1 and 2. The value of the quantity 
“2” for a given run is that corresponding to the best 
curve through the experimental attenuations, which are 
superimposed for interpolation on graphs (Figs. 2 and 3) 
computed from equation (2), with “a” equal to 0-81, 
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The results are plotted as normalized attenuations, obtained 


by dividing the attenuation or amplitude ratio at any 
period by that at infinite period. 

The results are correlated as a plot of 2D 1° h against 
uh hy as proposed by Foss (4, 5) whose results, which also 
followed the diffusion law, were obtained by determining 
the residence-time distributions of liquid on sieve trays. 
For comparison, the best line through Foss’ results is 
reproduced alongside the author's results in Fig. 6. 


03 x 


2 
| 
| 


a 0-006 Best ine thraugh 


+—4 Yer 4 
+ + + + 7 
Foss’ results with 
soe ——$}- — 


+ + 
du= u, om/sec 
Fic. 6. Correlation of results. 


Discussion 


That the observed aeration factors for the 
sieve tray are lower than those for the bubble-cap 
tray, indicating greater aeration on the sieve tray, 
is attributed to more effective dispersal of the 
air on the sieve tray and also to the inclusion of 
the unaerated layer in the calculation of aeration 
on the bubble-cap tray. 
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2h 


The difference between the bubble-cap and 
sieve plate lines in Fig. 6 may be due to differences 
in turbulent) mechanism and and the 
stratification into aerated- and unaerated layers 
on the bubble-cap plate, resulting in higher mean 
froth densities. No obvious explanation suggests 
itself of the slight deviation of the plotted sieve- 
plate points from the best line through Foss’ 
results. The lines fit the equation; 


2D 
wh 


scale 


= €, (pu) (3) 


where ¢, » 0-9, ¢, ~ 2-4 for the author’s bubble- 
cap plate and ¢, » 0-5, c, ~ 3-0 for Foss’ sieve 
plate. “4” is the ratio h hy and will be referred 
to as the froth “* density.” 

Hence, over the range of operating conditions 
investigated, for the author’s bubble-cap plate : 


D h 24 
(4) 
and for Foss’ sieve plate 
D | hy . 
0-25 lh = (5) 


The group wh, D is another form of the Péclet 
number, with “” as the characteristic linear 
dimension. It is not possible to write complete 
dimensionless equations because the relevant 
physical properties of the liquid such as density 
and viscosity were not varied in either set of 
experiments. It may be noted that Foss found 
D to be independent of u, and that in the present 
experiments with the bubble-cap plate it varied 
only as the 0-6 power of u. The effect of gas 
velocity on the diffusivity appears to manifest 
itself only in so far as it affects the froth “density” 
(h hy). The underlying physical mechanism 
responsible for longitudinal mixing is not yet 
clearly understood. 


Errect oN Puate EFFIENCY 


It will be assumed that the Murphree point 
efliciency, detined by : 


is constant over the plate, that the composition 
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of vapour entering a plate is constant over the 
column cross-section, that the equilibrium curve 
is a straight line over the concentration range 
involved, i.e 


that the molar volume 1, p, where ¢ — pr, and the 
clear-liquid height, “4,” are constant over the 
plate, and that liquid mixing can be characterized 
by an eddy diffusivity. For convenience, the 
following equations are derived in terms of 
“h” and “p,” the height and density of clear 
liquid, rather than the height and density of the 
froth. The products hp’ and hy dp’, where 6 — h hy, 
represent liquid loading and are equal. 

Making the above assumptions, a steady-state 
material balance on an element of thickness ds 
at distance s from the inlet gives : 


er L de _ Gy ka _ & (y. — 
ds? hpD'ds hpD hp D 
The boundary conditions, the same as_ those 


assumed in deducing equation (2), are as follows 
[ll, 36): 


(6) 


By a mass balance over the whole plate: 


(9) 


Yar 


and by solving for the arbitrary constants in 
the general solution of equation (6) it can be 
shown [30, 32] that the ratio of the Murphree 
overall vapour-phase plate efficiency to the point 
elliciency is given by : 


Nplate _ 
4p" By 
= 


{(1 + p'Pexp —2(1 — p’) — 
p’ exp —2(1 +p’) — 4p’.} (10) 


where p= (14 


and aZ (the modified Péclet number) — 2z. 


y Me pont 


Fig. 7. Correlation of Murphree plate vs. point efliciences. 
((The graph is due to S. P. S. Andrews and R. W. Mackay, 
L.C.1. Ltd., Billingham Division). 


Fig. 7 was constructed from equation (10). 
When D 0, this expression for the efliciency 
ratio reduces to the Lewts [1] equation : 


“ptate By — 1. 
Bn 


(11) 


it reduces to unity when D 


The Murphree overall liquid-phase plate efliciency 
can similarly be related to point efliciency, which 
can be calculated from physical data and the 
operating conditions [30]. 

Once “8” is known, it should be possible, 
therefore, to predict the overall plate efficiency 
from Fig. 7 by the following procedure: du is 
estimated from the operating conditions and 
plate dimensions using the graphs of aeration 
factors of Figs. 2 and 3 to obtain “4” and hence 
“u.” The relation due to Crozier [10] : 


—Ing = O715 py + O45 (12) 


which implies dependence on gas rate only, 
can be used to predict an approximate value of 
“4.” u/D which determines the value of the 
parameter zZ in Fig. 7, is now obtained from Fig. 6 
for either a bubble-cap or sieve plate. This 
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graph is valid provided the operating conditions 
do not necessitate extrapolation to du > 12 or 
ou <i. 

No conclusive tests of this proposed calculation 
have been possible with data reported for ex- 
perimental plates. It is intended to apply it 
to large plates working under typical conditions 
with binary systems in which both vapour— 
and liquid-phase resistances are important. 


CONCLUSIONS 

(a) Assuming a diffusion mechanism, and 
using frequency response, a relation has been 
obtained from which the value of a mixing 
parameter (a modified Péclet number, which 
represents the ratio of the transfer rates due 
to bulk flow and to diffusion) can be predicted from 
the operating conditions on bubble-cap and 
sieve plates. The mechanism of the dispersion 
process is not understood, but longitudinal eddy- 
diffusivity has been correlated with liquid velocity, 
liquid hold-up and froth * density,” which have 
themselves been related to operating conditions 
over a range typical of industrial practice. 

(b) Gas rate has no significant effect on the 
mixing parameter distinct from its effect on froth 
density and on hold-up. 

(c) Values of longitudinal eddy-diffusivity 
vary between 25 and about 120 cm* sec, over the 
range of operating conditions investigated on both 
plates. 

(d) The overall Murphree plate efliciency for a 
binary mixture can be related to the point 
efliciency through the modified Peéclet number. 
The effect of longitudinal liquid mixing on the 
overall efficiency can therefore be calculated 
given the relevant physical properties of the 
components, the plate dimensions and the 
operating conditions. The effect increases in 
importance with increasing length of liquid path. 

(e) It is suggested that for the more common 
organic and inorganic distillates the density, 
viscosity and surface tension of the liquid and 
vapour are likely to be of only secondary impor- 
tance compared with the operating conditions 
and that the relations found here for water can 
be expected to hold for other liquids. 

(f) It is probable that, as Foss’ relation for 


sieve plates was found to be independent of hole 
spacing, so the relation found here will hold for 
different designs of bubble-cap plate. 
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NOTATION 


a=a3/Z 

¢ — deviation of a sinusoidal signal of tracer coneen- 

tration from its mean value (g mol/em?*) 

constants in the correlation of u, D, h and ¢ 

D — longitudinal eddy diffusivity in the liquid 
(em? sec) 


D, = eddy diffusivity in the vertical plane due to gas 
bubbles 
J, aeration factor 
h 
hy + h* 


G — vapour flow rate per unit area of plate 
(g mol see em? plate area) 
g — liquid flow rate per unit weir width 
hr. ft weir width) 
h = clear-liquid height (after collapse of froth) (em) 
=V/wlt 
hy — height of froth 
h,, = weir height (em) 
hh.” = weir height (in) 
h* — crest over weir (cm) 
= 00585 
i=y-1 
k = slope of vapour liquid equilibrium line 
L = liquid flow rate per unit weir width 
(g mol em weir width) 
P = forcing period (sec) 


= (: *) 


Q = liquid volume flow rate (em? /see) 
s = distance along plate from inlet weir or first cell 
(em) 
T — mean liquid residence time (sec) 
= 
t = time 
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mean liquid velocity (em ‘sec) » ~ Vapour-phase Murphree point efficiency 
= elear-liquid volume hold-up (em*) y* — 
superficial air velocity, based on total area . 0 
between weirs (ft) sec) 
= plate width (em) 


complex concentration vector 
mol fraction in the liquid of the more volatile p ~ molar density of unaerated binary mixture 


(g mol 


Yplate — OVerall vapour-phase Murphree plate efliciency 


Sa — ¥ 


component of a binary mixture 
- mol fraction in the vapour entering the plate Pg ~ Vapour density (Ib ft*) 

of the more volatile component 2 = total effective plate length (em) 
= mol fraction—in the vapour leaving the plate oe & variance of the distribution curve about its mean 


of the more volatile component (sec?) 


value of y for vapour in equilibrium with liquid 2Ds ui 
of composition h hy (air-water froth density 


=uZ/2D = aZ/2 w 27 = frequency (cs) 
D = « Diem-') =1/P 
= ratio of slopes of equilibrium and operating lines Suabacripts 


¢ = exit (when s = J) 
i ~ inlet (upstream of the weir, not when s = 0) 
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The performance of rotary concentric cylinder fractionating columns—1! 
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(Received 16 December 1958) 


Abstract—A type of laboratory fractionating column is described, consisting essentially of two 
concentric glass cylinders of which the inner one rotates. This apparatus is of general use in 
preparative work, giving an exceptionally low hold-up, pressure drop and H.E.T.P., and having 
no bearings or shaft seals in the working section. It also offers a convenient means of studying 
the general properties of film-type rectifying systems ; it functions as a wetted-wall column 
in which the vapour flow regime may be changed in a predictable way without alteration of 
throughput by changing the form and speed of the rotor. 

A technique is described for determining liquid and vapour side mass transfer resistances. In 
experiments at atmospheric pressure with a mixture of n-heptane and methyl cyclohexane it is shown 
that the separating power of the column is controlled by the vapour phase mass-transfer resistance, 
the liquid side resistance being negligible. Attempts are made to relate the overall column charac- 
teristics to the fluid mechanical properties of the vapour stream under various conditions of 
operation. 


Résumé—_Les auteurs décrivent une colonne a fractionner de laboratoire, formée essentiellement 
de deux cylindres de verre concentriques dont le cylindre intérieur tourne, Cet appareil peut 
étre utilisé dune fagon générale pour les travaux pratiques car il donne une faible rétention, 
des pertes de charge et un H.E.T.P. trés bas et qu'il n'a pas de palier ou joints de ’axe dans la 
section de travail. Hl permet deffectuer commodement l'étude des propriétés générales des 
systémes de rectification de type film. Il fonetionne comme une colonne a paroi mouillée dans 
laquelle le régime d’écoulement de la vapeur peut étre modifié de fagon prévisible sans altération 
de débit, en changeant la forme et la vitesse du rotor, Les auteurs donnent une technique de 
détermination des résistances de transfert de masse du liquide et de la vapeur. Pour les expériences 
a la pression atmosphérique avec un mélange de n-heptane et methyl cyclohexane, le pouvoir 
de séparation de la colonne est determine par la résistance de transfert de masse de la phase vapeur, 
la résistance du coté liquide étant négligeable. Les auteurs ont essayé de rattacher les 
caractéristiques globales de la colonne aux propriétés mécaniques du courant de vapeur sous 
diverses conditions d’opération. 


Zusammenfassung — Es wird der Typ einer Laboratoriums-Fraktionierkolonne beschrieben, die 
hauptsichlich aus zwei konzentrischen Glaszylindern besteht, von denen der innere rotiert. Dieser 
Apparat ist im allgemeinen Gebrauch bei priiparativen Arbeiten und gibt einen besonders geringen 
Flissigkeitsinhalt, Druckabfall und H.E.T.P.; er hat weder Lager noch Wellendichtungen im 
Arbeitsteil. Man kann mit ihm bequem das allgemeine Verhalten von Rektifiziersystemen yom 
Filmtyp studieren; er arbeitet als Kolonne mit benetzter Wand, bei der man die Dampfstromung 
in definierter Weise bei gleichem Durchsatz veriindern kann durch Veriinderung von Form und 
Geschwindigkeit des Rotors. 


Zur Bestimmung der Stoffiibergangswiderstiinde auf der Fliissigkeits- und Dampfseite wird 
eine Versuchstechnik beschrieben. In Versuchen bei Atmosphirendruck mit einer Mischung 
von n-Heptan und Methyleyelohexan wird gezeigt, dass die Trennwirkung der Kolonne durch 
den Cbertragungswiderstand auf der Dampfseite bestimmt wird, waihrend der flissigkeitsseit ize 
Widerstand vernachlissigt werden kann. Es wird versucht, zwischen dem gesamten Verhalten 
der Kolonne und den Daten der DampfstrOmung unter verschiedenen Arbeitsbedingungen eine 
Beziehung herzustellen. 


*Present address : Chemical Engineering Laboratories, University of Edinburgh and Heriot Watt College, Edinburgh 1. 
+Present address: Plessey Ltd., Northampton. 
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The performance of rotary concentric cylinder fractionating columns—1 


INTRODUCTION 

In fractionating columns of the type discussed 
here the vapours from the stillpot are made to pass 
through the annular space between two vertical 
concentric cylinders, of which the inner cylinder 
rotates while the outer one is stationary. Reflux 
from a condenser at the column top returns as a 
falling film down the stationary outer wall, 
from which heat loss is as far as possible prevented 
by lagging and compensating heaters. This kind 
of column is of interest for two main reasons. 
First, in comparison with other columns of 
similar dimensions and throughput, the H.E.T.P., 
pressure drop and hold-up may be made small, 
so that the rotary concentric cylinder column 
promises to be of great practical use, particularly 
for low pressure batch operation. Second, this 
type of column is essentially a wetted wall 
system into which an additional degree of freedom 
is introduced by the presence of the rotor, whose 
effects, unlike those of a spinning band, for 
example, are sufliciently determinate for theore- 
tical treatment. The rotor surface can be assumed 
to be dry under ordinary running conditions, for 
any reflux deposited on it from the condenser would 
quickly be thrown off by centrifugal force at the 
operating speeds used. It may therefore be 
supposed that mass transfer between the liquid 
and vapour phases in the rectifying section occurs 
at the relatively stable and well-defined surface of 
the film of reflux on the stationary outer wall of 
the annulus. Results obtained with such apparatus 
would thus be expected to throw light on some 
general aspects of the mass transfer processes 
occurring in wetted wall columns. 

The first published account of a rotary con- 
centric cylinder fractionating column was that of 
Witiincuam et al. [1], who described in 1947 an 
extensive series of trials with such an apparatus 
in the U.S. National Bureau of Standards. This 
machine was not designed for use as a practical 
still, but was well adapted for investigations of the 
rectifying properties of the rotary concentric 
cylinder system. In the course of this work it 
was found that the separating power increased 
as the rotor speed was raised ; beyond a certain 
critical speed of rotation the rise in separating 
power was very abrupt. It was supposed that 


this critical speed corresponds to the onset of 
turbulence in the vapour, resulting in increased 
radial diffusion. If the resistance to mass transfer 
between the liquid and vapour phases is predo- 
minantly in the vapour phase, it would be 
expected that the sudden development of vapour 
turbulence would explain the observed sharp 
improvement in column performance at a critical 
speed of rotation. 

It was found that the critical rotor speed 
corresponds to a Reynolds number for the vapour 
stream (based on annular thickness and mean 
vapour velocity) near that which would be 
expected to mark the transition between laminar 
and turbulent regimes for flow between flat 
plates ;* Le. assuming that the surfaces of the 
rotor and outer cylinder (radii approx. 3-7 em 
and 3-8 em respectively) can be treated as parallel 
planes. However, an attempt to calculate the 
H.E.T.P. of the column at various rotor speeds 
on this basis was not very successful [1]. 

Several papers describing rotary concentric 
cylinder fractionating columns bave since appeared 
({2}. [3]. [4], [5]) but no further serious attempt 
seems to have been made to elucidate the theory 
of operation. In particular, it has not been 
established to what extent the mass transfer 
resistance in the system studied by WiLLiIncHuam 
et al. resides in the vapour phase. Evidently this 
must be known before the influence of rotor speed 
on the performance of the column can be inter- 
preted quantitatively in terms of the effect of the 
rotor on the vapour stream. Furthermore, the 
effect of the rotor on the character of the vapour 
flow has not been made clear. Whatever may be 
the criterion determining the transition to 
turbulence, it does not seem to have been pointed 
out that the motion of the annular vapour 
stream probably differs from the assumed Couette 
flow in at least one important respect ; viz., that 
at rotor speeds well below that for the onset of 
general turbulence Taylor vortices [6] would be 
expected to be present. Inasmuch as these 


*From [1], the critical value of V, bp » found experi- 
mentally is about 3000, 

Wittincuam et al. give the observed critical Reynolds 
number as 2500; they define Re as 0-84 Vibe ° 
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regular vortices would be expected to promote 
radial diffusion in the vapour one might expect 
their appearance to be marked by an abrupt 
improvement in column performance like that 
now ascribed to the onset of turbulence. No such 
low-speed effect is reported the literature, 
however. There is thus some uncertainty about 
the relation between the fluid mechanics of these 
columns and their rectifying properties. In the 
present work an attempt has been made to throw 
some light on this. 

The potential practical value of rotary concen- 
tric cylinder fractionating columns is evidently 
very great. The apparatus described — by 
WiLunGcuam ef al. was shown to give (at total 
reflux) an efficiency factor, defined [1] as; 

No. of plates per unit length x Throughput 


Hold-up per unit length 


of the order 10,000 plates hr, several times greater 
than that for any type of column previously 
reported. In spite of this, little progress seems to 
have been made in adapting the pioneer machine 
for general preparative purposes, The original 
apparatus and all its variants so far described 
have a rotor bearing at the lower end of the 
working section, immediately above the still pot, 
fully exposed to the hot vapour stream and the 
returning reflux. Even when this is a_ plain 
journal bearing of inert material, as in’ the 
column of and Bronce [2], its presence 
introduces lubrication difliculties at the relatively 
high rotor speeds required, and limits the range of 
substances which the column can handle. 

In the work to be described, two rotary concen- 
tric cylinder fractionating columns were employed. 
One of these was designed primarily for experi- 
mental purposes; the other, a larger machine, 
was designed later for general preparative work 
on a laboratory scale. This apparatus has an all 
glass working section containing no exposed 
bearings or shaft seals, and has proved entirely 
suitable for general use as a laboratory column of 
exceptionally high separating power, low hold-up 
and low pressure drop. 


EXPERIMENTAL 


Figures 1 and 2(a) show the general arrange- 


ment of the larger column, which will be referred 
to as Column 2. The rotor, R, is a Chance 
“ Veridia ” tube, specially selected by the makers 
for straightness and exterior uniformity, 2.6 em 
nominal outside diameter and approximately 
7T4em long. It is driven as shown in Fig. 2 (b) 
by means of a B.T.H. “ canned” squirrel cage 
motor, the squirrel cage (an aluminium casting 
with copper inserts) being actuated through the 


+A. 


Fic. 1. Column 2. General mechanical arrangement. 
Schematic. Horizontal scale much distorted and enlarged. 
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thin stainless steel wall of the stationary housing, 
C, by the field of an external stator, F, supplied 
from the a.c. mains. The stainless steel housing is 
connected to the column proper by means of a 
vacuum tight system of flanges and O-rings, so 
that the column may be operated at reduced 
pressure. 

The glass rotor and squirrel cage turn about a 
central fixed shaft, A, which extends downwards 
into the hollow closed rotor and is rigidly secured 
at its upper end in a massive wall clamp. The 
ball races of the squirrel cage serve as the upper 
bearings of the rotor; as these are located above 
the spiral water condenser they are relatively 
inaccessible to the column vapours. 

The lower bearing, Z, is a ball race mounted on 
the stationary inner shaft, A, at a point determined 
by calculation, about two-thirds the length of 
the shaft below the wall-clamp. As the lower end 
of the tubular rotor is closed, this bearing is not 
exposed to the column vapours. 

This arrangement in which the cylindrical rotor turns 
about an internal bearing on a fixed shaft overhung or 
cantilevered from the wall clamp above has, of course, a 
critical speed at which the motion becomes unstable. 
Fairly extensive preliminary work was carried out to 
determine the dynamical characteristics of this arrange- 
ment of rotor and cantilevered shaft ; using data obtained 
in these investigations the apparatus was designed so that 
the normal operating speed of the motor, 1470 rev/min, 
is far above the critical, and the rotor runs without 
perceptible vibration. 

The annular clearance between the rotor and the fixed 
precision-bore outer tube is about 0-1 em and the annular 
rectifying section is 61 em long. Since there are no bearings 
in the annular space to maintain this relatively small 
clearance between rotor and outer tube, or stator, these 
parts must be exactly located in their proper relation 
by the wall mountings. These demountable fixtures were 
made in accordance with Kinematic Design principles [8], 
and have the property that, once aligned, they allow of the 
ready dismantling and reassembly of the column without 
significant change of rotor-stator clearance. The mountings 
are also designed to allow of longitudinal thermal expan- 
sion of the column jacket without distortion or change of 
alignment. 

Near adiabatic operation, important for high perfor- 
mance wetted wall columns in general, is ensured by the 
following automatic system. The rectifying section is 
enclosed by a tubular copper jacket, surrounding which are 
three sets of electrical heaters, W. The electrical input to 
these is controlled by relays operated by the out of balance 
current in a Wheatstone network, of which the ratio arms 
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are connected in sequence by a motor-driven rotary switch 
with each of three pairs of matched thermistors mounted 
on the outer glass tube of the column and its copper jacket. 
When one of the thermistor pairs is unbalanced the current 
in the appropriate jacket heater is increased, Using this 
system the temperature difference across the air space 
between the column and jacket is found not to exceed 
02 C, as measured by thermocouples located at the mid- 
points of the three heated sections, 

The column is designed for intermittent take-off, and 
has «a magnetically operated condensate valve, 5S, 
immediately below the condenser. In ordinary use the 
valve solenoid is energized periodically by a process timer. 
In the experiments to be described, however, the column 
was operated at total reflux and the valve was kept 
closed except when a condensate sample was required for 
analysis. 

Still pot samples were withdrawn through a syphon of 
fine-bore glass tubing fitted by means of a ground glass 
joint to a side arm of the still pot flask, 

The boil up rate was controlled manually by means of a 
variable transformer in the circuit of the still pot heating 
mantle. A boil-up rate meter [9] was interposed between 
the still pot and the column base, with which the flow rate 
of returning reflux could be measured to approximately 
+ 25 per cent. 

The arrangements in the experiments with the smaller 
apparatus, Column 1, were substantially the same. The 
column dimensions were: rotor diameter 1-6.em, annular 
clearance O-Lem, length of rectifying section, 60.em*. 
This column was so designed that the rotor speed could be 
varied readily ; the rotor was driven by a } h.p. induction 
motor through a variable gear. The cantilevered shaft 
for this early column was designed before the mechanical 
data used in the construction of Column 2 were obtained, 
and it was found in operation that the desired range of 
working speeds lay close to the critical for the rotor- 
cantilever system. An additional constraint or guard 
bearing was therefore provided below the rotor at the base 
of the working section of the column. This took the form 
of a Teflon bush supported and centered in the outer 
tube by a glass spider ground to fit the bore ; in a con- 
centric hole drilled in the bush ran a short pivot of glass 
rod about 0-3 cm diameter, fused to the lower end of the 
rotor and ground concentrie with the external surface. 
The annular clearance between pivot and bush was of the 
order 0-01-0-02 em. This guard bearing or check ring, 
used in conjunction with the internal cantilevered bearing, 
gave little trouble at speeds up to 3000 rev min, though 
exposed to the hot vapours and reflux in the rectifying 
section. 

In experiments with Column 1 the rotor speed was 
measured by means of a stroboscopic lamp and checked 
by means of a mechanical tachometer having a friction 


*A 45 em rectifying section was used in some of the work 
at higher rotor speeds. 
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wheel which could be applied directly to the surface of the 
upper part of the rotor. For Column 2 the stroboscopic 
lamp alone was used. 


In all experiments the test mixture contained 
approximately 20 mole per cent n-heptane, 
80 mole per cent methyl cyclohexane. The 
component hydrocarbons, of the purest commer- 
cial grade, had the refractive index values given 
in [10]. Head and still pot samples from the 
column were analysed refractometrically at 20°C. 
with an Abbé refractometer having an accuracy 
of 2-3 units in the fourth decimal place. This 
allowed plate values to be determined with an 
error not exceeding 4 per cent under the operating 
conditions used. 

The columns usually reached thermal equili- 
brium, as indicated by the column and jacket 
thermocouples, 4-5 hr after start up; after this, 
material equilibrium was assumed to have been 
reached when two consecutive head samples, one 
taken not less than } hr after the other, had the 
same refractive index. The plate values of the 
columns were obtained from these equilibrium 
refractive index readings by construction from a 
large-scale diagram prepared from the data of 
Bromitey and QuIGGLE [10}. 


Resutts anp Discussion 


Column 1 gave some trouble, particularly at 
the higher speeds, due to bowing of the rotor, 
and its performance was sometimes erratic. 
Results obtained with it are recorded in Table 1. 

Column 2, the later of the two machines to be 
built, proved to be remarkably reliable and trouble 


free in operation. No source of variable frequency 
a.c, was available, and it was therefore operated 
in all experiments at the fixed design speed, 
1470 rev, min (or at zero rev min). Curve A of 
Fig. 3 gives its performance under atmospheric 
pressure at various throughputs at total reflux. 
It can be seen that at the lowest boil up rate 
used, viz. 130ml hr, the column gives about 
57 plates in its 61cm working length. This 


60) 


LENGTH 


PLATES / 61cm 


“360 400 500 600 700 S00 
BOIL-UP RATE, ml /hr 
Fic. 3. A’ Westhaver theory for glass rotor 
A (©) Glass rotor (0-12 cm annulus) 
B (©) Copper rotor (0-33 em annulus) 
B (A) Copper rotor with sweeping wires 
C (9) Copper retor with nichrome spirals 
D (>) Copper rotor with wire fins (16 5.W.G.) 
Rotor speed = 1470 rev/min, 


Table 1. No. of theoretical plates for column 1. 


Boil-up rate 
(ml hr) 
0 


oo 
125 
18-0 
20:5 
2-5 


Rotor speed (rev min) 


2.900 


1,500 2.500 


50 

Ino 
240 
27°55 


140 
17-5 
20-0 
56-0 


25-0 


Length of rectifying section 45cm. Diameter of inner cylinder 1-6 em. 


Diameter of outer cylinder 1-5 em, 
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separating power is greatly superior to that to be 
expected from, e.g., a spinning band column of 
the same length and overall diameter working 
under the same conditions, The hold-up, measured 
by allowing the apparatus to come to equili- 
brium, removing the still pot and allowing the 
column to drain into a graduated vessel, is not 
greater than 1-5 ml over the whole boil-up range 
investigated. At 130ml hr therefore, the 
efliciency factor [1] referred to above, becomes : — 

61 15 
This value compares very favourably with that for 
columns of other types, but is only half the best 
value obtained by Witiincuam ef al. [1] for their 
rotary concentric cylinder machine, probably 
(as will be shown below) because the rotor speed 
here is not great enough to excite turbulence in 
the vapour. 

As the working section of the apparatus is 
entirely of glass, the very favourable characteristics 
of this column should recommend it for general 
laboratory purposes. 


Liquid and vapour phase resistances 


Experiments were undertaken with modified 
versions of Column 2 to determine the relative 
significance of the liquid and vapour phase mass 
transfer resistances. In these, the glass rotor was 
first replaced by a copper one of smaller diameter, 
(annular clearance 0-33 em), and the separating 
power of the new system was measured for 
various boil-up rates at the normal rotor speed, 
These measurements are recorded on curve B 
of Fig. 3. 

Next the copper rotor was furnished with two 
“ wipers,” straight lengths of fine piano wire, 
(0-035 em diameter), parallel to the axis, attached 
to the rotor at either end by spring strips so 
arranged as to cause the wires to bear against 
the inner surface of the fixed outer tube of the 
column, It was found that when the rotor revolved 
this simple arrangement quickly dislodged a fine 
layer of tale dusted on to the inside surface of the 
outer tube, and it was therefore concluded that 
the wires would disturb a film of reflux there. 
The separating power of the column fitted with 
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this rotor was now measured at various boil-up 
rates, as before. The results recorded in Fig. 3 
fall almost exactly on curve B, the curve obtained 
with the bare copper rotor. Apparently the 
stirring of the reflux film does not significantly 
improve the column efliciency. It thus appears 
that the liquid phase resistance in this system is 
negligible. 

Finally, the copper rotor was provided with two 
longitudinal fins of 16 gauge copper wire soldered 
to the surface along their length. These did not 
touch the surface of the film of reflux at the outer 
wall but did presumably stir the vapour stream 
when the rotor revolved. From the results 
recorded on curve D it is evident that their 
presence considerably enhances the separating 
power of the column. Taken together with the 
previous result, this confirms that the mass 
transfer resistance here resides in the vapour 
phase. 

Further confirmation of this conclusion is 
provided by the following result. Two long finely 
wound nichrome helices were attached to the 
rotor in such a way that when the latter revolved 
the wire coils touched the interior surface of the 
outer glass tube. Such an arrangement has been 
recommended [11] as constituting an exceptionally 
effective form of “spinning band.” In the 
present experiments, the wire helices are found to 
improve the performance of the bare rotor but 
to give poorer results than the solid fins of 16 
gauge wire, as appears from curve C of Fig, 3. 
This would be expected on the basis of the 
previous conclusion, since, though they disturb 
the liquid film as the solid fins do not, the fine wire 
coils are probably less effective in the important 
function of agitating the vapour stream. 

These results offer direct evidence that resistance 
to radial diffusion within the liquid film has 
no signiicant effect on column performance in 
the present circumstances*. This conclusion 

*It is possible that the sweeping wires cause some 
longitudinal mixing in the reflux film, such as might mask 
an improvement in columa performance brought about 
by the improved radial mixing in the liquid ; but since the 
net effect remains undetectable over a wide range of 
experimental conditions it may be concluded that both 
liquid mixing effects are of very small importance. 
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bears on a matter of controversy. Thus Narrapze 
[12] concludes from observations on annular and 
open tube distilling columns that the liquid phase 
resistance makes a considerable contribution to the 
total resistance to mass transfer between the 
phases. Peck and Wacnen [13], on the other 
hand, conclude that the liquid phase diffusion 
resistance is in general negligible in wetted wall 
columns, a view which is supported by the work of 
Zuiperwes [15], Maryusov et al. [16] and others, 

However in view of the fact that the work of 
Cuani and Srorrow [14] suggests that the relative 
importance of liquid and vapour resistances in 
rectification is a function of the nature and 
composition of the test mixture, the general 
question cannot be decided on the result of a 
single experiment. But it appears that the 
experimental technique described is of general 
application in the study of liquid and vapour 
mass transfer resistances, and it is hoped to 
extend these enquiries to other systems. 


Significance of Taylor vortices 


In the present case it is now clear that the 
observed performance of the Columns 1 and 2 
can be related to the conditions in the vapour 
stream alone. Thus, the improvement in separa- 
ting power obtained by increasing the rotor speed, 
shown for Column 1 in Table 1, is attributable to 
the increase in radial mixing in the vapour stream, 
and is not (for example) to any significant extent 
due to the disturbance of the film of reflux by 
the whirling vapour. It would therefore be 
expected that observed changes in the overall 
performance of these columns would reflect 
changes in the radial diffusion coeflicient in the 
vapour stream; and consequently, that par- 
ticularly marked changes in column performance 
would accompany changes in the vapour flow 
regime. 

In view of this it is somewhat surprising that 
the performance of such columns is apparently 
little improved at rotor speeds at which Taylor 
vortices would be expected to appear in the 
vapour current. Such vortices [6], oceupying 
toroidal regions coaxial with the annulus, cause 
circulation of the fluid between the inner and 
outer cylinders in the plane of their common 


axis ; this motion has been shown [7], [20] greatly 
to enhance the transfer of heat from a cylindrical 
wall into an annular gas stream, and it would 
be expected to promote radial mass transfer also. 
However, it appears from the results given in 
Table I that the column performance is relatively 
little improved at the critical speed for vortex 
formation in a stagnant fluid. Using Tayior’s 
approximate equations 5.48 and 7.11 of [6], and 
assuming the value for the kinematic viscosity 
of the column vapour given in [15], it would be 
expected that vortices would appear in the annulus 
of Column 1 at about 350 rev, min; Table 1 
shows that although there is a rise in the separating 
power of the column (from 9 to 15 plates) between 
rotor speeds of 0 and 410 rev ‘min, there is no 
further significant change in performance as the 
rotor speed is increased.* Even the small improve- 
ment observed might be accounted for, not by 
the effect of vortices, but by the improved 
distribution of reflux obtained when the rotor is 
driven fast enough to throw liquid from_ its 
surface onto the outer wall of the annulus. 
Similarly, according to Tayior’s equations, 
vortices would be expected to appear at a rotor 
speed of approximately 130 rev ‘min in the column 
of Wituincua™ eft al. [1]. These authors made no 
measurements near this speed ; but between 0 and 
1750 rev min their results show an increase of 
separating power of a few plates only. In view of 
the results given in Table 1 it seems possible that 
this increase does not occur progressively (as is 
assumed in [1]), but abruptly, near the Taylor 
critical speed. But in any case this is an insignifi- 
cant change compared with that observed at 
higher rotor speeds, attributable to the onset of 
general turbulence. Again, for Column 2 furnished 
with the copper rotor the appearance of the 
vortices would be expected at about 50 rev, min. 
The fact that the attachment of fins to the rotor 
materially improves the performance at 1470 
rev min indicates that the radial mixing induced 
by the fins is a substantial improvement on that 


*It is interesting to note from Fig. 4 that the breakdown 
of Couette streamline flow in this system at a rotor speed 
of 410 rev/min would be in accordance with the results 
of Face [17], but not with those of Cornisu [18] or Kaye 
and ExGar [7]. 
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which obtains when the smooth rotor operates 
well above the Taylor critical speed. 

These facts suggest, cither (a) that Taylor 
vortices are not formed at or near the Taylor 
critical speed in these columns ; or (b) that they 
are not effective in promoting mass transfer from 
the outer wall of the annulus to the vapour 
stream ; or (c) that they do enhance such radial 
mass transfer but at the same time give rise to 
other effects which counter the expected improve- 
ment in column performance. 

There is no direct evidence for the presence of 
vortices in these experiments; and Tayior’s 
criterion for their formation, quoted above, 
applies to systems in which there is no axial flow, 
not strictly to the present case. However, the 
occurrence of Taylor vortices in annuli in which 
there is axial flow is apparently established by 
the work of Kaye and Excar [7]. It is known 
that the critical rotor speed for the change of 
regime from purely laminar flow is higher when 
there is axial flow than when the fluid is stagnant, 
The precise criterion for this transition in the case 
of simultaneous axial and circumferential flow is 
still in doubt, no successful theoretical treatment 
of the problem having yet appeared and some of 
the experimental investigations having given 
somewhat conflicting results [17], [18]; but on 
the basis of what is known at present it appears 
that at the axial vapour velocities used in the 
distillation experiments the critical speed for 
the development of vortices should differ very 
little from that calculated from TayLor’s equa- 
tions in all the cases considered. Thus, unless some 
feature peculiar to fractionating column operation 
inhibits vortex formation, it is reasonable to 
suppose that vortices are present at the low rotor 
speeds given above. It then appears from the 
results quoted that they may have an effect on 
column performance manifested at or near the 
critical speed for their formation, but which does 
not then increase significantly as the rotor speed 
is raised ; and that this effect is in any case small 
compared with that of the change (presumably 
the development of turbulence) found to occur 
[1] in the vapour at relatively high rotor speeds. 

It is thought that these conclusions, which must 
remain tentative until the existence of Taylor-type 


vortices under column conditions is definitely 
established, might be reconciled with the observa- 
tion [7] that the vortices are effective in promoting 
radial heat transfer in the following way. 
Kramers et al. [19] have shown that the Taylor 
vortices are strikingly effective in promoting 
longitudinal, as well as radial, mixing. In the 
heat transfer experiments described in [7] only 
radial mixing effects were involved ; the overall 
performance of wetted-wall fractionating columns, 
however, is sensitive to both radial and longitu- 
dinal diffusion effects. If both increase together, 
the improvement in column performance to be 
expected from improved radial mixing will be to 
some extent offset by increased back-mixing, 
which impairs the counter-current character of 
the vapour-liquid contacting process. Studies of 
these radial and longitudinal diffusion effects in 
annuli, designed to throw light on this matter, 
are at present in progress. 


Conditions for turbulence 

The signiticance of the high critical rotor speed 
observed by Witiincuam et al. [1] may now be 
considered, It appears from the experimental 
work of Kaye and Eicar [7] that the annular 
vapour stream in these columns could be in any 
one of four flow regimes under different condi- 
tions of axial and circumferential motion. The 
stability limits for these regimes are represented 
by the curves bounding the lettered areas in 
Fig. 4. which relates to the flow of air in an 
annulus, and is reproduced (with additions) 
from [7]. Of these four, it seems that three 
regimes are of account under ordinary conditions 
of column operation, viz. laminar, (/), laminar 
+ vortex, (1 + v), and turbulent + vortex, 
(t + v). The fourth possible flow regime described 
by Kaye and Excar, in which the fluid is purely 
turbulent with no regular vortices present, is 
apparently not stable in the range of axial flow 
rates used in the columns described here. It thus 
appears probable that the critical rotor speed 
observed by Wittincuam et al. at about 2300 
rev ‘min corresponds to the transition between the 
vapour flow regimes indicated in Fig. 4 as 
“laminar + vortex” and “turbulent + vortex.” 

Unfortunately the criterion for this transition 
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stream. Curves of Kaye and Excar, [7]. 


A. Laminar flow. 
B. Turbulent flow. 


I. Column I. 410 rev/min, 450 ml. hr-! reflux, 6/r,, 
IL. Column of ef al. 1500 ml. hero! reflux, 


Il. 2,000 rev min: no turbulence. 


Vand V. Column 2. 
IV. 430ml. hero? reflux, 0-26. 
for the column of Wituincuam et al. cannot be 
obtained from the results of Kaye and Excar. 
The latter workers used two annuli in their 
experiments, having 6/r,, values very much 
larger than that for the column of Wituincuam 
et al. (0-198 and 0-307 as compared with 0-03 in 
Witiincuam’s work). Although the diagrams 
given in [7] might be supposed to include the 
effect of b r,, changes, and the authors do in 
fact use them in attempts to correlate results 
obtained at other b/r,, values, it seems probable 
that the curve representing the (/ + v)/(t + v) 
boundary should be plotted with b/r,, as a 
parameter, as indicated in Fig. 4. This appears 
likely from the following argument. Six variables 
are here involved in specifying the flow regime, 
viz. w, p and ». The largest number of 
these that will not form a dimensionless product 
is three ; hence, by the /7 theorem, three dimen- 
sionless groups are required to characterize the 
flow, not two only as used by Kayr and Excar. 
Otherwise expressed, the dimensionless groups 
used by these authors are (1) a Reynolds number 
relating to axial flow, (2) the dimensionless 


(copper rotor). 


C. Turbulent flow + vortices. 

D. Laminar flow + vortices. 

O11. Vortices. 

= 0-027 
2.400 rev min: turbulence. 

1470 rev, min. 

120 ml. hr-! reflux, b/r,, = 0-27. 


Il. 


product of three groups, viz. a Reynolds number 
relating to the circumferential flow, the group 
(b/r,,)'*, and a geometrical factor, 1/Fg. The 
value of this product (referred to in [7] as the 
Modified Taylor number) can be shown from 
TayLor’s work to characterize the transition 
from laminar to vortex motion, independently of 
b/r,, in systems where there is no axial flow ; 
but there is no justification for supposing that it 
has a like significance in characterizing the 
transition to turbulence, in determining which 
the quantity 4 /7,, is probably involved in some 
functionally different way. It seems therefore 
that Kaye and Exvcar’s curves can be used only 
to correlate results obtained in experiments at 
the b/r,, values for which the curves are plotted. 
The differences between their published curves 
obtained at two different values of this ratio 
bear this out. 

It is not therefore surprising that Fig. 4 does 
not fit the data of WiLuincuas et al., indicating for 
their column a strong dependence of critical rotor 
speed on axial vapour velocity which the experi- 
mental results [1] do not support. Points on Fig. 4, 
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representing conditions in the column of 
Witiincuam et al. at rotor speeds of 2,000 and 
2,400 rev min and a boil-up rate of 1,500 ml hr. 
lie close to the (extrapolated) nearly horizontal 
right hand branch of the curve marking the 
(1 +v) (t + v) transition in Kaye and Excar’s 
system. This implies a critical dependence of 
vapour regime on axial vapour rate, and indicates 
that the rotor speed at which the vapour becomes 
turbulent should depend very markedly on the 
axial velocity. But Fig. 7 of [1] shows that 
conditions in the column are in fact little in- 
fluenced by boil-up rate; it appears from this 
figure that the vapour becomes turbulent at a 
relatively well-defined and constant rotor speed, 
irrespective of boil-up rate. Evidently for a value 
of br, as small as that for the column of 
Wituincuam et al. the curve representing the 
boundary of the ¢ + 0 region is more steeply 
inclined to the axis at low Reynolds number 
than those shown. Indeed where the parameter 
b r,, becomes very small the following considera- 
indicate that the bounding the 
region may approximate to elliptical ares 


tions curves 


having the major and minor axes on the axes of 
co-ordinates, and so meeting the horizontal axis 


at right angles. 


This result follows if it is assumed that when 6 1, is 
very small the eriterion for the 0) v) transition 
(4c. the eritical Reynolds number based on mean vapour 
Velocity) is nearly the same as that for the Laminar-turbulent 
transition for pressure flow between stationary tlat parallel 
planes. Since the 
shearing flow (for which the turbulence criterion is not 
known) this is clearly a questionable assumption ; but 
for zero axial velocity it does lead to the turbulence 
criterion used in [1], for whose validity there seems to be 
experimental evidence from at least one other source [3). 
Where there is both axial and circumferential flow, the 
transition from one regime to the other would on this 
basis be expected to occur at some critical value of a 
Reynolds number based on the Pythagorean sum of mean 
axial and circumferential vapour velocities. 


motion due to rotation is in fact a 


The groups plotted in Fig. 4 may be written ; 
V. 2b Pp” ! ry 


Then if 
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and KV. Ke X 
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represent a corresponding pair of critical axial and circum- 
ferential vapour velocities for the change of regime from 
vtot+ v, we have: 


Ke 12h en 1 Re, 
2 "2 2 
y x bp 


AK depends on the radial velocity gradient, and may be 
taken as approximately constant in (1), at least where 
the axial flow rate is not large. In that case, ie. near the 
vanxis, N and Y and the quantities on the right-hand 
side of (1) will be constant for a given thuid and a given 
pair of evlinders, Under these conditions (1) is the equation 
of an ellipse having the major and minor axes on the axes 
of and 

From (1), the length of the semi-minor axis is 2Re,. 
The intercept on the c-axis is 


A Fe Um 


which reduces to 


Re, 
kK 
when 6 7,, is very small, Thus for this limiting case, on 
the assumptions given, the (1 > v) (f+ 0) boundary curves 
on Fig. 4 should be a family of elliptical ares having a 


common minor axis, and semi-major axes (a-intereepts) 
proportional to yb 


It would be expected from the above treatment 
that for columns having very small values of 
br, the simple turbulence criterion proposed in 
[1] would be, in effect, valid ; for under ordinary 
operating conditions the axial vapour velocity 
in these columns is small compared with the 
whirling velocity due to the rotor, and so makes 
a very small contribution to the vector sum of the 
two. In this connection it is noteworthy that 
the observed critical rotor speed corresponds to a 
value for Lb p » of approximately 3,000, not 
only for the column deseribed in [1], but for the 
larger columns of [3] and [4] also, which similarly 
have small b 7,, values (0-019 and 0-016 respec- 
tively)*. 

Column 2 fitted with the copper rotor, however, 
has a b/r,, value (0°27) intermediate between 


*The authors of [4] conclude that rotor surface speed 
alone is the criterion governing the abrupt change of 
column performance. This appears unlikely on the basis 
of the present discussion, and the results given do not 
establish this view. 
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LENGTH 


LATES/6Icm 


POINT OBTAINED WITH 
PPER ROTOR AND 
SWEEPING WIRES 


200 300 406 500 666 760 800 S00 


BOIL-UP RATE, mi/hr 


Fic. 5. B Smooth copper rotor 
D (>) Copper rotor with wire fins 
Westhaver theory for laminar vapour 
F Westhaver theory for turbulent vapour 
Rotor speed = 1470 rev min, 


the two values to which the curves in Fig. 4 
relate; this diagram should therefore furnish 
a turbulence criterion for the system. The fixed 
operating speed (1470 rev, min) for the copper 
rotor corresponds to a modified Taylor number 
of 1100 —-1200*; as shown in Fig. 4 this implies 
a very sensitive dependence of flow regime on 
vapour rate. Thus at low boil-up rates the 
vapour stream would be expected to be laminar, 
but at high boil-up rates the flow would be 
expected to become turbulent according to 
Fig. 4, although the rotor speed is here not high 
enough to satisfy the simple criterion for tur- 
bulence proposed in [1]. This effect would tend 
to flatten, or make more horizontal, the curve of 
separating power versus throughput for this 
column operating at 1470 rev min in comparison 
with that for the same column fitted with the 
glass rotor, for which the 6 r,, value is much 


*The value of } used is in all cases that for the annular 
vapour stream; ie. difference of column cylinder radii 
less reflux film thickness. Hence the Taylor number is a 
function of boil-up rate as well as rotor speed. 


smaller and the dependence of vapour regime on 
axial flow rate correspondingly less marked. 
For, as the boil-up rate is increased the separating 
power in any case tends to fall as the radial 
vapour-tiquid equilibrium is progressively less 
able to establish itself; but in the case of the 
copper rotor system it might be expected that 
this fall would be to some extent offset by the 
improved radial mixing at higher vapour rates 
due to the development of turbulence. This 
introduces the interesting possibility of designing 
rotary concentric cylinder columns giving good 
separation at high vapour rates (i.e. large through- 
puts for a compact unit) at conveniently low 
rotor speeds; much smaller 
dependence of separating power on throughput 
than is customary with film type columns, 
Figure 3 shows that the experimental separating 
power-throughput curve for the copper rotor is 
indeed flat compared with that for the glass rotor, 
and in Fig. 5 this curve is compared with the 
theoretical curve calculated from WestHaver’s 
equation [1] on the assumption of laminar 
vapour flow. At low vapour rates the performance 
of the column is evidently not very different from 
that to be expected on the latter assumption ; 
but at higher vapour rates the plate values are 
very much greater than those predicted, and the 


and showing a 


experimental curve agrees much more closely 
with the Westhaver curve calculated on the 
assumption of vapour turbulence*. Two points 
corresponding respectively to the lowest vapour 
rate used and to a vapour rate at which the 
experimental curve of Fig. 5 lies closer to the 


*The theoretical curves used here are based on moditied 
forms of WerstHaven’s equations, in which the liquid 
phase resistance is assumed negligible, in accordance with 


the experimental findings reviewed above. However, it 
must be pointed out that in view of the possible presence 
of vortices in the vapour stream, of which Wesrnaver’s 
treatment takes no account, their applicability here is 
questionable. A comparison of the curves A and A’ of 
Fig. 3 for the glass rotor system indicates the extent to 
which the Wesruaver laminar theory fails in a case for 
which the above criteria indicate a laminar + vortex 
regime. On this basis the discrepancy between the 
experimental and laminar theory curves (B and E of Fig. 5) 
for the copper rotor is apparently greater than would be 
accounted for by the presence of the vortices alone. 
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upper (turbulent) Westhaver curve, are plotted 
on Fig. 4. Though the curves of that figure do 
not extend to high enough Taylor numbers to 
establish the conclusion, it seems that these points 
might well lie above and below the boundary 
between the turbulent + vortex and laminar + 
vortex flow regimes. It thus appears that the 
overall performance of this column can be related 
to the general fluid mechanical observations of 
Kaye and Ercan. 


CONCLUSIONS 


The form of rotary concentric cylinder frac- 
tionating column described (column 2) is’ of 
practical use as a laboratory column, having no 
exposed bearings, lending itself to all-glass con- 
struction and having an exceptionally high 
performance. The principles of operation of 
such columns are thus of practical interest ; 
and as the laws governing their performance can 
be related very directly to hydrodynamic laws, 
their behaviour throws light on questions of 
general interest in fluid mechanics. 

The results obtained with the columns described 
are apparently substantially in accordance with 
the flow regime diagram of Kaye and E1car 
relating to the technically important case of fluid 
flow in the annulus between a rotor and stator ; 
but consideration of the results of other experi- 
ments with rotary concentric  fractionating 
columns suggests an important modification of 
this diagram. 

The apparatus described affords a direct means 
of investigating the relative importance of liquid 
and vapour phase mass transfer resistances in 
rectification. For the mixture of n-heptane and 
methyl cyclohexane studied, the liquid phase 
resistance at atmospheric pressure is shown to be 
negligible, and the column performance 


b = annular width 
= mean radius of annulus 
=~ abscissa 


Tm 


ordinate in Pig. 4 


Fg ~ geometrical factor in Taylor's equations 


| (; 1 0-652 [(b/r,,)/(1 — b2r,,)) | 
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NOTATION 


governed solely by limitations of radial mass 
transfer in the vapour. Thus, stirring the reflux 
film does not affect the column ; but the addition 
of fins to the rotor improves the separating 
power of the system, as does the development of 
general turbulence in the vapour in the experi- 
ments of Wintincuam et al. There is, however, 
little change in column performance at the rotor 
speeds at which vortices would be expected to 
appear. Possible causes of this anomaly are 
under investigation. 
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Key to Figs. 1, 2(a) and 2 (b). 


A Fixed central shaft 

BB Upper flange of motor housing 

'Thin-walled stainless steel motor housing 
1) Upper ball-race of squirrel cage rotor 

Squirrel cage rotor 

F Motor stator 

G Lower ball-race of squirrel cage rotor 


It Hollow shaft of squirrel cage motor 

1 Ball race retainer and locking ring 

Drive pin 

K Lower flange of motor housing 

1. Driving member of universal coupling 

M Inner driven member of universal coupling 

Cork-lined clamping ring 

0 Body of glass-metal O-ring union 

O-ring sealing plate 

Q = Spiral water condenser 

Glass column rotor 

S Take-off valve solenoid 

TT Cork lined clamping ring supporting outer column 
tube 

Centering screws 

Glass outer column tube 


W Compensating heaters 

X = Thermistors mounted on metal column jacket 
Y= Thermistors mounted on outer column tube 
Z = Lower bearing assembly 
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K mean circumferential vapour velocity /V, 
Re, = critical Reynolds number for turbulence for fluid flow between flat parallel planes, based 
on mean velocity of fluid 
Re, = critical Reynolds number for turbulence for circumferential shearing flow between concentric 
cylindrical surfaces, based on velocity V, of moving convex surface 
V mean linear axial velocity of fluid 
Vo surface speed of rotor 
Y = 2b py 
» = fluid viscosity 
p thuid density 
w = angular velocity of rotor 
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Abstract 


The droplet collection efficieney of wire gauzes exposed normal to an airstream 


has been measured and also computed using a rather erude theoretical model of the processes 


occurring. A spinning-top sprayer was used to produce a uniform cloud of dyed oil droplets 


in a wind tunnel; the cloud concentration was found by exposing a cascade impactor and measur- 


ing the deposit colorimetrically. 


The theory predicted the variations of collection efficiency 


with droplet size and closeness of wires, but afforded no explanation of the anomalous low experi- 


mental values found with 20 « diameter droplets 


Résumé 


Les toiles meétalliques en fer exposées normalement a un courant dair ont une 


capacité pour rassembler les particules qui a été non seulement mesurée, mais calculée a Vaide 


une theorie quelque peu grossiere des processus reels. Un pulveérisateur a téte tournante a servi 


a former un nuage uniforme de gouttelettes Qhuile colorées dans un tunnel. La concentration 
en gouttelettes a été trouvee par mesure colorimétrique du dépot sur un séparateur a 


impacts successifs. La théorie prévoit les variations de leflicacité de rassemblement en fonction 


de la dimension des particules et des dimensions des mailles, mais n'explique pas les valeurs 


experimentales anormalement basses trouveées avec des particules de de diameétre. 


Zusammenfassung 


Die Wirksamkeit des Tropfenfanges von Drahtsieben, die senkrecht zu 


einem Luftstrom stehen, wurde gemessen und unter Benutzung einer rohen Modellvorstellung 


berechnet. 


Kin retierender Zerstiuber erzeugte eine gleichfOrmige Wolke von angefiirbten 


Oltropfen in einem Windkanal; die Konzentration in der Wolke wurde mit Hilfe eines Kaskaden- 
Auffiingers kolorimetrisch bestimmt. Die Theorie konnte die Wirksamkeit des Tropfenfanges 
als Funktion der Tropfengrésse und Maschenweite voraussagen, nicht aber die anomal niedrigen 


Versuchsergebnisse im Bereich von Tropfendurchmesser. 


INTRODUCTION 

Wine gauzes are very useful for sampling small 
droplets drifting with the wind, as. unlike the 
cascade impactor [10], no elaborate subsidiary 
equipment such as and batteries are 
needed. They have been used to sample insee- 
ticide sprays [1, 2] ete. and natural clouds and 
fog [3, 4}. 

When a piece of gauze is placed in the atmos- 
phere some of the air originally flowing towards 
it flows around it and some through the holes, 
Particles in the atmosphere tend to follow’ the 
airflow, but because of their greater inertia they 
do not deflect as readily as the air. 

Therefore the air reaching the gauze contains 


pumps 


*Present Address: 164 Tremona Road, Southampton 


less particles than are found in the same cross 
section upstream, because some of the particles 
are deflected around the whole gauze network. 
The concentration of particles in the air-stream 
which eventually passes through the gauze is 
increased as the smaller number of particles are 
mixed with an even smaller proportion of the air. 
As this aerosol comes very close to the gauze, 
the air flows around the cylinders (through the 
holes) and again some particles take a straighter 
path and hit the cylinders. While particles of 
some kinds may bounce off, oil droplets deposit 
on metal surfaces. The relation between the 
deposit on a plane gauze normal to the airstream 
and the upstream aerosol concentration has been 
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measured experimentally in a wind tunnel, and 
also calculated using a simple, approximate 
theory. 

The results are discussed in terms of the collec- 
tion efficiency (FE) which is the deposit on the 
gauze as a percentage of the liquid which would 
have deposited if the droplets had not been 
deflected by the air currents. It is expressed 
mathematically by the equation 


100 = 100 Q 
A Heth (1) 


the volume of droplets deposited 

the cross-sectional area of the gauze 
in the plane normal to the airstream 
actual (projected) frontal area of 
material 

area of holes per unit area of gauze 
the concentration of droplets, in 
volume of droplets per unit volume of 
air 

the free stream velocity 

the time the object is exposed to the 
droplet cloud 


EXPERIMENTAL METHOD 


The collection efliciency was measured by 
comparing the deposit on wire gauzes with the 
deposit on the first stage slide of a cascade 
impactor when each were exposed in turn to a 
fairly uniform cloud of droplets in a horizontal 
wind tunnel. Using equation (1) for each object 
we have, 

(2) 
where the subscript ¢ refers to the cascade 
impactor. All the ratios on the R.HLS. were 
measured, The collection etliciency of the 
cascade impactor (£,) was obtained from some 
recent measurements by May [9] assuming the 
results obtained when sampling at 17} 1/min 
ina 3 miles hr wind would apply when sampling 
at 10 | min in a 1m see wind. 

A perspex wind tunnel with a 1 ft diameter 
working section was used. The spinning top 
sprayer [11] situated 3 ft upwind of the working 
section produced a cloud of uniformly sized 


droplets, but disturbed the airflow somewhat. 
To minimize this a coarse wire gauze was placed 
immediately downwind of the sprayer. The 
characteristics of the airflow in the working 
section are shown by the profiles in Table 1, 
which were obtained using a shielded hot wire 
anemometer [12]. 

Table 1. Wind speed (em/sec) — viewed from 

upwind 


Left Zin. Lin. Centre lin.  2in. Right 


2in. up 105 115 105 | 102 
Centre 105 | 107 107 
Zin, down 1 105 102 


The variation of the droplet cloud concentration 
in space, was found by exposing the cascade 
impactor at various positions across the working 
section, and measuring the dyed spray deposit 
colorimetrically using a Unicam Absorptiometer. 
The results of one of these tests using droplets 
of 35 u dia. is shown in Table 2. The deposit 
shadow is associated with the wind shadow,” 
partially because less air carried less droplets. 
The concentration also varied in time, ranging 
» 11 per cent about the mean value. 


Table 2. Concentration variations (arbitrary units) 


lin. Left Centre — Lin. Right 


} in. above centre 1s 14 20 
} in. below centre 20 17 19 


The liquids, kerosene and Kerosene Extract 
Bottoms [K.E.B.] dyed with 1 per cent and 
8 per cent grasol blue respectively, were fed to 
the spinning top at 1-2cm* min. Three wire 
samples (2 x 2cem, 3 x 3em and 4 x 4 em) 
were exposed to the spray alternately with the 
cascade impactor for intervals of 1-9 min. The 
deposits were washed off with petroleum ether 
and the quantity estimated colorimetrically 
using the Unicam Absorptiometer. This was only 
roughly linear, so the exposures were arranged to 


where Q - 
IL. 
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produce similar pointer deflections. The size of 
the droplets being used was found by exposing 
an MgO coated slide [10, 13] in the cascade 
impactor and measuring the diameters of the pits 
produced under a low powered microscope. 


ResuLtTs or THE EXPERIMENTS 


The wind speed was restricted to lms and 
the variation of the collection efliciency with 
droplet size and the closeness of the wire mesh 
was studied. The results are shown as plotted 
points on Figs. 1, 2 and 3 and are tabled in the 
appendix. We see that the collection efliciency 


increased with droplet size and decreased with 
the closeness of the wire mesh. 


Curves —caiculated 
Points -measured 


8 


64 66 08 to 
Blockage coefficient, 3 


— 


The variation of collection efficiency with 8. 
(8 » dia. K.E.B. droplets). 


Curves -calculcted 
Points - measured 


Collection efficiency, 


Blockage coefficient, 


Fig. 2. The variation of collection efficiency with 8. 
(30 dia. kerosene droplets). 


Curves -calculated 
Points-measured 


Collection efficiency, 


20} 


1°) 10 20 30 40 50 60 
Droplet diameter, j 
Fig. 3. The variation of collection eflicieney with droplet 


THEORETICAL Stupy 

To summarize these results and to make 
interpolations from them easier, a brief study of 
the theory of deposition of windborne droplets 
onto objects was made. Analytical methods of 
calculating the percentage of droplets depositing, 
using first principles, did not seem promising ; 
some of the dilliculties have been described by 
Ranz and Wonc [7]. However numerical methods 
have been used successfully to calculate the 
trajectories of droplets in an airstream flowing 
around isolated, infinite cylinders by Davies 
and Prerz [8] and others, and the fraction of the 
airflow that goes through a gauze has been 
established empirically [6]. Therefore it) was 
possible to postulate possible models of the 
processes occurring, and to calculate how the 
collection efficiency varied with droplet size, 
wire size and spacing, and wind speed in these 
theoretical cases. In each case it would be 
necessary to make a number of assumptions that 
could not be justified a priori. In the model 
described here the simplest and most obvious 
assumptions were made and the collection etlicien- 
cies calculated from it simulate the variation with 
closeness of wire, and with a notable exception 
the variation with droplet size. This simple 
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and unsophisticated analysis of the problem 
may throw some light on the physical processes 
occurring. It also shows the accuracy obtained 
when such calculations are used to assess the 
deposition occuring on complex objects; this 
is often necessary in practice. 

It has been found that the fraction of the air 
that flows through the gauze (u/U) is given by 
the equation 


u/U (1 + (3) 


where k is the resistance coeflicient of the gauze 
[5]. Simmons and Cowprey [6] measured k 
and found that with a square mesh of cylinders, 
and « U < 10m sec to a first approximation, 


(4) 


In the computation, the complex airflow and 
particle trajectories near the square mesh gauzes 
were assumed to approximate to an aerosol of 
concentration (c) flowing at the reduced velocity 
(u) towards isolated infinite cylinders of the same 


diameter as the wire. In this case, 


where E’ is the collection efliciency of an isolated, 
infinite cylinder in a wind of speed (1). 

For comparison with the experimental results 
E’ was computed using the results of Davies and 
Prrrz [8], who expressed the collection efliciency 
of an isolated, infinite cylinder in terms of the 
dimensionless parameter P, the Reynolds number 
(Re) and the ratio of the particle radius (} Dp) to 
the cylinder radius (4 De). P for spherical particles, 
is defined by the equation 


_p 
9» 


(6) 
where p is particle density and » the coeflicient 
of viscosity of air. The Reynolds number of the 
flow around an infinite cylinder is given by the 
equation 


Re = uD, Pa 


where p, is the density of the air, and wu its 
velocity. Sedimentation of droplets onto hori- 
zontal cylinder was not considered by Davies 


(7) 


and Prrrz [8] but using the methods of Ranz 
and Wone [7] it was shown that negligible 
sedimentation onto the horizontal cylinders of 
the mesh occurred. Values of P and Re in the 
various experiments were computed, taking 
the coellicient of viscosity of air to be 1-79 « 10-3 
poise, and the density of air, kerosene droplets 
and Kerosene Extract Bottoms droplets to be 
10-5, 0-82 and 0-92 g/cm’ respectively. It 
was found that Re varied from 6 to 18, so the 
results of Davies and Prrrz [8] for Re — 10 
were used in the computations. The collection 
efliciencies computed in this way are compared 
with the experimental results in Fig. 1, 2 and 3, 
and are tabulated in the appendix. 


Discussion oF RESULTS 


Figures 1, 2 and 3 show that the simple theory 
tends to underestimate the collection elliciency, 
and the mean square deviation between calculated 
and measured collection efliciencies is 13 per cent. 
The underestimate probably arises because the 
droplet inertia and concentration are under- 
estimated ; the droplets decelerate in the slower 
airstream near the gauze, but probably not 
to as low a speed as the airflow (u) and so more 
droplets deposit than in the quasi-equilibrium 
cease used in the calculations. The variations 
with closeness of wire and droplet size are well 
simulated by the theory, apart from the un- 
expected anomalous experimental results with 
approximately 30, diameter droplets which 
may be due merely to an extreme undetected 
experimental error; a possible, but unlikely, 
cause is that for some unexplained reason the 
30 « diameter droplets adhere less firmly to the 
wires than other oil droplets. 

When applying these results to droplets of 
other liquids it should be remembered that they 
may not adhere in the same way as oils [14]. 
Also in winds exceeding 10 m/sec the airflow 
through the gauze is different [6], and this fast 
turbulent airflow might increase the collection 
elliciency [3] so strictly speaking, these results do 
not apply to those conditions ; though they may 
be the best guide available to the field worker ! 
In the non-turbulent conditions (wind < 10 m_ sec) 
the collection efliciency is independent of the size 
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of the sheet of gauze exposed normal to the 
wind, and so these results can be applied to large 
sheets of wire gauze. In the open, atmospheric 
turbulence might effect the collection etliciency, 
but no information is to hand on that, and these 
results have been useful for interpreting aerial 
spraying trials of insecticides and other liquids 
(15, 16). 

The correspondence between the results 
obtained by computation from a very unsophis- 
ticated theoretical model and the measurements 
shows that simple calculations using results 
obtained for the flow around simple smooth 
objects, may be quite useful, though not highly 
accurate, when used to predict the deposition of 
droplets onto various complex objects in- the 
field. 


Acknowledgements—-The writer would like to thank 
Mr. P. T. Kine for assistance with the measurements, 
and Drs. D. J. Tuomas and R. C. Pankuurst for calling 
his attention to References 5 and 6, which were so useful 
in the theoretical study. He is grateful to Begg, Cousland 
& Co, Ltd, and N. Greening & Sons Ltd. for kindly 
providing the variety of wire gauzes used. 


NOTATION 


A = cross-sectional area of gauze 

= intake area of cascade impactor 
frontal area of actual gauze wires 

- concentration of droplet cloud 
wire diameter 
droplet diameter 
collection efliciency of gauze (defined in equation 
1) 
collection efliciency of cascade impactor 
collection efliciency of isolated infinite cylinder 
resistance coefficient of gauze (defined in 
equation 3) 

= inertia parameter (defined in equation 6) 

= volume of droplets deposited on gauze 

= volume of droplets collected by cascade impactor 


Re = the Reynolds Number of airflow around an 
isolated cylinder (equation 7) 
t — period of exposure of gauze 
t. = period of exposure of gauze cascade impactor 
U = free stream velocity 
u = average wind speed through gauze 
B = area of holes per unit area of gauze 
p = particle density 
Pq = air density 
» ~ coeflicient of viscosity of air 


APPENDIX 
Table of the measured and calculated collection efficiences 


Collection efficiency 
Wire No. of Liquid Droplet (%5) 
size) wires dia, 
used 


(SUWG ) in. (yu) Measured Calculated 


37 Is | kerosene . | 
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Abstract An interferometric technique is described by which the absorption of carbon dioxide 
into quiescent water has been examined, There is no measurable resistance at the interface to 
solution of the gas into pure water, the rate of absorption into the liquid being governed by 
molecular diffusion. A measurable interfacial resistance is found when carbon dioxide is absorbed 


by solutions of surface-active agents. 


Résume Les auteurs decrivent une methode miterferometrique qui permet detudier Padsorption 
du CO, dans Peau au repos, La résistance a Tinterface de la solution du gaz et de Peau pure est 
trop faible pour étre mesurable, la vitesse dadsorption dans le liquide étant déterminée par la 
diffusion moléculaire, Les auteurs trouvent une résistance interfaciale mesurable quand le CO, 
est absorbé par des solutions dagents tensio-actifs. 


Zur Bestimmung der Absorption von Kohlendioxyd durch ruhendes 
Da bei der Liisung des Gases 


— 
Wasser wird eine interferometrische Versuchstechnik beschrieben, 
in reinem Wasser kein messharer Grenzflichenwiderstand vorhanden ist, wird die Geschwindigkeit 


suned 


der Absorption in die Flissigkeit durch die molekulare Diffusion bestimmt. 
\bsorption ven Kohlendioxvd durch Lésungen von 


Grenztliichenwiderstand wird bei der 
obertlichenaktiven Stoffen gefunden. 


INTRODUCTION 


A vaniery of techniques have been devized for 
the purpose of detecting and measuring interfacial 
resistance to the absorption of a gas by a liquid 
2. 3. 5, 6. 8, 11, 13). In all cases the experi- 
mental arrangements have been such that liquid 
and gas have been moved relative to each other. 
At the time when the present work was 
commenced, in no case was the exact nature of 
the motion of the liquid known and its effect 
upon conditions at the interface could not, 
therefore, be reliably assessed. Also the rate of 
absorption was determined from the total 
quantity of gas dissolved over a relatively long 
exposure time and the size of any interfacial 
resistance could not be caleulated from = such 
data with accuracy. For these reasons it seemed 
desirable to obtain some data on the absorption 
of a pure gas by a stationary liquid ; for by using 
a pure gas it is ensured that there is no resistance 
to transfer in the gas phase, and by using a 


messbarer 


stationary liquid the unknown influence of the 
motion of this phase is thereby eliminated, In 
the work to be described the pressure of carbon 
dioxide over the surface of pure water and over 
solutions of surface-active agents was suddenly 
altered. The non-steady absorption of the gas 
into the liquid, resulting from this change of 
conditions, was followed by an interferometric 
method. It will be shown that even a small 
interfacial resistance (0-25 see/em) can be 
measured with accuracy from data obtained in 
this way. 

Several of the earlier workers [3, 5, 8] found 
an apparent interfacial resistance to the absorp- 
tion of gas by water. More recent work [1, 2, 6, 
11, 13], however, has shown that the resistance 
must be very small. In particular Lynn et al. 
[11] and also CULLEN and Davipson [1] measured 
the quantity of gas absorbed by a liquid flowing 
over the surface of a sphere: by assuming that 
the liquid was inviscid, the form of the motion 
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was deduced. Diffusion coeflicients of carbon 
dioxide into water, calculated from data obtained 
by Cutten and Davipson on the assumption that 
there is no interfacial resistance, have the same 
value as diffusion coeflicients measured by other 
means. When absorbing carbon dioxide into 
solutions of surface-active agents CULLEN and 
Davipson did measure a resistance at the inter- 
face at low concentrations, but at high concentra- 
tions of surface-active agents they found no 
resistance. Measurements on similar systems 
have been included in this work. 


Basis or THE INTERFEROMETRIC 
Metruop 


The method depends upon the measurement of 
the shift of interference fringes near the surface 
of a liquid, caused by a change in gas concentra- 
tion in the liquid, For this purpose a Mach- 
Zehnder interferometer has been used. This 
instrument has four optical plates (Fig. 1), of 
which two (M2 and M3) reflect all the incident 
light. The other two plates (M1 and M4) are 
semi-reflecting ; that is, about one half of the 
incident light is reflected and the rest is trans- 
mitted. When the instrument has been adjusted 
correctly the surfaces of all four plates are 
parallel and their centre points lie at the corners 
of a parallelogram. Path M1, M3, M4 should 
equal path M1, M2, M4 to within a few wave- 
lengths of the light to be used. 

Light from a high pressure mercury vapour 
lamp is directed by a lens of short focal length 
onto a slit, which acts effectively as the light 


source. The wedge of light formed at the slit 
is collimated by a second lens to yield a parallel 
beam of light which is passed into the interfero- 
meter, falling on plate M1. The light reflected 
at this plate is reflected again at plate M2 and 
transmitted through plate M4. The light trans- 
mitted at plate MI is reflected at plate M8 and 
again at plate M4, where it recombines with the 
light from the first path. A pair of compensating 
plates are introduced into one of the light paths. 
These provide an optical balance to the cell in 
which gas absorption is carried out and which is 
placed in the other path. On the emergent side 
of the interferometer a lens of long focal length 
forms an image of the absorption cell. This 
image is magnified by a microscope objective and 
thrown onto a screen or photographed. The 
plates are positioned so that the emergent beams 
intersect, virtually, in the plane of the absorption 
cell. Interference fringes, from which all measure- 
ments are taken, are formed in this plane, and a 
magnified image of the absorption cell and the 
fringes are obtained together. The constructional 
details of the instrument used in this work and 
the method by which it was adjusted have been 
described in full elsewhere [7]. 

The position of the interference fringes depends 
directly upon the difference between the optical 
path lengths of the two light beams, and the 
optical path length is directly proportional to the 
refractive index of the medium containing the 
light. If a gas is absorbed by a liquid placed in 
one of the two light beams, a non-uniform 
distribution of gas in the liquid will cause a 


\\Wosted 
2 M4 \ 


Screen 


Compensating 
plotes 


M3 


. Optical system and interferometer. 
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E, A, 
non-uniform distribution of the refractive index in 
the solution ; and this will appear as a distortion 
of the fringes. 

Consider fringe, initially 
perpendicular to the liquid surface, indicating 
that the liquid has a uniform concentration and 
temperature. At some time after the commence- 
ment of absorption of gas through the interface 
this fringe will be distorted, as shown in Fig. 2. 
The optical path length through the liquid is the 
physical path length multiplied by the refractive 
index, and the physical path length remains 
With the gases used in this work, at 


one straight) and 


constant, 


Surface 


Fig. 


the concentrations encountered, the refractive 
index of the solutions is linearly dependent upon 
concentration (mass volume), and therefore for 
an isothermal absorption process 


Gas entry 


Hanvey and W. 


(1) 


where ¢ is concentration and 
NX is fringe shift in the liquid (Fig. 2). 
Subscript s refers to the surface of the liquid and 
0 to the bulk of the liquid (where the 
concentration and temperature re- 
main unchanged from their initial 
value). 


If the solution of gas at the surface of the 
liquid is in equilibrium with the gas phase, 
the surface has the saturated concentration c*, 
which is known in terms of the gas pressure over 
the liquid. The fringe shift Y*, given by a 
concentration change ¢* — cs, can be determined 
from a photograph of fringes taken after a long 
period of contact, when the surface will be 
saturated whatever the interfacial resistance. 
From the measured value of V* and the known 
solubility the relationship between and 
can be calculated, 


or ture ABsorrTION 


The cell in which gas was absorbed into liquid 
was a horizontal cylinder of stainless steel (Fig. 3). 
Optical flats formed windows which were placed 
in recesses at each end of the cell. The windows 
were sealed against an “O°” ring of synthetic 


Fic. 3. Section t 


hrough the absorption cell. 
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rubber, Screwed over each end of the cell was a 
locking device, carrying three screws bearing on 
the windows opposite the “O” ring. By adjust- 
ing these screws the cell windows were aligned 
parallel, the distance between the windows being 
measured by a comparator capable of comparing 
distances to within 1/25,000in. The cell was 
supported in the light path by means of a cradle 
which could be adjusted laterally and vertically 
and could also be rotated through a small angle. 
By altering the position of the cell until an 
image of a small light source reflected by the 
cell windows coincided with that source, it was 
ensured that the light path was parallel to the 
liquid and perpendicular to the cell windows. 

A boss diverted the gas, on admission, parallel 
to the axis of the cell. A stainless steel gauze 
(200 mesh in) was placed between the boss and 
the liquid surface, with the object of smoothing 
the movement of the gas towards the liquid. 
The effect on the liquid of a sudden admission of 
gas into the cell, from a pressure of 5 em Hg. Abs. 
to 1 atm, was investigated by observing an image 
of the surface magnified one hundred times and 
projected onto a screen. On admission of the gas 
lycopodium and aluminium powder floating on 
the surface of the liquid suffered no_ visible 
displacement, no swirling of the liquid was 
observed and the liquid surface remained at the 
same level. Before the admission of the gas the 
interference fringes were perpendicular to the 
liquid surface and stationary, showing that 
there were no temperature differences or convec- 
tion currents within the liquid. After admission 
of the gas the fringes in the bulk of the liquid 
remained perpendicular to the liquid surface, 
showing that no convection currents between 
the surface and the bulk of the liquid had been 
set up. It was concluded that gas entered the cell 
without moving the liquid and that distortion 
of the interference fringes was caused solely 
by the absorption process. Experiments in which 
gas was released from the absorption cell also 
showed no hydrodynamic disturbance. 


PHOTOGRAPHY AND MAGNIFICATION 


The interference fringes were photographed 
with a 16mm cine-camera which operated at 
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32 frames/sec. A microscope objective of 16 mm 
focal length, fitted to the camera, gave a magni- 
fication of about ten diameters at the photo- 
graphic emulsion, The width of emulsion exposed 
was lem and therefore a fringe spacing of 
3/20 mm enabled six fringes to be viewed. The 
magnification given by the optical system on the 
emergent side of the interferometer, before the 
camera, was measured by photographing a 
slit of known size in the plane of the image of the 
absorption cell and the fringes. The magnification 
between this plane and the absorption cell was 
determined by measuring the size of the image 
of the end of the absorption cell. In the course 
of an experiment the camera was operated 
continuously for 1 see before admission or release 
of gas, in order to ensure that the camera had 
accelerated to working speed, and then for 4 sec 
afterwards. Thereafter sets of two or three 
frames were taken at 7, 10, 15 and 20 see after 
commencement of the experiment. 


CALIBRATION 


Apart from distortion of the fringes near the 
surface of the liquid, caused by a change in the 
concentration of dissolved gas, a sudden altera- 
tion of the pressure of gas in the cell moves all 
the fringes in the bulk of the liquid rapidly sideways. 
This is due to a change in the refractive index of 
water with pressure. The fringes, however, remain 
parallel to each other and in the same relative 
position, and the movement of the fringes is 
completed in less than the interval between two 
frames of the camera (0-08 sec). 

At the beginning of absorption, just after gas 
has been admitted to the cell, a movement of 
the fringes commences near the surface of the 
liquid, rises to a maximum value in about } sec 
and then dies away within 5 sec, leaving a fringe 
shift evidently caused by change in the concentra- 
tion of gas in the liquid. This preliminary 
movement of the fringes is due to the conduction 
of heat through the liquid. An alteration of the 
total pressure in the cell alters the partial pressure 
of the water vapour above the liquid and brings 
about condensation or evaporation at the liquid 
surface with the liberation or absorption of latent 
heat. 
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Heat conduction is a more rapid process than 
diffusion in a liquid, and the temperature in the 
liquid becomes uniform much more rapidly than 
does the concentration. However, for the same 
initial pressure change, the same amount of heat 
will be transferred whatever gas is being used, 
because the amount of water vapour condensing 
or evaporating at the surface will be the same. 
The resultant temperature distribution will also 
be the same, since the thermal conductivity of 
water is changed very little by low concentrations 
of gas. Air is only slightly soluble in water and 
has no observable effect on the refractive index 
of water when left in contact for periods which 
are long compared with the duration of the 
experiments. Thus, by subtracting the fringe 
shift caused by heat effects after admission of air 
from that observed during the absorption of carbon 
dioxide, each measured at the same time after 
the commencement of an experiment, the fringe 
shift due solely to the change in the concentration 
of carbon dioxide in the water is obtained. 


INTERFACIAL Resistance 

ABsorrPTION OF CARBON DIOXIDE 

An interfacial resistance (1/k,) is usually 
defined by the relationsihp 


(9Q/20), (e* — ¢,) (2) 


where k, is the interfacial mass transfer coeffi- 
cient, 
Q is the mass of gas absorbed through 
unit area of interface in a time @ and 
x is distance normal to the interface. 


For uni-directional diffusion in the non-steady 
state 
= D (3) 


where D is the diffusion coefficient. 

The solution of equation (3) for diffusion into a 
“ semi-infinite ” liquid with initially uniform 
concentration ¢y, using equation (2) as a boundary 
condition, is 


Values of k, are most easily computed from a 
knowledge of the change with time of the con- 
centration of the gas at the surface of the liquid ; 
and this is determined from the variation of 
fringe shift with time at the surface of the liquid 
(X,) recorded on the interferograms. A_ theore- 
tical relationship for the surface concentration in 
terms of time is found by putting 2 — 0 in 
equation (4). 


—_ 


This cannot be converted to an explicit equation 
for k,, but provided that k,\/@/D is greater 
than about 3, equation (5) can be written 


- 
= 1— (6) 
In terms of fringe shifts 


NX, N* — X* (1k) y/(D/) (1/0) (7) 


A plot of X, versus 1/4/@ should therefore be a 
straight line, for sufficiently small values of 1/1/0; 
this line will cross the X, axis at Y*. From the 
slope of this line and equation (7), k, can be 
approximately determined. Using this value of 
k,, X, can be plotted against 1/\/0@ according 
to equation (5), and it can be checked whether 
or not the measured values fall about this curve for 
high values of @. When there is no interfacial 
resistance, that is, k, is infinite, X, will be constant 
at the value X*. 

For the absorption of carbon dioxide into water 
X, was found to have a substantially constant 
value, corresponding to the saturated concentra- 
tion, X* (Fig. 4 and Table 1). Within the time 
interval of one frame of the cine-camera (about 
1/50 sec) the interference lines assumed that form 
which would be caused by molecular diffusion 
of the gas through the liquid with the surface 
of the liquid saturated. This means that k, must 
be at least about 4 cm/sec, the highest value which 
could be detected with certainty. (There must 
always be some “ resistance ” at the interface, for 
the highest rate at which gas can arrive there is 
fixed by the value of the kinetic collision number. 
This limiting value of k, for carbon dioxide at 
1 atm. and 0°C is 5240 cm/sec). 
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Fro. 4. Carbon dioxide absorbing into distilled water. 
Initial pressure 75-1 em Hg Abs. 
Final pressure (i) 34-1 em Hg Abs. 
(ii) Hg Abs. 


3 


4 1 
sec 


Fic. 5. Carbon dioxide absorbing into 0-41% lissapol-N, 
k, = 0-028 cm /sec. 


Table 1. Transfer of carbon dioxide into distilled 
water. 


The following experiments were carried out. In 
no case was a resistance detected. 


| Final 
pressure 
(em Hg Abs) 


Initial 
pressure 
(em Hg Abs) 


change 
(em Hg) 


1-7 
19 
21 
751 
46-9 
47-3 
28-8 
200 
29-0 | 


32-3 
741 
55-4 
341 

44 
74-4 
28-2 
10-4 
743 


For carbon dioxide dissolving into aqueous 
solutions of surface-active agents X, was not 
independent of 1/,/@ and the analysis described 
above was used to calculate k, (Fig. 5 and Table 2). 
The interfacial resistance was measured at a low 
and at a high concentration of two surface-active 
agents. It was found that the resistance was 
substantially the same in all four cases. 

In order to be able to compare this work with 
the work of CULLEN and Davinson [1], about the 
same concentrations of the same surface-active 
agents were used. At the low concentration 
CuLLEN and Davipson found an_ interfacial 
resistance of the same magnitude as that found 
here, but at the high concentration they found no 
resistance. In this work the resistance was the 


Table 2. Transfer of carbon diowide into solutions of surface-active agents 


Cone. (% by wt.) 


Initial pressure Final pressure ke, (cm /ece) 


Lissapol 


agent 


Teepol 


141 82-1 


13-7 
14-0 | 


30-4 
30°3 
30-3 


(em Hg Abs) (em Hg Abs) 
| 


140 


— TO-7 10-8 
10 “6 4 5 : 10 — 36-9 100 
59 sec” 45-5 | 76 
29 | 746 
| 
| 
| Eqn. Eqn 
— 
| 
” | 56 | 0-083 | 0-023 
8-6 0-286 0-028 
- 6-2 | 1-00 | 0-033 
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same at both concentrations. The difference 
between the two sets of results may be due to the 
fact that the liquid was motionless in our experi- 
ments, whereas in the experiments of CULLEN 
and Davinson the liquid was moving over the 
surface of a sphere. The same kind of result 
has been noted for solvent extraction. Thus, for 
example, LinpLanp and Tersesen [10] and 
Me.uus and Tersesen [12] extracted iodine from 
water by carbon tetrachloride in the presence of 
oleyl-p-anisidinesulphonate in a drop extractor 
and measured an interfacial resistance which 
rose rapidly with concentration of the surface- 
active agent to a constant value. Hoim and 


Tersesen [9] who carried out the same experi- 
ments in a stirred vessel, where the phases are in a 
disturbed condition, found that the rate of transfer 
of iodine passed through a maximum value with 
concentration of surface-active agent." 
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NOTATION 
concentration of gas in liquid (mass / volume) 
distance moved by fringe from its initial position 


mass of gas absorbed in unit time area of inter- 
face 


time after admission of gas 
distance perpendicular to the interface 
interfacial mass transfer coeflicient (distance time) 
diffusion coetlicient 

Subscript s surface of liquid, subscript 


condition of liquid, superseript * 
corresponding to gas pressure, 


initial 
equilibrium value 
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Abstract—A mathematical justification is given for the use of the “ Danckwerts * boundary 
conditions for continuous flow reactors, It is shown that the apparent indeterminacy, which 


DANCKWER?Ts resolves intuitively, is caused by the use of a discontinuous coeflicient of diffusion. 
By treating this as the limit of a continuous function and imposing continuity of the reactant 
concentration as the physically relevant boundary condition, the Danckwerts solution is obtained 


in the limit. 


Résumé— L’auteur donne une justification mathématiques des conditions aux limites utilisées 
par Danckwerrs dans le cas d'un réacteur 4 écoulement continu, Tl montre que Vindétermination 


apparente résolue intuitivement par DaNcKkwerts est inhérente & un coeflicient de diffusion 


discontinu. 


I] considére ce coeflicient comme la limite d'une fonction continue et il impose une 


continuitée a la concentration du réactant comme étant la condition limite physiquement correcte : 


la solution de DANCKWeERTs est alors obtenue a la limite. 


Zusammenfassung -Die Grenzbedingungen nach Danckwerrs fiir den kontinuierlich durch- 


stromten Reaktor werden mathematisch gerechtfertigt. Die anscheinende Unbestimmtheit, die 
Danckwerrs intuitiv auflést, ist durch die Verwendung eines diskontinuierlichen Diffusions- 


koeflizienten verursacht. 


Behandelt man diesen als Grenzfall einer kontinuierlichen Funktion 


und setzt die Stetigkeit der Konzentration des Reaktanten als die physikalisch entscheidende 


1. INTRODUCTION 


In an oft-quoted paper (Chem. Engng. Sci. 1953 2 
1), Danckwerrts has considered the steady state 
flow of reactant through a packed tubular reaction 
vessel in terms of a second order ordinary differen- 
tial equation (equation 30 in loc. cit.). This 
equation for the steady state concentration, c, 
of reactant in a first order reaction supposes that 
¢ is a function of one space variable only, y, the 
distance down the tube. The streaming velocity, 
u, the rate constant, k, and the (constant) coefli- 
cient of diffusion, D, enter as parameters into 
the diffusion equation, which may thus be written 

D- de _de 


(1) 


The boundary condition at the entry to the tube, 


Grenzbedingung fest, so erhilt man die Danckwerts-LésuNG im Grenzfall. 


+ Danckwenrts’ notation is retained for ease of comparison. 
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y = 0, where the diffusion coefficient discon- 
tinuously changes from zero, is obtained by a 
consideration of mass balance, and is 

de 


uc*=-uc—D-—, 


(2) 


where c* is the concentration in the entering 
stream. 

A similar relation is obtained for ¢ at the exit 
from the tube, y = L, where D again changes 
discontinuously, but is replaced by DancKWERTS, 
on intuitive grounds, by the stronger condition 


(8) 


Conditions (2) and (3) lead to a unique solution 
for ¢; however this solution, as presented in 


0 
de 
— =0, = 
| | 
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DaNcKWerts’ paper, appears to rely for its 
uniqueness on the acceptance of an intuitive 
boundary condition. If only because this boun- 
dary condition has not been universally accepted, 
it seems desirable to investigate a little more 
closely the formal mathematical implications of 
the idealization represented by equation (1) and 
the boundary conditions (2) and (3). 

First of all, we observe that in the general 
solution to (1), using boundary condition (2) at 
y = 0 and a similar condition at y = L, a discon- 
tinuity in c, at either y = 0 or y = L, or both, is 
necessarily consequent upon the imposed discon- 
tinuities in D. Within the reactor, i.e. where 
0 < y < L, we naturally consider only continuous 
solutions for ¢. This we do for sound physical 
reasons, but not for mathematical ones, for it 
may easily be verified that solutions exist for ¢ 
that are discontinuous within the region 0 < y < L. 
This suggests, as a next step, that we should 
eliminate, if possible, the discontinuity in D, and 
seck solutions for ¢ that are continuous everywhere, 
and are therefore physically acceptable. We may 
do this in a variety of ways. In the following 
section, we shall deliberately choose a form for 
D, the diffusion coefficient, which is continuous 
everywhere and which, by means of an obvious 
limiting process, may be made to approach as 
closely as we please to the “ top hat " distribution 
of Danckwexts. It will be shown that a unique 
solution for ¢ is in fact obtained without resort 
to intuition, and that this solution tends, in the 
limit, to the solution obtained by Danckwerts. 
In other words the apparently intuitive boundary 
condition (3) becomes the inevitable result of the 
limiting process applied to the continuous 
solution. 


2. MatnematicaL Treatment 
We replace the constant coeflicient of diffusion 
used by Danckwerts by the continuous form 
(i) D= when 0 < y < 1,4 (4.1) 
(ii) D=D, when 1 A cy <L— 
— (1 A) (4.1) 


(iii) D = AD,(L — y)when L — (1, A) 
<y<L (4.iii) 


where D, and A> 1 are constants. This is 
shown diagrammatically in Fig. 1. The diffusion 


Do 


Va 


Fic. 1. Coefficient of diffusion, D, as a function of y. 


equation which, for a non-constant coeflicient of 
diffusion, must be written 


D- 
dy dy 

becomes, in the three regions (i), (ii) and (iii), 
ty u k 

dy AD, 
dy Dy dy De 
AD,| dz 
where we have written z= L — y, (7) 


and ¢,, ¢, and c, refer to the solution for ¢ within 
the three regions. The boundary conditions that 
we must apply to this set of three equations are 


(i) = 0, (6.1) 


(ii) c = 0, (6.ii) 


(iii) {1+ 


= 0, (6.iii) 


at y = 0, 


aty = 1/4, | 
(8) 


aty=L 
— (1/4), (2 = 1/4), | 
since mass balance at y = 1/A andaty = L —1/A 


requires that de/dy as well as ¢ be continuous 


there. We have specified two boundary condi- 
tions for cy and ¢,, but only one for ¢,. Thus, at 
first sight, the system is apparently indeter- 
minate. 
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We now consider the general solution for each 
of the three equations (6). That for (6.ii) can be 
written directly as 


+- (gg EXp Ea 


and being arbitrary constants. This 
involves two linearly independent solutions which 
are finite in the region (ii). Those for (6i) and 
(6iii) are best considered in terms of ascending 
power series expansions about the points y = 0 
and z = 0. (It will be sufficient to do this because 
we will later apply the limiting process A — ©). 
Thus, if we write 

(ot 
and substitute into equation [6(i)] we find that 


r or AD,’ 
giving the two linearly independent solutions 
(k AD, — u)y + O(y*)...... 
and ¢,°) = by 4?) [1 + (k/ AD, + u)y + 
+ O(y*)....] 
that are bounded within the region (i). 
At y = 0, ¢," = ag, ¢, = 0. 
At y = 1/A, ¢,") = ag + 0(1/A?), 
= by (1/ A) 4?" + 0(1/A), 


e, = + 


(12) 


0 
(2) 
= + 0(1/ A), 
0 


since 1/A is small. 
Similarly if we write 

= 2 (dg +d,z+.....- d,z" + ...), (18) 
and substitute into equation (6iii) we find that 


s=00 


AD,’ 


giving only one solution 


A note on the “ Danckwerts "’ boundary conditions for continuous flow reactors 


= dg (1 + (k/ AD, + u)z + O(2*)... 
that is bounded in the region (iii). 


-) (14) 


At =z @, = dy. 


At z = 1/A, c,'") = dy + 0(1/A?), 
= ADy + u) + 0(1/ 


We now see that the five boundary conditions 
(8) are in fact suflicient to define a unique solu- 
tion for c, in terms of c*, since there are only 
five constants, [from (9)], [from (11)], 
by [from (12)] and d, [from (14)] to be determined. 

We need not evaluate explicitly these constants 
in terms of c*, but proceed at once to the limit 
A -»@. This corresponds, in the limit, to the 
discontinuous form for D chosen by Danckwerts. 
As 1/A -> 0, we see that 


(1/A) (0), 


but that ¢,"*) (1/4) +b, while (0) = 0. In 
other words the solution c," represents, in the 
limit, a discontinuity in ¢. Thus (2) is the correct 
limiting boundary condition. Similarly, as 
1/A +0, 
(1, A) (0) 

while (1/A) +0. Thus the strong boundary 
condition (8), implying continuity of ¢ at y = L, 
is recovered. 


8. Discussion 


It might be argued that the choice (4i), (4ii) and 
(4iii) for D is in some sense particular, and that 
a different limiting solution would be obtained by 
using a different continuous distribution for D 
within the regions (i) and (iii), This is not so, and 
the same limiting solution can in fact be demon- 
strated for a very wide class of functions D. 
What we have done in section 2 is to give but an 
illustrative example of the effect of the sign of 
u on solutions for ¢ in equation (5) in the neigh- 
bourhood of points where D becomes equal to 
zero. This purely mathematical dependence on the 


sign of u can be interpreted physically, as 


J. R.A. Prarson 


Danckwerts has done, in terms of different 
boundary conditions to be applied at the entry 
to or at the exit from a packed tubular reactor, 
where the diffusion coeflicient changes rapidly 
from zero. 

The same procedure can be adopted to obtain 
relevant boundary conditions in more complicated 
situations involving diffusion where there are 
several reactants, and where temperature varies. 


It can readily be shown that the boundary con- 


dition (3), which implies continuity of the 
variables at the exit, is recovered in every case. 
It need hardly be pointed out that the use of a 
discontinuous diffusion coeflicient so simplifies 
the solution of the diffusion equation that it 
enables analytical results to be obtained in many 
cases where a more physically acceptable form 
for D would make exact solution impossible. To 
this extent, the justification given above becomes 
of practical value. 
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Errata 


Rivs A., Oreno J. L. and Macarroén A. Equilibres liquide — vapeur de mélanges binaire donnant une réaction 
chimique : systémes méthanol — acide acétique ; éthanol — acide acétique ; n-propanel acide acétique ; n-butanol — acide 
acétique., Chem. Engng. Sci, 1959 10 105. 

On p. 110° Les Tableaux numériques 3 a 7” should read Les Tableaux numériques 3 6." 

On p. 111 in the Notation “a — volatilité relative = Y,. X, X,. 7," should read “ « 
Y,. Ny Vy.” 

Tables 1 — 6 were omitted from the text and are given below, 


Tableau 1. Constantes physiques des produits utilisés 


volatilité relative - 


Ethanol | n-Propanol n-Butanol Acide acétique 


Méthanol 


Index de 
refraction 


Observe 
Bibliographic 


(20°C) 13286 
(20 1.3288 
(5) 


(20,5 C) 1.3612 
(20,5 ©) 1.3610 
[3] 


(20°70) L856 
(20°C) 1.3854 
(5) 


(20-C) 1.3990 
(20°C) 1.3901 
(5) 


(22,9°C) 1.3714 
(22,9°C) 1.8715 
(5) 


Densite 

(g ‘ml) 
Observe 
Bibliographic 


0.8077 
O.8082 


[3] 


(20 ©) 0.7809 
(20 ©) O.7805 
(3) 


(188°C) 0.8044 
(188°C) 08047 (22. 
(9) 


(13-C) 0.7958 
(13 €) 0.79577 


[9] 


| 


Température 
débullition | 
(697 6mm) 115.1 

(6097, 6mm) 115,11 


[10] 


(710.6rmm) 95.3 115.8 C 
(710,6rmm) 95,32 (708, 4mm) 115,65 
(9) | (9) 


(710.600) 76,8 
(710.6mm) 76,78 


{9} 


(7038.3mm) 62.0°C 
(708. 62.85 


(9) 


Observe 
Bibliographie 


Température de 

congeélation 

16,58 — 16,64°C 
16.7 


{9] 


| 


| 
Systeme éthanol-eau (Figs, 3 et 4) 


| 
Xx, 


Tableau 2. 


Y2 


1,000 
0.998 
0,905 
1,008 
1015 
1,016 
1,035 
1,046 
1,204 
1,285 
1,341 
1,573 


0,088 
0.165 
0,282 
0,229 
0,390 
O419 
0,471 
0,572 
O.597 
0,618 
| 0,682 


0.008 
0.018 
0,026 
0.028 
0.070 
0,086 
0,105 
0,120 
0,291 
0,361 
0,407 
0,561 


| 
| | 
FH | VOL. 
- 10 
| | 1959 
} 
| 
| 
|_| || 
| | = 
1 709.1 95.9 7,7 5.35 
2 708.6 93.8 95.5 4.82 
3 709.1 M8 03,6 5.02 
4 709,1 1.6 4,63 
| 709,1 864 88,2 3.82 
6 709,1 85.6 874 3,68 
7 709,10 84.0 85,8 | 3,26 
8 709,1 83.4 85,2 | 3.01 
707.3 79.9 81,7 1,717 
10 709,1 81.2 1,473 
11 706.3 78.8 80,7 1,378 
12 706.3 | 77,7 79.6 1,154 
288 


Errata 


Tableau 3. Systeme méthanol—acide aceétique (Figs. 5, 6, 10 et 14) 
] 


P t 


710.5 115.8 0,000 0,000 — 

709.3 112.0 0,037 0,109 555 1,060 3.18 
709. 109.9 0,058 0.165 5 1,082 3.21 
709.3 105.0 0,107 0,303 1,113 3,62 
102.8 0.137 0.353 1,163 3.44 
709.3 O74 0,208 0491 1,195 3,67 
709.3 O48 0,250 0.556 0,705 1,205 3,75 
705.9 92,7 0,280 0,602 0.720 1,200 3.91 
705.9 91,7 0.306 0,630 0.726 1,210 386 
86,8 0,387 O73 0,798 1,180 4,39 
709.3 84,7 0427 0.776 0.817 1,150 4,65 
81,7 0.492 0,831 0,843 1,100 5,08 
705.0 79.9 0,853 O.875 1,075 5,42 
705.0 0,601 0.911 0,920 0,910 6.79 
703.3 78,7 0.657 0.934 0,940 0.858 7.39 
703.3 0.713 0.951 0,955 0,828 7.81 
703.3 69,1 O.785 0.971 0.960 0,725 9,17 
704.7 67,3 0.839 0.983 0,985 0,627 11,09 
65.9 O.882 0,989 0,996 0,571 12,02 
703.3 63.8 0.959 0,996 1,000 0,660 10,64 
703. 63,0 0.995 0,999 1,000 0,703 5.02 
703.3 62.9 1,000 10000 


Tableau 4. Systéme éthanol-acide acétique (Figs. 5, 7, 11 et 15) 


| | | 
P X, 
710.5 0,000 0,000 
708.1 0,055 2,06 
707.5 0.073 0.144 0,560 2,14 
O.105 O97 0,576 2.14 
700.8 107.7 0.133 0.274 0,664 2.46 
207.5 107.2 0,166 0.312 0,622 1,097 2,28 
708,1 104.7 0,207 0,393 0,686 1,106 2,48 
710.0 103.5 0.233 0,487 0,710 1,108 2.55 
689.7 99.9 0,282 0.526 0,764 1,100 2,82 
705,5 98,8 0.347 0,749 1,094 2,79 
706.0 93.0 0.460 0.750 0,870 1,000 3,53 
690,6 11 0.516 0.793 0,860 0,970 3.59 
707.2 88,0 0.587 O.854 0,903 0,886 4.11 
705.4 86,5 0,659 0,900 0,924 0,804 4,66 
697.9 83,8 0,728 0,934 0,947 0,730 5,29 
700.8 81.3 | 0,816 0.966 0966 0610 6,41 
17 706,0 78,3 0.924 0,990 0,996 0,492 8,14 


18 710.6 76.8 1,000 1,000 | 


Remarque. Les points entre parenthése () correspondent a des mesures realisées avec un autre appareil d'équilibre 
quelque peu différent, mais fonctionnant aussi selon le méme principe. Il n’avait qu'un seul réservoir, chargé avec un 
mélange binaire de produits & 99,5 pour cent. 
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Tableau 5. Systéme n-propanol-acide acétique (Figs. 5, 8, 12 et 16) 


| 


0.000 0,000 

0,126 0,159 1,023 

0.132 1,040 

0,205 0,256 1,023 

O212 0.397 1,025 

0,493 | 0,990 

0.437 OS812 0,950 

0.488 0,638 O.825 0.920 

107.2 0.500 0,651 0,920 

105.5 0,745 O.864 0.861 

103.1 0.669 O.831 0.935 O.775 

101.0 O.757 0,900 0.950 

0,940 0.975 

O84 O.868 0,960 0.973 

0.908 0,976 0.995 
95.5 1,000 1,000 — 


Tableau 6, Systéme n-butanol-acide acétique ( Figs. 5, 9, 13 et 17) 


P | t | XxX, | y, | 

115.8 0,000 - 

702.2 115.7 0.018 0,009 1,008 

702.2 116.4 0.060 0.039 O.625 1,000 

TOOT 116.9 0.105 0.078 0.706 0,989 

706.9 0.085 0.600 0,998 

O.161 O.124 O.718 0.984 

118.6 0.720 0.973 

706.7 118.8 0,200 0.968 

11a 0.250 O.741 0.000 

119.8 O.355 O.755 0.945 

706.8 120.3 0.459 0.439 0,907 

707 3 120.2 0.485 O.831 0,890 

7073 120.5 0,555 O.865 0.856 

707.9 120.1 0.590 O64 O.880 O.828 

7079 119.8 0,657 0,702 0,922 0,769 

707.9 194 O73) 0,790 0,946 0,700 

705.6 118.2 0,889 0.980 

706.4 0.920 0,960 

116.6 0.054 0.980 0,900 0.413 

115.8 1.000 1,000 

| 

Remarque. Pour construire la Fig. 9 et pour calculer les coordonnées du point azéotropique, nous avons corrigé le. 
températures 4 706 mm Hg (pression moyenne des déterminations), en employant le coéfficient At/ AP = 0,042°C /mm Hgs 
Les coordonnées du point azéotropique, calculées selon la méthode de Kaiser (13), sont: P = 706 mm Hg, tf = 120,3°C, 
X, = = 0518, 


| 
No. | P t 
1 710.5 
2 704.2 1,31 
3 706.8 | 1,20 
708.2 | 1,32 
706.8 1,45 
706.8 1,61 
701.8 1,76 
706.8 1,85 
| 1,86 
704.2 2,06 
701.8 2.43 
12 702.9 2.89 
13 702.9 3,37 
4 704.2 3,65 VOL. 
15 703,7 4,12 10 
16 710.6 | 1959 
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Heertyes P. M., van Mens M. H. and Buraye M. The absorption and chemisorption of 
ammonia in a jelly. Chem. Engng. Sci. 1959 10 47. 


The right hand side of equation (4) on page 48 should read : 


Cy* (Dy)! exp (— By) 
erfe (2/Bp)} 


IL. 
LO 
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